地震勘探——相关概念(一)

news2024/7/6 17:45:24
地震波的基本介绍

波前:波在同一时刻所到达的点所构成的面,这个面上构成的相位是相同的。波前的形状取决于传播介质的物理性质。我们可以用地震波动方程模拟波前变化(波场快照)。
在这里插入图片描述
射线(Ray):是一条假想的“线”,垂直于波前。射线指明了波的传播方向,但不能反映地震波的能量变化。射线的形状取决于地下岩层的性质。
均匀介质
在这里插入图片描述
在这里插入图片描述


惠更斯原理和费马原理

地震波在传播过程中遵循惠更斯原理费马原理
惠更斯原理:对于任何一种波,其波前上的任何一点都可以作为子波的波源,各个子波波源波前的包络面就是下一个波的波前。(一个波阵面的每个点可各看做是一个产生球面子波的次级球面波的中心波源,次级波源的波速与频率等于初级波的波速与频率;而且,以后任何时刻波阵面的位置是所有这种子波的包络面。这个结果实质上是“光学平行”面的一个作图法则,所以有时叫做惠更斯作图法。)
包络面:在波的传播过程中,总可以找到同相位各点的几何位置,这些点的轨迹是一个等相位面,叫做波面(即包络面)。
在这里插入图片描述
可以利用惠更斯原理证明折射波的传播路径
在这里插入图片描述

费马原理:波所走的两点之间的实际路径是传播时间最短的那一条。
在这里插入图片描述
可以用费马原理证明反射定律和投射定律
在这里插入图片描述


地震波的主要参数

描述地震波的主要参数:振幅、频率、相位、周期、波长、波数、时间、速度、波峰、伯父、主频、频带宽度等
振动曲线:在某一个固定的接收点,质点位移随时间变化的曲线。(在时间域里研究地震信号的特征)
在这里插入图片描述
波剖面:在同一时刻质点唯一随接收点位置变化的曲线。(在空间域里研究地震信号的特征)
在这里插入图片描述
描述振动曲线的参数

  • 振幅(amplitude) A:在振动曲线上,质点唯一离开平衡位置的距离。
  • 波峰和波谷:振动曲线上,正的极大值称为波峰,负的极大值称为波谷
  • 周期(period) T T T:在振动曲线上相邻极大值间的间隔,单位为秒或毫秒。
  • 频率(frequency) f f f:单位时间内完成周期性变化的次数,是描述周期运动频繁程度的量。
    • 频率与周期互为倒数。: T = 1 f T = \frac{1}{f} T=f1 f = 1 T f= \frac{1}{T} f=T1
  • 角频率 ω \omega ω ω = 2 π f \omega = 2\pi f ω=2πf
  • 相位(phase):对于一个波,特定的时刻在它循环中的位置,单位是度,最大值是360度。
  • 波速(velocicy) v v v:单位时间内振动所传播的距离,单位是m/s或km/s。
    在这里插入图片描述
    描述波剖面的参数
  • 波长(wave lenth) λ \lambda λ:波在一个振动周期内传播的距离,单位是米。
  • 波数(wave number) k k k:空间频率,指的是单位距离内的波长数。
    • 波长与波数互为倒数: λ = 1 k \lambda= \frac{1}{k} λ=k1 k = 1 λ k= \frac{1}{\lambda} k=λ1
      在这里插入图片描述
      波速 v v v,频率 f f f,周期 T T T,波长 λ \lambda λ,波数 k k k的关系:
  • λ = v T = v / f \lambda = vT = v / f λ=vT=v/f
  • v = λ / T = λ f = f / k v= \lambda / T= \lambda f = f / k v=λ/T=λf=f/k

像地震波这样的非周期信号,可以利用傅里叶变换提取波的频率、相位等特征参数。
傅里叶变换基本思想:几个不同振幅、相位、频率的信号叠加可以得到一个复杂的非周期信号。
在这里插入图片描述
若查看非周期的复杂信号的不同频率、相位、振幅等信息组成,叫做信号的分解。
地震波的频谱的两个主要参数:

  • 主频 f 0 f_0 f0:频谱曲线极大值对应的频率。
  • 频带宽度:频谱曲线上地震波绝大部分能量集中在那个频率范围内。
    在这里插入图片描述
    不同地震信号的频带范围具有较大差异:
    在这里插入图片描述

