共查询到20条相似文献，搜索用时 187 毫秒
1.
2.
传统的Kirchhoff传播算子结构简洁,适用于描述横向均匀介质中波的传播.Ray-Kirchhoff传播算子较为精确地描述了波在非均匀介质中传播的运动学特征,其理论上的先天不足依赖于介质的复杂性.本文通过Born序列逼近波在非均匀介质中传播的大角度波分量,提出一种Born-Kirchhoff传播算子,将传统Kirchhoff传播算子的适用范围扩展至非均匀介质,同时描述波的运动学和动力学特征,其精度取决于Born序列逼近的阶数.利用Born序列频散方程,可以精确分析各阶Born-Kirchhoff传播算子对波长、传播角和非均质性的尺度依赖特征,其中,一阶Born-Kirchhoff传播算子的精度高于传统的相屏传播算子.波数域的Born-Kirchhoff传播算子对于高波数波是奇异的,导致波数域数值计算发散,但其空间域版本是非奇异的,无条件数值稳定,可通过Kirchhoff求和数值实施.本文给出各阶Born-Kirchhoff传播算子及其频散方程,可用于不同程度非均匀介质中的波传播模拟,复杂构造地震成像和速度估计.本文利用零阶和一阶Born-Kirchhoff传播算子计算简单二维模型的合成地震图,并与边界元法进行了比较. 相似文献
3.
映射动力无穷元及其特性研究 总被引：6，自引：0，他引：6
本文提出了一种映射动力无穷元,这种单元不仅可以表征在均匀介质中传播的不同波数的波,而且可以模拟地基中不同介质内的波的传播特性。文中通过定义无穷元在局部坐标系中的波数(名义波数)β和在整体坐标系中的波数(介质中的实际波数)β_G,使得:1.整体坐标系中不同波数的无穷元均可用在局部坐标系中具有同一名义波数β的无穷元来模拟。这一模拟是通过改变无穷元边中点的整体坐标以及引入适当的过渡单元来实现的。2.适当地改变有限元与无穷元连接边界的法线方向与无穷元各节点在整体坐标系中的相对位置,还可以改变波在整体坐标系中的传播方向,从而为模拟不同方向传播的波提供了一种简便的技术。文中最后给出了无穷元在整体坐标系中常用的几种单元形态及其传波特性,并用算例对其精度进行了验证。 相似文献
4.
在海洋可控源电磁法勘探中,接收站常置于海底.在进行海洋电磁场模拟时,由于海水和海底介质存在显著电性差异,这给海底接收点处场值的求取带来困难.本文提出一种新的接收点插值算法,该算法考虑到海底电场法向分量不连续性问题,用法向电流分量进行插值以准确求取海底任意接收点处电磁场值.本文利用交错网格有限差分法实现了二维介质中频率域海洋可控源法(CSEM)正演.对构造走向做傅里叶变换,将三维电磁模拟问题转换为波数域2.5维问题,即三维场源激励下针对二维地电模型的电磁模拟问题.使用交错网格有限差分法,基于一次场/二次场分离方法导出波数域二次电场离散形式,并进一步求得波数域电磁场.采用本文提出的改进的插值算法可求得海底任意接收点处波数域电磁场,采用傅里叶逆变换对波数域电磁场进行积分可得到接收点处空间域电磁场.模型算例表明,与常规的线性插值和严格插值算法相比,本文提出的改进的插值算法具有更高的精度. 相似文献
5.
常见弹性波动理论的建立是基于介质均匀这一基本假设,实际介质的非均匀性非常普遍.为研究连续介质中波的传播特征,本文从弹性力学中建立弹性波动方程的三个基本方程出发,考虑连续介质弹性参数的空变特征,建立非均匀介质的弹性波动方程,利用Alkhalifah声学近似思想建立位移表征的纵波波动方程,利用本征值问题求解方法建立标量波频率-波数域传播算子,从而建立描述纵波传播的标量波方程,其中波函数为纵波位移的散度,不同于均匀介质标量波方程的波函数为位移势.随后推导含PML边界波动方程差分格式并建立不同模型数值模拟进行数值试算,与均匀假设标量波方程和变密度方程对比证明本方法的准确性和稳定性. 相似文献
6.
现有人工源频率域电磁法的研究大多仅针对某种具体的方法,而较少将问题综合起来分析.本文综合多种方法的共同点提出了层矩阵法,它采取了源置于层间的模型进行公式的推导,理论上可以计算任意层状介质中任意位置的任意场源在空间中任意位置产生的场强,可适用于多种电磁法的正演模拟计算.层矩阵法的核心是对空间域的变量x, y, z中的x和y变量进行傅氏变换后转换到波数域kx和ky中,在波数域利用边界条件,用层矩阵建立起各层的关系后计算得到各层的波数域电磁场值,然后经过二维反傅氏变换最终得到空间域中任意位置的场值.因为文中定义的层矩阵是建立层关系的关键,所以称此方法为层矩阵法.本文以水平电偶源为例独立推导了层状介质中人工源频率域电磁场解的理论公式.为了验证方法的正确性,文中建立了多种模型,利用自行编排的程序将层矩阵法与现有文献的各种解析公式的解进行了对比,结果表明本文提出的层矩阵法是灵活的、可靠的. 相似文献
7.
提出了位场向下延拓的波数域迭代法. 对水平面上的位场观测值进行Fourier变换,得到其波谱. 根据第一类Fredholm积分方程的空间域迭代解法,推导出计算向下延拓水平面上位场波谱的波数域迭代公式. 在波数域中进行迭代,一直进行到相继两次迭代近似解的差值最大绝对值小于给定的精度,或迭代达到给定的最大迭代次数. 对这种迭代近似解进行Fourier逆变换,得到向下延拓的位场. 数值计算结果表明：与空间域迭代法比较,这种波数域迭代法简单、快速,并有同样好的向下延拓效果. 本文还证明了这种迭代法是收敛的,并给出了它的收敛特性和滤波特性. 相似文献
8.
9.
单程波近似实际上是一种多次前向散射和单次后向散射近似.利用单程波近似来描述波传播可以极大地节省地震数值模拟的计算时间和内存,实现地震波长距离传播模拟和三维地震模拟快速计算.本文基于单程波近似和波动积分方程的分离变量逼近,从广义Lippmann-Schwinger波动积分方程推导出耦合反射/透射系数的单程波传播算子.该算子由两部分构成：分离变量Fourier单程波传播算子和薄板间的反射/透射系数表达.前者将常规的Fourier分裂步单程波传播算子(SSF)推广适应横向强速度变化介质和大角度传播波场.后者是利用垂直波数来表示反射/透射系数,自然耦合到波场传播的计算过程中,其为地质界面倾角的隐式表达,精确描述振幅随入射角的变化,能适应任意复杂的模型.通过两个数值算例和一个实际地质模型的计算,本文将该方法和边界元法进行了比较,结果表明：在算例给出的介质横向速度变化情况下,本文提出的方法在相位和振幅方面与全波数值方法基本吻合. 相似文献
10.
《地震地质》2020,(3)
完全匹配层吸收边界条件通常可以很好地吸收模型边界的地震反射波,但对于横向各向同性介质的模拟效果欠佳,且尚在发展阶段。为此,文中推导了横向各向同性介质中弹性动力学波动方程,给出了施加完全匹配层(PML)吸收边界条件下时间域二阶、空间域十阶精度的高阶交错网格的有限差分形式,并分别建立了均匀的垂直向对称轴的横向各向同性介质(VTI介质)和倾斜向对称轴的横向各向同性介质(TTI介质)模型。计算结果表明,对于对称轴为任意角度的横向各向同性介质,当PML边界层厚度达到一定的数值时,可以很好地抑制人工边界所产生的地震波反射效应,且PML的吸收效果不会被入射角与入射波频率影响。 相似文献
11.
声波测井中的相速度与群速度讨论(英文) 总被引：3，自引：2，他引：1
声波测井过程中获取的速度到底是相速度还是群速度,目前仍存在一些争议,本文从理论分析和数值模拟的角度,使用三种模型对这一问题进行了研究。首先,构造一个相速度与群速度可调的稳态声波传播模型——不同声速的两个平面波叠加模型,利用慢度时间相关(STC)方法提取声波波速,数值模拟结果表明,无论相速度较大或是群速度较大,STC方法提取出来的波速都是相速度;其次,通过频散分析和割线积分得到刚性壁圆柱流体模型中的频散曲线与分波波形,使用STC方法得到的速度与相速度的频散曲线吻合较好,而直接读取波至获得的速度与群速度的频散曲线趋势一致;最后,利用频散分析和实轴积分方法,获得偶极子在慢地层中激发的模式及全波波形,得到的结果再次验证了刚性壁圆柱流体模型中的结论。 相似文献
12.
本文采用三维应力-速度有限差分(SV-FD)方法,数值模拟了横向各向同性(TI)地层对称主轴与井轴斜交情况下正交偶极子声源激发的井孔声场.主要解决了与倾斜角有关的三维空间弹性模量矩阵的推导,柱坐标系下应力-速度有限差分方程组的建立,井轴上场点奇异性与内边界处理等几个关键问题,提高了计算精度.在横向各向同性地层对称主轴与井轴平行的情况下,与实轴积分法所得结果进行了对比,验证了本文方法的正确性.计算了不同倾角情况下xx和yy两分量的弯曲波,并用频域加权相似法提取了弯曲波频散曲线,结果显示了横向各向同性介质中不同方向偏振的弯曲波传播的分裂现象,其频散曲线在低频段分裂,随频率增大而逐渐重合.弯曲波低频截止频率处速度与理论公式得到的横波速度基本符合. 相似文献
13.
14.
A wideband dipole signal is required for dipole dispersion correction and near-borehole imaging. To obtain the broadband flexural wave dispersion, we use a nonlinear frequency modulation (NLFM) signal and propose a segment linear frequency modulation (SLFM) signal as the dipole excitation signal to compensate for the excitation intensity. The signal-to-noise ratio (SNR) of the signal over the entire frequency band is increased. The finite-difference method is used to simulate the responses from a Ricker wavelet, a linear frequency modulation (LFM) signal, an NLFM signal, and an SLFM signal in two borehole models of a homogeneously hard formation and a radially stratified formation. The dispersion and radial tomography results at low SNR of the sound source signals are compared. Numerical modeling suggests that the energy of the flexural waves excited by the Ricker wavelet source is concentrated near the Airy phase. In this case, the dispersion is incomplete and information regarding the formation near or far from the borehole cannot be obtained. The LFM signal yields dispersion information near the Airy phase and the high-frequency range but not in the low-frequency range. Moreover, the information regarding the formation far from the borehole is not accurate. The NLFM signal extends the frequency range of the flexural waves by compensating for the excitation intensity and yields information regarding the formation information, but it is not easy to obtain. The SLFM signal yields the same results as the NLFM signal and is easier to implement. Consequently, the dipole detection range expands and the S-wave velocity calculation accuracy improves. 相似文献
15.
研究声波在倾斜充液井孔中的传播对于声波测井数据处理和解释具有重要意义.应用三维交错网格有限差分方法模拟了处于倾斜各向同性分层地层中的井孔声场.首先,针对均匀地层中单极子声源在裸眼井中激发的声场,将有限差分的结果和实轴积分法的结果进行对比验证.然后,采用单极子和偶极子两种声源,针对地层分界面和井轴间的不同倾角,计算了相应的声场分布和井轴上的接收波形.数值计算的结果表明,当声源处于倾斜分界面以下,即处于快速（下方）地层,接收器处于倾斜分界面以上（慢速）地层时,随着地层倾斜角度的加大,测得的慢度值从接近上方慢速地层值逐渐减小直至接近下方快速地层的值.任何源距情况下测得的首波慢度均小于上方地层实际的纵波慢度.并且,慢度与源距的关系曲线随源距的加大逐渐平缓.用偶极子声源激发得到的横波慢度与纵波结果相同,并表现得比纵波对倾角的改变更敏感.上述结论在本文中用声场快照和利用合成接收波列的慢度计算得以清楚显示,并且用射线声学理论验证. 相似文献
16.
在随钻测井条件下,由于钻铤占据了井孔内的大部分空间,充液井孔中沿着井轴方向传播的模式波的特性与电缆测井非常不同.本文建立了随钻正交偶极子测井声学模型,采用三维有限差分方法模拟了偶极子声源在随钻条件下各向异性地层井孔内激发的声场,研究了地层的声学各向异性在随钻正交偶极子声波测井中的响应特征.数值模拟结果表明,在随钻测井条件下,对于井轴同TI地层对称轴垂直的情况,弯曲波分裂现象仍然存在,通过正交偶极子测量方式和合适的反演算法能够准确有效地确定地层的快横波方位角,可以考虑采用同正演理论相结合的反演算法来获得地层的快、慢横波速度及声学各向异性信息;对于井轴同介质对称轴呈一定夹角的TI地层井孔,情况变得非常复杂,不同井斜倾角下弯曲波的速度的变化趋势并非同对应的地层横波速度的变化趋势完全一致,不过在一定的频段内,地层横波速度仍然是弯曲波的最主要控制因素.对于本文研究的模型,当井轴同介质对称轴的夹角大于大于60°时,此时获得的弯曲波的各向异性值基本能够反映对应角度下地层横波速度的各向异性信息. 相似文献
17.
偶极横波远探测技术作为近年来发展起来的一门地球物理测井技术,一个重要应用前景就是在裸眼井和套管井测井时寻找井外地层中的油气构造.本文从井中偶极声场积分表达式的远场渐近解出发,对裸眼井及套管井中偶极声源的辐射场进行了对比分析,首次得到了偶极声源在不同声源频率和不同地层类型下的辐射特征.并深入研究了套管井中偶极声源的SH横波远场辐射特性,给出了不同耦合情况下的SH横波远场辐射特性的变化规律,对比了套前套后SH横波远场辐射的异同.结果表明：井中偶极声源的远场辐射特性受声源频率和地层类型影响较大;在套管井中,偶极声源辐射到地层中的横波能量也可以进行远探测,尽管与裸眼井情况相比,辐射能量会有所降低.本文的结果为偶极横波远探测技术的应用提供了理论支持. 相似文献
18.
地-井瞬变电磁(TEM)方法是在地表发射,探头沿钻孔(井)测量瞬变响应的一种井中物探方法.由于接收探头沿钻孔深人地下,因此能够获得更加可靠的地下目标体信息,尤其当存在低阻覆盖层、浅部矿化等地质干扰,或者勘查深部规模不大之良导矿体时,地-井TEM方法的优势更加明显.相对于磁性源系统,电性源系统发射功率大,辐射面积广,更适合地形复杂地区之深部找矿.本文的目的是研究电性源地-井TEM的近似解释方法,首先给出了层状介质电偶极子在地下的TEM响应,进而通过电偶极子叠加的方式获得了电性源条件下的TEM响应.利用感应磁场与均匀半空间电阻率的单调关系通过反函数定理进行了全域视电阻率定义.理论模型的视电阻率计算结果显示,不同深度测点的视电阻率曲线首支不同,而尾支基本一致.这说明了全域视电阻率是测点周围有限范围内介质的综合反映,所以不同深度测点的视电阻率曲线首支所反映的范围是不同的.绘制了不同时间道视电阻率剖面曲线,以期显示地下电性分布规律.计算结果表明全域视电阻率定义能够基本可靠地反映地下信息,从而为该技术在矿区深部找矿的应用提供有力支持. 相似文献
19.
偶极声源的辐射特性及仪器的四分量正交采集系统使得偶极横波远探测可在一维井筒中对井外构造进行方位探测.基于反射横波幅度随声源偏振方向呈余弦平方变化这一特征,本文提出一种利用函数拟合快速反演井周构造走向的方法.根据这一规律,该方法仅用4个方位剖面数据即可准确反演井周构造走向,同时得到井周构造的离井距离和倾角.数值模拟和物理实验模拟均验证了幅度随方位的变化以及反演方法的可行性.根据这一方法做的现场实例也进一步证明了该方法的适用性.本文工作为利用四分量偶极数据进行井周构造的反射横波成像提供了实用而有效的处理方法. 相似文献
20.
现有的随钻声波测井通常采用三种测量方式:单极子、偶极子及四极子测量模式.分别采用单极子源、偶极子源和四极子源激发单极子、偶极子和四极子信号对地层速度信息进行测量.由于随钻测井环境比较复杂,钻杆在井孔中占据大部分空间,钻杆速度高于地层速度,导致接收到的信号中仪器波能量占主导,无法有效地提取地层的纵横波速度信息.本文探讨一种采用非对称偏心的点声源,一次激发,利用方位分布的四个接收器接收,经过合成后同时分别得到单极子、偶极子和四极子模式波信号的测量方法,为建立同时利用多极模式波信号测量地层纵横波速度的随钻声波测井实用技术提供基础.本文利用1:12的缩比模型在实验室对该方法进行了考察,实验结果证明,可以通过非对称点声源得到有效的地层单极子、偶极子和四极子模式波信号.通过分析和对比实测数据和理论模拟的频散曲线,本文进一步提出利用三种模式波的频散曲线联合反演确定地层横波速度的方法. 相似文献