共查询到20条相似文献，搜索用时 125 毫秒
1.
本文介绍一种电磁波衰减层析成像的方法。这种成象方法利用电磁波振幅谱重心频率的下移作为投影函数来重建本征衰减系数（目标函数）的空间分布。作者已经将帝种方法应用于处理地质雷达的预测资料。对于大多数天然介质而言。当电磁波信号的频率不高于过滤频率（介质电导率与电容率之比）时，电磁波的衰减会随着频率的增高而增强。因而，接收到的电磁波信号的振幅谱重心频率将发生向低频方向的飘移。这种下移正比于本征衰减系数沿射线 相似文献
2.
3.
Lg波的QLg值是描述区域地壳结构及介质衰减特性的重要参数之一,QLg层析成像被广泛应用于地壳衰减结构横向不均匀性的研究中.但是,对QLg在垂直方向上的不均匀性的研究较少.当发现一个地区发生QLg值低时,我们希望进一步了解这种介质衰减深度.本文通过数值模拟合成地震图的方法,通过分析不同震源深度和介质模型计算QLg,判断Lg波发生主要衰减的深度.研究结果显示：（1）当震源深度浅时,上地壳介质衰减对于Lg波Q值的影响明显比下地壳的影响更大.而震源深度深时,两者贡献基本相当;（2）如果上地壳介质衰减强,则随着震源深度的增加,η值逐渐减小;如果下地壳介质衰减强,η值先增大后减小;但如果下地壳存在低速层,则η值持续增大;（3）如果上地壳介质衰减强,则随着震源深度的增加,Q0逐渐增大;反之如果下地壳介质衰减强,则Q0逐渐减小.因此,我们可以通过调查某一区域不同深度地震的Lg波衰减规律,来判断这一地区的地壳介质衰减情况. 相似文献
4.
地震波在非弹性介质中的衰减效应常用品质因子Q度量.相对准确的Q模型对提高强衰减介质中地震波成像的质量至关重要.本文提出了黏声介质反射波形反演(QRWI)方法来重建地下宏观Q模型.在缺乏大偏移距和低频地震数据时,该方法以黏声波方程为波场传播引擎,利用反射波核函数对模型中深部的敏感性去提取背景Q值.当速度高、低波数成分均已知时,基于波形拟合的QRWI可以获得较高分辨率的反演结果.由于地下介质速度的高波数扰动很难准确估计,本文通过引入峰值频移目标函数,极大地降低了QRWI对速度高波数成分的依赖.理论合成数据实验结果表明,本文方法反演得到的宏观Q模型可以满足衰减补偿逆时偏移成像的要求. 相似文献
5.
自旋回波作为测量T2时间的重要脉冲序列,由该脉冲系列交变磁场激发而获得的核磁共振信号不受磁场不均匀性的影响,这对于以地磁场作为稳定磁场的SNMR方法来说,可以减小地磁场不均匀对NMR信号的影响,得到真实、可信的T2值,提高信噪比。将自旋回波信号应用到地面核磁共振找水方法中求取T2分布,在国内尚属首例,本文借鉴核磁共振测井中的多指数反演方法,采用奇异值分解算法编程,反演得到T2分布,并与SNMR反演软件的结果T2*对比。结果证明,用SE信号提取T2时间,可以减少地磁场不均匀对NMR信号的影响,更真实地反映T2分布,提高信噪比。 相似文献
6.
传统基于地震背景噪声互相关的层析成像技术,主要利用了互相关函数的相位信息,即利用走时来反演地球介质的速度和各向异性特征.最近,开始有研究人员利用互相关函数的振幅信息提取介质的衰减.主要是将二维弹性情形下,互相关函数正比于第一类零阶贝塞尔函数J_0(k_0r)这一结论,通过引入指数衰减系数,直接推广到衰减介质情形,令衰减介质中背景噪声互相关函数类比于J_0(k_0r)·e~(-α(ω)r),以此来反演介质的衰减系数α.然而,在衰减介质中,互相关受源的方位平均的影响,这种简单的推广,可能无法提取可靠的衰减系数.本文基于平面波的叠加模型,研究不同的互相关定义和坐标选择下,衰减介质中两点间互相关函数的理论表达式.结果表明,在平面波叠加模型下,互相关函数的表达形式随着坐标原点的选择,并因而随着源分布的变化而变化,对不同的归一化因子,表达式也不尽相同.利用J_0(k_0r)·e~(-α(ω)r)的形式拟合背景噪声的观测数据得到的衰减比实际值偏小. 相似文献
7.
地震走时层析成像是反演地层各向异性参数分布的有效方法,但是关于地震各向异性介质走时层析成像的研究并不多,其技术远远没有达到成熟的阶段.在野外数据采集时,地表反射波观测方式相对井间和垂直地震剖面观测方式的成本更低,利用qP反射波走时反演各向异性参数具有更加广泛的实用价值.本文实现的TI介质地震走时层析成像方法结合了TI介质反射波射线追踪算法、走时扰动方程和非线性共轭梯度算法,它可以对任意强度的TI介质模型进行反演,文中尝试利用qP反射波走时重建TI介质模型的参数图像.利用qP反射波对层状介质模型和块状异常体模型进行走时反演,由于qP波相速度对弹性模量参数和Thomsen参数的偏微分不同,所以可以分别反演弹性模量参数和Thomsen参数.数值模拟结果表明:利用qP反射波可以反演出TI介质模型的弹性模量参数与Thomsen参数,不同模型的走时迭代反演达到了较好的收敛效果,与各向同性介质走时反演结果相比较,各向异性介质走时反演结果具有较好的识别能力. 相似文献
8.
在利用地震波数据进行地球物理反演时,地层对地震波的吸收衰减效应会对地层物性参数的准确反演产生较大的影响,因此利用黏弹性声波方程进行反演更符合实际情形.本文在考虑地层衰减效应进行频率空间域正演模拟的基础上,提出基于黏弹性声波方程的频率域逆散射反演算法并对地震波传播速度进行反演重建,在反演过程中分别用地震波传播复速度和实速度来表征是否考虑地层吸收衰减效应.基于反演参数总变差的正则化处理使反演更加稳定,在反演中将低频反演速度模型作为高频反演的背景模型进行逐频反演,由于单频反演过程中背景模型保持不变,故该方法不需要在每次迭代中重新构造正演算子,具有较高的反演效率;此外本文在反演过程中采用了基于MPI的并行计算策略,进一步提高了反演计算的效率.在二维算例中分别对是否考虑地层吸收衰减效应进行了地震波速度反演,反演结果表明考虑衰减效应可以得到与真实模型更加接近的速度分布结果,相反则无法得到正确的地震波速度重建结果.本文算法对复杂地质模型中浅层可以反演得到分辨率较高的速度模型,为其他地震数据处理提供比较准确的速度信息,在地层深部由于地震波能量衰减导致反演分辨率不太理想. 相似文献
9.
以分层半空间内部含有一层孔隙介质为物理模型进行数值计算,研究半空间表面瑞利波的传播和衰减特性.为更加接近实际,结合瑞利波的激发特性,确定了瑞利波的主衰减曲线,并主要以此进行规律分析.针对速度递增和含低速层这两种典型的地质模型,讨论了瑞利波的传播衰减在不同地质模型下的特性,并分析了各自的规律.结果表明,在这两种模型下瑞利波的主衰减曲线都受孔隙介质所处空间位置影响产生比较明显的变化,但衰减系数极大值对应的波长与模型的表层厚度存在较明显的线性对应关系,利用这一关系,可以在实际勘探中快速得到表层介质厚度.另外,通过对比分析还可以看到,瑞利波主衰减曲线随孔隙介质的孔隙度和渗透率的变化都强于主频散曲线的变化,表明衰减曲线对孔隙度和渗透率的变化更加敏感,理论上更加适合进行介质参数反演工作.综合对比结果,我们认为瑞利波主衰减曲线中包含了更丰富的介质参数信息,如果能够有效利用,将可以提高瑞利波勘探的准确性和应用范围. 相似文献
10.
1997年11月8日西藏Mw7.5级玛尼地震是干涉雷达技术应用于地震观测以来的一次重要事件.在第一部分中，我们应用广泛使用的Okada线弹性位错模型，假设断层的各个分段滑动量均匀，反演得到断层各个分段的几何参数和均匀滑动量.本部分的反演进一步去除滑动均匀假设，并利用更能反映断层真实状态的角形元位错模型（线弹性），在第一部分反演得到断层几何的基础上，反演断层面的静态位错分布.反演结果表明，线弹性滑动分布模型能够更好地解释观测数据，进一步提高反演的数据拟合程度.最终得到了断层面上的走滑和倾滑位错分布.首次得到的断层面滑动分布显示断层面滑动在浅部（0～12 km）比较集中，地震破裂长度约170 km，最大左旋走滑位移达4.8 m；反演结果还表明局部段落存在较大倾滑位移，量值达到1.9 m，这在断层模型中是不能忽略的，它可能是断层两侧形变不对称的原因之一；反演得到的标量地震矩为2.18×1020 N·m，相当于矩震级Mw7.5，与Velasco等利用地震波形反演得到的结果一致. 相似文献
11.
In this paper, we develop a methodology for early detection of potential CO2 leakage from geological storage formations using pressure and surface-deformation anomalies. The basic idea is based on the fact that leakage-induced pressure signals travel much faster than the migrating CO2; thus such anomalies may be detected early enough for risk management measures taking effect in avoiding substantial CO2 leaks. The early detection methodology involves automatic inversion of anomalous brine leakage signals with efficient forward pressure and surface-deformation modeling tools to estimate the location and permeability of leaky features in the caprock. We conduct a global sensitivity analysis to better understand under which conditions pressure anomalies can be clearly identified as leakage signals, and evaluate signal detectability for a broad parameter range considering different detection limits and levels of data noise. The inverse methodology is then applied to two synthetic examples of idealized two-aquifer-and-one aquitard storage systems, with an injection well and a leaky well, for different monitoring scenarios. In Example 1, only pressure data at the monitoring and injection wells are used for leakage detection. Our results show that the accuracy of leakage detection greatly depends on the level of pressure data noise. In Example 2, joint inversion of pressure and surface-deformation measurements significantly improves the speed of convergence toward the true solution of the leakage parameters and enables early leakage detection. In both examples, successful detection is achieved when two monitoring wells are appropriately placed within up to 4 km from the leaky well. 相似文献
12.
13.
In this paper, we present a new method for seismic stratigraphic absorption compensation based on the adaptive molecular decomposition.
Using this method, we can remove most of the effects resulting from wavelets truncation and interference which usually exist
in the common time-frequency absorption compensation method. Based on the assumption that the amplitude spectrum of the source
wavelet is smooth, we first construct a set of adaptive Gabor frames based on the time-variant properties of the seismic signal
to transform the signal into the time-frequency domain and then extract the slowly varying component (the wavelet’s time-varying
amplitude spectrum) in each window in the time-frequency domain. Then we invert the absorption compensation filter parameters
with an objective function defined using the correlation coefficients in each window to get the corresponding compensation
filters. Finally, we use these filters to compensate the time-frequency spectrum in each window and then transform the time-frequency
spectrum to the time domain to obtain the absorption-compensated signal. By using adaptive molecular decomposition, this method
can adapt to isolated and overlapped seismic signals from the complex layers in the inhomogeneous viscoelastic medium. The
viability of the method is verified by synthetic and real data sets. 相似文献
14.
PRINCIPAL COMPONENT ANALYSIS AND LOCAL CORRELA-TION TRACKING AS TOOLS FOR REVEALING AND ANALYZING SEISMO-ELECTROMAGNETIC SIGNAL OF EARTHQAUKE 
下载免费PDF全文
下载免费PDF全文 This study provides new seismo-electromagnetic data processing methods to extract the anomalous signals by combining the principal component analysis(PCA)and local correlation tracking(LCT)methods. The PCA method can separate signals of different frequencies by projecting them to different axes according to their energy. So it can solve the problems of identifying the relatively weak signals in strong interference background to a certain extent. The LCT method is more suitable for non-stationary signal processing compared with classical cross-correlation method. This method is based on the good spatial correlation between the magnetic field components of different ELF stations to pick up the correlation coefficient, so as to achieve the purpose of weak anomalies signals identification. As a case study of the M4.6 Jinggu earthquake in Yunnan, China, we investigated the electromagnetic data observed by ELF stations near the epicenter. First, we applied the PCA method to the magnetic-filed data and got the temporal variation of percentage of each principal component. The results indicate that the contribution of the second principal components, which may relate with the earthquake, increased significantly about a week before the earthquake. Then, we applied the LCT method to the magnetic-filed data processing as well, and the results of both north-south and east-west magnetic field components showed that local correlation coefficient saw anomalies about a week before the earthquake, which had a good corresponding relationship with the former results by PCA method. Both means for the same earthquake case got a consistent conclusion, which not only enhanced the reliability of the results, but also confirmed the effectiveness of two methods applied in the earthquake-related electromagnetic anomalies extraction. We also discussed the possible relationship between these anomalies and the earthquake. Although there are no direct evidence and supporting theories in terms of the relationship between abnormal electromagnetic signals and earthquakes at present, the studies in this paper may strengthen the understanding of seismoelectromagnetic phenomena and promote further research. 相似文献
15.
本文试图采用卫星重力资料和一种新的反演方法来研究地幔的横向密度异常分布.先将密度异常△（r,,（?））在一个三维正交函数系下进行展开,其展开系数待定.然后,根据密度异常与重力扰动位之间的关系建立观测方程组,其中未知向量由密度异常展开系数组成,重力扰动观测向量由 GEM10B 重力模型中的位系数计算而得,并通过适当选取重力位系数的阶数,对观测向量进行滤波.最后,就下地幔（670km——CM 界面）作了实际计算.计算中,重力扰动位阶数取为2——11阶,密度异常展开式的截断阶数取为 K=4和 L=6,求解观测方程组时采用阻尼最小二乘法.结果表明:密度扰动值在670km 不连续面及核幔界面处达到极大值,且在环太平洋地区存在一高密度带,太平洋中部对应于一低密度区,这些特征与 Dziewonski 得到的下地幔三维波速异常分布特征相一致.但是,在南极地区、大西洋及印度洋部分地区,所得的密度异常分布与三维波速异常分布呈负相关,文章就其原因作了初步分析. 相似文献
16.
17.
Noise suppression for microseismic data by non‐subsampled shearlet transform based on singular value decomposition 下载免费PDF全文
下载免费PDF全文 The existence of strong random noise in surface microseismic data may decrease the utility of these data. Non‐subsampled shearlet transform can effectively suppress noise by properly setting a threshold to the non‐subsampled shearlet transform coefficients. However, when the signal‐to‐noise ratio of data is low, the coefficients related to the noise are very close to the coefficients associated with signals in the non‐subsampled shearlet transform domain that the coefficients related to the noise will be retained and be treated as signals. Therefore, we need to minimise the overlapping coefficients before thresholding. In this paper, a singular value decomposition algorithm is introduced to the non‐subsampled shearlet transform coefficients, and low‐rank approximation reconstructs each non‐subsampled shearlet transform coefficient matrix in the singular value decomposition domain. The non‐subsampled shearlet transform coefficients of signals have bigger singular values than those of the random noise, which implies that the non‐subsampled shearlet transform coefficients can be well estimated by taking only a few largest singular values. Therefore, those properties of singular value decomposition may significantly help minimise overlapping of noise and signals coefficients in the non‐subsampled shearlet transform domain. Finally, the denoised microseismic data are obtained easily by giving a simple threshold to the reconstructed coefficient matrix. The performance of the proposed method is evaluated on both synthetic and field microseismic data. The experimental results illustrate that the proposed method can eliminate random noise and preserve signals of interest more effectively. 相似文献
18.
本文从麦氏方程出发,运用电流磁场的等效理论阐明了在磁电勘探中,当矿体的二次极化电流具有垂直对称轴分布时,在地表面及以上空间中任一点的磁异常场（H）均等于零。因而,具有任何垂直对称状电化学矿产模式（石油、天然气、地下水及金属矿床等）的磁电勘探法异常场线积分均不存在。此外,由于对各种形状和任何极化倾角的矿体情况,地表面上均无电流密度垂直分量（in）存在,故地表面上异常场线积分均为零（∮L H·dL=0）。因此,磁电勘探法不具备在地面上寻找此种矿产的物理前提。 本文还对倾斜极化球体和柱体的磁电勘探异常进行了计算,给出了主剖面上理论曲线的异常分布规律,指出了磁电异常的数量级。同时,还对井中磁电异常的分布规律给出了计算资料。为了对比,也对相应条件下的自然电场异常做了计算,指出了电场异常与磁场异常之间的不同特点及其对实际找矿的意义。 相似文献
19.
地幔低速异常体的精确探测对于研究热点/地幔柱的成因至关重要.但是,由于观测手段以及波前愈合现象的影响,远震走时层析成像在探测地球深部的低速异常体时具有极大的限制.因此,越来越多的学者利用地震波穿过低速体后产生的绕射波和直达波.干涉造成的波形畸变信息来研究地幔低速体.本文通过三维地震数值模拟来拟合长白山、黄石公园地区宽频地震记录,分析绕射波和直达波相互干扰造成的波形畸变现象.模拟结果显示,绕射波/波形畸变对于探测地幔深处的低速体具有较高分辨率;目前常用的基于互相关测量的有限频层析成像方法对于深部低速体速度异常低估现象明显.如果将有限频层析成像结果在深部的速度异常值放大到原来的2倍,数值模拟可以很好的拟合观测到的地震波形畸变现象.通过对长白山地区、黄石公园地区记录的地震波形畸变现象的拟合,我们进一步确认了长白山、黄石地区下方分别在400~600 km、700~1000 km处存在地幔低速体,并对低速体速度异常值进行了矫正. 相似文献
20.
地震波在各向异性介质中以一个准P波(qP)和两个准S波(qS1和qS2)的形式传播.研究三种波的相速度、群速度以及偏振方向等传播性质能够为各向异性介质中的正反演问题提供有效支撑.具有比横向各向同性(TI)介质更一般对称性的正交各向异性介质通常需要9个独立参数对其进行描述,这使得对传播特征的计算更为复杂.当两个准S波速度相近时具有耦合性,从而令慢度的计算产生奇异性.因此,奇异点(慢度面的鞍点和交叉点)附近的反射与透射(R/T)系数的求解不稳定,会导致波场振幅不准确.本文首次通过结合耦合S波射线理论和基于迭代的各向异性相速度与偏振矢量的高阶近似解,得到了适用于正交各向异性介质以qP波入射所产生的二阶R/T系数的计算方法.与基于一阶近似的结果相比,基于二阶近似的方法提高了qP波R/T系数的精度,能得到一阶耦合近似无法表达的准确的qP-qS转换波的R/T系数解,且方法适用于较强的各向异性介质. 相似文献