地震波速度

地震波速度:地震波在地下岩层中传播的距离与时间的比值,单位是m/s或km/s。地震波速度在空气中的传播速度是340m/s;近地表速度大约是数百米/秒;地震波的纵波速度在地下2885千米是达到最大,可以达到13.7km/s。
影像地震波速度的主要因素:
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
(面波的速度随频率变化而发生变化,有频散现象)
墨霍界面: 莫霍面是划分地壳与地幔的界面,位于地表以下数千米到30—40km的深度,纵波速度从7.0 km/s跳跃增加到8.1 km/s左右,横波速度从4.2 km/s增加到4.4 km/s左右。
古登堡面: 古登堡面是具有高密度的固体地幔与具有液态性质的外核之间的界面(核幔界面)。位于地下2900km深度,横波到达这一界面即告消失;纵波能够通过,但通过后其速度明显减慢。同时,纵波到达该界面还明显地发生反射与折射,出现纵波阴影带,即纵波在地表的一定区间不能被接受的地带。
在这里插入图片描述


地震子波 seimic wavelet

地震子波是组成地震记录的基本元素。
现代数学中,wavelet被称为子波。地震勘探会用到wavelet transformer(小波变换)这个数学工具,但是地震子波和小波没有任何关系。
地震子波的形成:炸药产生一个延续时间极短的尖脉冲,爆炸点附近的介质以冲击波的形式传播,当爆炸脉冲向外传播一定距离以后,地层产生的弹性形变再向外传播,由于介质对高频成分的吸收,波形发生明显的变化,直到传播更大的距离以后,波形逐渐稳定,形成一个具有两到三个相位、有一定延续时间的地震波,我们称这个波为地震子波。
在这里插入图片描述
地震子波的特点:有确定的起始时间(因果性);有限的能量;有一定延续时间的信号。
在这里插入图片描述
在这里插入图片描述
地震子波的关键属性:振幅、频率、相位
通过傅里叶变化对地震子波进行频谱分析,我们从而可以了解地震子波的频率特征和相位特征。
不同时间延续的地震子波对应的频谱中,我们可以看到子波越尖锐,则频带越宽,地震勘探分辨率也就越高,如下图所示:
在这里插入图片描述
根据地震子波的相位特征的不同,可以将子波分为以下四类:
在这里插入图片描述
不同类型地震源激发的地震子波特征不同,如图:
在这里插入图片描述理论地震子波
雷克子波零相位子波在地震正演模拟、合成地震记录制作、反演中常用到:
在这里插入图片描述
Klauder零相位子波在地震可控震源研究中适用:
在这里插入图片描述
Ormsby零相位子波
在这里插入图片描述
Ormsby最小相位子波
在这里插入图片描述
随时间发生变化的子波叫时变子波;随空间发生变化的子波叫空变子波
地震子波的极性
在这里插入图片描述
地震子波的作用
在这里插入图片描述
在这里插入图片描述


地震合成记录

相关介绍:地震褶积模型:
在这里插入图片描述
反射褶积模型认为——地震记录 g ( t ) g(t) g(t)是真元地震子波 w ( t ) w(t) w(t)和地下地层反射系数 e ( t ) e(t) e(t)褶积的结果:

褶积模型的离散形式:
g k = ∑ i = 0 m e i w k − i g_k = \sum_{i = 0}^me_iw_{k-i} gk=i=0meiwki
在这里插入图片描述
制作合成地震记录:对某一口井进行声波测井、密度测井,可以得到声波测井曲线和密度测井曲线,将速度和密度对应相乘,得到波阻抗曲线,然后根据地震垂直反射系数方程计算反射系数序列;给定一个震源子波,地震记录上看到的反射波波形就是地震子波在地下各反射截面上发生时反射形成的。
在这里插入图片描述
在这里插入图片描述
不同频率的地震子波与反射系数褶积,得到的合成地震记录是有差异的——地震子波主频越高,地震纵向分辨率越高。如下图所示:
在这里插入图片描述
在这里插入图片描述
实际的地震记录是地震子波、反射系数以及各种影像因素共同褶积的结果:
在这里插入图片描述
人工合成地震记录制作中,井的纵坐标是深度,合成记录的纵坐标是时间,所以我们要用合适的平均速度进行时深转换:
在这里插入图片描述
地震合成记录可以对层位进行准确的标定:
在这里插入图片描述
地震反演地震正演
在这里插入图片描述

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/1069638.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

