共查询到19条相似文献，搜索用时 171 毫秒
1.
在利用实际地震资料反演地球内部结构时,群速度往往比相速度更有意义.本文在水平层状均匀介质模型中,用快速广义反射透射系数方法（Pei,2008）计算了基阶瑞雷面波的群速度,并推导了Chen(1993) 4×4 层矩阵E的逆的解析表达式.文中用数值实验的方法合成了两个典型模型的群速度频散曲线,并与复合矩阵方法计算的结果进行了对比,结果表明本文的方法是切实可行的,与广义反射透射系数方法相比明显地提高了计算效率,尤其是在高频、层厚度大和层数多的情况下表现得更为突出. 相似文献
2.
层状介质的声波波动方程反演 总被引：4，自引：3，他引：1
基于广义反射透射系数矩阵正演方法 ,讨论了层状介质的声波波动方程反问题 .推导出波数频率域中的雅可比矩阵的解析表达式 ,其计算在正演过程中求出 .采用最小二方法可得到层介质参数 .数值结果表明反演方法的正确有效性 . 相似文献
3.
4.
为认识地震波诱导的电磁场的特性,本文研究地震波在孔隙介质中由于动电效应引起的电磁场.基于Pride弹性-电磁耦合方程组推导了双力偶震源对应的位移-应力-电磁场间断向量的表达式,模拟了双力偶源激发的震电波场.作为比较,还模拟了爆炸点源激发的震电波场.结果表明：存在伴随纵波的电场,其各分量的波形与固相位移对应分量的波形相似,但相位相反;存在伴随横波的磁场,其波形与固相位移波形相似;双力偶震源还激发出了独立传播的辐射电磁场,其速度比纵波至少高一个数量级,几乎是“瞬间同时”到达了每个接收器,但是其强度比伴随电磁场小得多,且随着源距增大而迅速减小.本文研究还表明：伴随纵波的电场强度不仅与地震纵波幅度和动电耦合系数有关,还与由介质孔隙结构决定的流-固两相动力协调性有关,存在一种动力协调介质,纵波在这种介质中不引起电场. 相似文献
5.
孔隙介质弹性波传播理论在地球物理勘探、地震工程和岩土动力学等领域有着广泛的应用.而孔隙介质中的弹性波受孔隙度、渗透率、流体黏滞系数等参数的影响,因此研究波场的传播特征将有助于分析和提取这些信息.本文在Biot理论的基础上,针对三维层状孔隙介质模型,利用在合成理论地震图的研究中已经被证实具有稳定、高效且适用范围较广的Luco-Apsel-Chen（LAC）广义反透射方法,给出了弹性波场的一种积分形式的半解析解,可通过数值方法高效、准确地计算层状孔隙介质中的理论波场,所以该积分形式的半解析解可为三维层状孔隙介质波场传播特征的理论数值模拟研究提供一种新的途径和手段. 相似文献
6.
为了更好地理解和解释土壤结构、多相流混合物和流体矿化度对震电测量的影响,本文基于Pride和Revil分别提出的动电耦合控制方程组,采用LAC广义反射与透射系数法对均匀、横向各向同性的多层含流体孔隙介质进行数值模拟研究.本文分别研究了包气带、不同土壤结构的地层和污染物侵入含水层等三种不同地质背景的震电信号.数值结果表明：含水饱和度剧烈变化比均匀变化能激发出更强的动电效应;按照含水饱和度划分的地层层厚与主导的地震波长一致时,震电信号振幅最大;由于S波比P波波长短,所以S波转换的磁场分辨“薄层”能力更强;Pride与Revil的理论模型在地层富含黏土且孔隙中流体部分饱和时存在差异;当土壤为三孔的形式且富含黏土矿物时,基本上不满足“薄”双电层假设条件;黏土含量越高,砂土含量越低,电场总场及界面响应振幅越大;和依照饱和度划分的地层对比,震电信号不仅可以探测到含非水相液体(NAPL)地层的存在,而且分辨率远高于前者;非水相液体的密度和黏度会对震电信号产生明显影响;震电信号同样对流体矿化度变化敏感,非同震电磁信号对流体矿化度变化敏感主要体现于矿化度变化界面附近的隐失波,自由表面附近接收的震电信号无... 相似文献
7.
基于前一篇文章中得到的关于三维层状孔隙介质中弹性波场的积分形式半解析解,本文通过离散波数法开展了数值模拟.将全空间均匀孔隙介质中单力点源和爆炸点源作用下弹性波场的解析解和我们的数值模拟结果进行对比,发现两者是完全一致的.而在一个两层半空间模型下的数值模拟,验证了固相位移Green函数的9组空间互易性情况.通过以上两种对比检验,验证了半解析解理论公式、数值模拟方法以及相应程序代码的正确性和可靠性.随后利用敏感度分析研究了不同的介质参数变化对爆炸点源在界面上会产生的反射波场的影响.通过垂直地震剖面模型的数值模拟,发现弹性波场能很好地反映孔隙介质物理性质的变化,同时也讨论了动力协调这一孔隙介质中的特殊现象.我们发展的基于半解析解的数值模拟方法可以为三维层状孔隙介质中弹性波传播特征的研究提供一种可供选择的有效工具和手段. 相似文献
8.
本文利用最近开发的孔隙介质震电波场算法开展了数值模拟研究,重点探讨不同结构模型中双力偶点源以及有限断层面源作用下产生的地震波及其伴随的电磁信号的时域和频域特征,考察了不同震源时间函数类型对震电波场数值模拟结果的影响,并以Ricker子波为例讨论了震电波场的频率响应特征.结果表明:地震电磁信号的波形、振幅和持续时间等都不同程度受到模型结构的影响,结构越复杂,波形也越复杂;不同震源(双力偶点源与有限断层面源、震源时间函数类型、Ricker子波峰值频率)产生的位移场和电磁场在波形和频谱响应方面整体上呈现相似的特征,但在波形和频谱的幅度和复杂程度方面存在一定差异. 相似文献
9.
地震发生时产生的同震位移、应变以及应力变化,特别是同震应力的变化,在地震触发等问题的研究中有着重要的意义.本文进一步发展了基于均匀弹性水平层状介质,利用广义反射透射系数矩阵和离散波数计算同震位移的方法,使之可以计算相应的应变、应力以及同震库仑应力变化.可适用于多种情况,接收点可以位于地表以及地表以下,震源类型可以是剪切位错源以及拉张位错源.通过与半无限介质的解析解相比较,结果一致,验证了方法的可靠性以及计算精度,可以用于计算地震之后库仑应力变化,为判断余震分布提供参考.在计算同震位移时,使用了梯形积分与Filon积分相结合的积分方式,即提高了同震位移计算的速度,又保证了计算精度,有利于反演问题研究. 相似文献
10.
用Biot介质参数说明了波达波夫震电波方程组中弹性动力学 参数的含义，解释了第一类和第二类震电效应的意义，在忽略第一类震电效应条件下将该方 程组与Pride方程组进行比较，说明了二者在描述第二类震电效应方面的异同点. 同时指出 ：波达波夫方程组忽略了流体与固体的耦合质量；方程中的黏性耗散项丢掉了一个孔隙度因 子，依据该方程组计算出的弹性波和转换电场的幅度将偏大；边界条件之一存在错误，会影 响对波在界面上的反射透射规律的描述. 相似文献
11.
Using numerical modelling, we investigate the evolution of seismoelectric effects induced by seismic excitation in spatially confined lithological units. Typical geometries represent clay lenses embedded in an aquifer or petroleum deposits in a host rock. In fluid‐saturated rocks, seismic waves can generate electromagnetic fields due to electrokinetic coupling mechanisms associated with such processes in the vicinity of the fluid‐mineral interface. Two seismoelectric phenomena are investigated: (1) the co‐seismic field associated with the seismic displacement at each point in a subsurface and (2) the interface response generated at layer boundaries. Our modelling uses a simplified time‐domain formulation of the coupled problem and an efficient 2D finite‐element implementation. To gain insight into the morphogenetic field behaviour of the seismoelectric effects, several numerical simulations for various target geometries were treated. Accordingly, we varied both the thickness of the confined units and the value of the electrical bulk conductivity in porous media. Analysis of these effects shows differences between interface responses for electrically conductive versus resistive units. So the pertinent contrast in electrical bulk conductivity controls the shape and structure of these seismoelectric conversion patterns. Moreover, the seismoelectric interface response captures both the petrophysical and geometrical characteristics of the geological unit. These models demonstrate the value of using seismoelectric interface response for reservoir characterization in either hydrogeological or hydrocarbon exploration studies. 相似文献
12.
地震波在地下含流体孔隙介质中传播时,会引起中观尺度的"局域流",进而产生地震波震电效应.基于Biot(1941)固结理论的准静态方程,在频率域中采用空间有限差分方法,正演模拟虚岩石物理岩样的地震波衰减和震电效应.与时间域虚岩石物理方法相比,该方法既可以直接求取任一频率下的地震波衰减和电势,便于应用于实际岩样的预测分析,也避免了讨论岩样外表面施加的力源函数表达式及时间剖分稳定性条件等问题.首先利用周期性层状介质模型验证了本文所描述方法的有效性,并进一步求取分析了周期性层状介质两种不同特征单元的渗流电流密度及电势,数值模拟结果表明由中观尺度"局域流"引起的震电效应电势振幅数量级在实验室测量范围之内,随后,分析研究了四种不同高渗介质占比值的地震衰减及震电效应特征.最后,将本文提出的震电效应数值计算方法推广至二维,并求取了二维斑块饱和模型的地震波衰减、速度频散、电势的振幅和相位角数值结果. 相似文献
13.
14.
基于孔隙介质动电效应的测量方法在油气勘探领域有着巨大的潜在应用价值.本文针对理论预测的随钻动电测井可削弱钻铤波对地层声波干扰的问题,设计了室内随钻动电测井实验装置和测量系统,并在小尺寸砂岩模型井中开展了随钻动电测井实验研究.为了突出随钻动电测井技术的特点,本文先在砂岩模型井中进行了电缆声波测井和电缆动电测井实验测量,获得了地层声波的传播速度.然后又在同一模型井中进行了随钻声波测井和随钻动电测井实验,记录了这两种情况下的测井全波波形,比较了电缆声波/动电测井和随钻声波/动电测井的差异性,并进一步分析了钻铤波和井中动电转换信号的传播特性.本文还通过随钻动电测井的时域波形提取到被测地层的纵横波速度,从实验角度验证了随钻动电测井技术测量地层波速的可行性.此外,本文实验结果对随钻动电测井的仪器设计及现场测量也有参考意义. 相似文献
15.
A new numerical technique for simulating the coupled seismic and electromagnetic waves in layered porous media 
下载免费PDF全文
下载免费PDF全文 Chen's technique of computing synthetic seismograms,which decomposes every vector with a set of basis of orthogonality and completeness before applying the Luco-Apsel-Chen(LAC)generalized reflection and transmission coefficients method,is confirmed to be efficient in dealing with elastic waves in multi-layered media and accurate in any frequency range.In this article,we extend Chen's technique to the computation of coupled seismic and electromagnetic(EM)waves in layered porous media.Expanding the involved mechanical and electromagnetic fields by a set of scalar and vector wave-function basis,we obtain the fundamental equations which are subsequently solved by using a recently developed version of the LAC generalized reflection and transmission coefficients method.Our approach and corresponding program is validated by reciprocity tests.We also show a numerical example of a two-layer model with an explosion source.The P-to-EM conversion waves radiated from the interface may have potential application. 相似文献
16.
Munirdin Tohti Yibo Wang Evert Slob Yikang Zheng Xu Chang Yi Yao 《Geophysical Prospecting》2020,68(6):1927-1943
Seismoelectric coupling in an electric isotropic and elastic anisotropic medium is developed using a primary–secondary formulation. The anisotropy is of vertical transverse isotropic type and concerns only the poroelastic parameters. Based on our finite difference time domain algorithm, we solve the seismoelectric response to an explosive source. The seismic wavefields are computed as the primary field. The electric field is then obtained as a secondary field by solving the Poisson equation for the electric potential. To test our numerical algorithm, we compared our seismoelectric numerical results with analytical results obtained from Pride's equation. The comparison shows that the numerical solution gives a good approximation to the analytical solution. We then simulate the seismoelectric wavefields in different models. Simulated results show that four types of seismic waves are generated in anisotropic poroelastic medium. These are the fast and slow longitudinal waves and two separable transverse waves. All of these seismic waves generate coseismic electric fields in a homogenous anisotropic poroelastic medium. The tortuosity has an effect on the propagation of the slow longitudinal wave. The snapshot of the slow longitudinal wave has an oval shape when the tortuosity is anisotropic, whereas it has a circular shape when the tortuosity is isotropic. In terms of the Thomsen parameters, the radiation anisotropy of the fast longitudinal wave is more sensitive to the value of ε, while the radiation anisotropy of the transverse wave is more sensitive to the value of δ. 相似文献
17.
对于多孔介质中发生物理化学反应的溶质运移现象,可以用带有非线性源（汇）项作用的一维对流弥散-反应扩散方程来描述,但方程中反映溶质吸附/解吸附能力的源(汇)项系数往往是未知的. 本文讨论了基于出流端浓度观测数据的源项系数反演问题. 根据Tikhonov正则化和矩阵的奇异值分解系统,建立了一种离散的迭代正则化算法,给出了算法实现步骤.数值模拟结果表明该算法不仅精度高,而且对于数据的随机扰动具有稳定性.最后应用建立的算法反演计算了一个具体的土柱试验源项系数,数值结果也表明文中所构造算法的有效性. 相似文献
18.
In real strata anisotropy and viscosity extensively exists. They degraded waveforms in amplitude, resulting in which reducing of image resolution. To obtain high-precision imaging of deep reservoirs, we extended the separated viscous and anisotropic reverse time migration (RTM) to a stable viscoacoustic anisotropic RTM for vertical transverse isotropic (VTI) media, based on single generalized standard and linear solid (GSLS) media theory.. We used a pseudo-spectral method to develop the numerical simulation. By introducing a regularization operator to eliminate the high-frequency instability problem, we built a stable inverse propagator and achieved viscoacoustic VTI media RTM. High-resolution imaging results were obtained after correcting for the effects of anisotropy and viscosity. Synthetic tests verify the validity and accuracy of algorithm. 相似文献
19.
《Advances in water resources》1998,21(6):477-485
A numerical approach for approximating statistical moments of hydraulic heads of variably saturated flows in multi-dimensional porous media is developed. The approximation relies on a first-order Taylor series expansion of a finite element flow model and an adjoint state numerical method for variably saturated flows to evaluate sensitivities. This approach can be employed to analyze uncertainties associated with predictions of head of steady-state or transient flows in variably saturated porous media, with any type of boundary and initial conditions. Limitations of stochastic analytical methods such as spectral/perturbation approaches and the time-consuming Monte Carlo simulation technique are thus alleviated. An example is given to demonstrate the utility of the approach and to investigate the temporal evolution of head variances in a variably saturated flow regime. Results show that the fluctuation of the water table can have significant impacts on the propagation of the head variance. 相似文献