java API 包装类Integer,int和String互相转换

integer类 package daysreplace;import com.sun.jdi.IntegerValue;import java.util.Arrays;public class Test {public static void main(String[] args) { //Integer a new Integer(29);//过时Integer b new Integer("30");//过时System.out.println(a);System.o…

matplotlib制图初级篇

做任何的报表分析,最后是达到可视化的目的。 全部都是一堆数字,那肯定不是一个合格的数据报表制作分析人员: 需求:根据excel表,生成折线图、柱状图和饼状图 1、pandas读取数据 说明:本机的运行环境为mac…

碰撞检测 Neon优化

1、碰撞检测 碰撞检测的原理,两个圆心距离d <= r1 + r2即发生碰撞, 对于下图这种,d的计算方法即为余弦定理 2、Neon优化 这段代码过于简单直接使用neon intrinsic编写 #include <arm_neon.h> #include <stdio.h>struct circle {

雷电模拟器在打开“指针位置“后,无效,没有指针xy轴坐标显示?(解决方法)

解决方法&#xff1a;雷电模拟器我甘霖娘***(不是 1. 打开"指针位置" 2. 右击雷电模拟器图标 - 打开文件所在位置 - 找到vms文件目录并进入 3. 新建名为debug的txt文件 4. 重启雷电模拟器 5. 已解决

Day08-面向对象

1. 类和对象 1.1 类和对象的理解 客观存在的事物皆为对象 &#xff0c;所以我们也常常说万物皆对象。 类 类的理解 类是对现实生活中一类具有共同属性和行为的事物的抽象 类是对象的数据类型&#xff0c;类是具有相同属性和行为的一组对象的集合 简单理解&#xff1a;类就是…

集睿致远CS5518国产MIPI转LVDS 点屏方案芯片可pin to pin替代国腾GM8775C

CS5518是一款MIPI DSI输入与1或2Port LVDS输出转换芯片。Pin to Pin替换GM8775C&#xff01;MIPI DSI最多支持4Lane&#xff0c;1Lane最大运行速率为1Gbps。LVDS支持18或24位像素&#xff0c;25MHz至154MHz&#xff0c;支持VESA或JEIDA格式。单路1.8V供电方式&#xff0c;可选配…

使用VisualStudio制作上位机(补充)

使用VisualStudio制作上位机(补充) Author:YAL 这一部分,主要是补充CAN发送的实现。大概的思路是在主界面上,通过点击“VIN读写”然后弹出发送界面,用户在该界面填上信息后,点击就发送一次。当然这仅仅是一个例程,大家可以根据自己的需求来做。 先新建一个窗体,用于填…

诺贝尔奖(Nobel Prize)2023——生理/医学奖

诺贝尔奖&#x1f3c6;&#xff08;Nobel Prize&#xff09;2023——生理/医学奖 2023年10月2日颁布了诺贝尔生理学/医学奖&#xff08;Nobel Prize in Physiology or Medicine&#xff09;&#xff0c;Katalin Karik&#xff08;卡塔琳考里科&#xff09;和Drew Weissman&…

如何实现 Es 全文检索、高亮文本略缩处理

如何实现 Es 全文检索、高亮文本略缩处理 前言技术选型JAVA 常用语法说明全文检索开发高亮开发Es Map 转对象使用核心代码 Trans 接口&#xff08;支持父类属性的复杂映射&#xff09;Trans 接口的不足真实项目落地效果 前言 最近手上在做 Es 全文检索的需求&#xff0c;类似于…

postman使用旧版本报错version mismatch detected后如何恢复使用

postman下载旧版本 目前作者使用10.6.0版本可以忽略登陆&#xff0c;所以可以下载此版本。 但是安装打开后会提示version mismatch detected&#xff0c;所以需要解决该办法&#xff1b; 修改文件名称 修改 C:\Users\XXX\AppData\Roaming\Postman 为 C:\Users\XXX\AppData\R…

Python接口自动化测试实战(完整版)

接口自动化测试是指通过编写程序来模拟用户的行为&#xff0c;对接口进行自动化测试。Python是一种流行的编程语言&#xff0c;它在接口自动化测试中得到了广泛应用。下面详细介绍Python接口自动化测试实战。 1、接口自动化测试框架 在Python接口自动化测试中&#xff0c;我们…

Compose 编译器版本和Kotlin版本对应关系

使用了最新的kotlin版本&#xff0c;compose报错&#xff0c;不兼容&#xff0c;在这里记录一下版本对应关系 值得注意的是Compose Kotlin 编译器扩展 (androidx.compose.compiler) 未关联到 Compose 库版本。相反&#xff0c;它会关联到 Kotlin 编译器插件的版本&#xff0c;…

软件测试银行项目网上支付接口调用测试实例

公司最近有一个网站商城项目要开始开发了&#xff0c;这几天老板和几个同事一起开着需求会议&#xff0c; 讨论了接下来的业务规划和需求策略&#xff0c;等技术需求一下来还要讨论技术需求&#xff0c; 确认后再慢慢的进入开发阶段&#xff0c;趁着闲暇时间新造的人想总结一…

每日Leecode算法题:1337.矩阵中战斗力最弱的k行

方法1&#xff1a;暴力破解 class Solution:def kWeakestRows(self, mat: List[List[int]], k: int) -> List[int]:return sorted([i for i in range(len(mat))], keylambda x:sum(mat[x]))[:k] 方法2&#xff1a;二分查找排序 class Solution:def kWeakestRows(self, ma…

Windows配置ADB工具

一、目的 在进行嵌入式开发时&#xff0c;我们经常使用ADB工具登录到开发板上进行命令操作&#xff0c;本篇我们介绍如何在windows平台配置ADB环境。 二、实战 1.下载adb工具包​​​​​​​https://developer.android.com/studio/releases/platform-tools?hlzh-cnhttps://d…

【Git】Git下载安装环境配置 下载速度慢的解决方案

这里写自定义目录标题 介绍一、下载官网下载镜像站 二、安装安装成功 三、Git三种界面介绍Git cmd界面展示git bash界面展示git GUI界面展示 四、环境配置配置流程1、打开环境变量界面2、添加环境变量 /删除环境变量3、在变量中找到Git\cmd的值就表示配置成功4、没有找到点击新…

在硅云上主机搭建wordpress并使用Astra主题和avada主题

目录 前言 准备 操作 DNS解析域名 云主机绑定域名 安装wordpress网站程序 网站内Astra主题设计操作 安装主题 网站内avada主题安装 上传插件 上传主题 选择网站主题 前言 一开始以为云虚拟主机和云服务器是一个东西&#xff0c;只不过前者是虚拟的后者是不是虚拟的…

编程每日一练(多语言实现)基础篇:阳阳买苹果

文章目录 一、实例描述二、技术要点三、代码实现3.1 C 语言实现3.2 Python 语言实现3.3 Java 语言实现3.4 JavaScript 语言实现3.5 Go 语言实现 一、实例描述 阳阳买苹果&#xff0c;每个苹果0.8元&#xff0c;阳阳第一天买2个苹果&#xff0c;第二天开始每天买前一天的2倍&am…

双端队列--二叉树 Z 字层序遍历

力扣103题----二叉树的锯齿形层序遍历 给你二叉树的根节点 root &#xff0c;返回其节点值的 锯齿形层序遍历 。&#xff08;即先从左往右&#xff0c;再从右往左进行下一层遍历&#xff0c;以此类推&#xff0c;层与层之间交替进行&#xff09;。 代码&#xff1a; public L…

遥感数据与作物模型同化:遥感数据、PROSAIL模型、DSSAT模型、参数敏感性分析、数据同化算法、模型耦合

查看原文>>>遥感数据与作物模型同化实践技术应用 基于过程的作物生长模拟模型DSSAT是现代农业系统研究的有力工具&#xff0c;可以定量描述作物生长发育和产量形成过程及其与气候因子、土壤环境、品种类型和技术措施之间的关系&#xff0c;为不同条件下作物生长发育及…