基于FFT-MA的随机介质建模方法          下载免费PDF全文

杨修伟  朱培民  毛宁波  徐志伟  肖丹 《地球物理学报》2018,61(12):5007-5018

随机介质是描述地球介质小尺度非均匀性的有效模型，在地震散射波场分析、储层描述等领域具有广泛应用，快速准确的随机介质建模方法是开展相关研究的基础与前提.本文首次将FFT-MA（Fast Fourier Transform Moving Average）算法引入到随机介质建模研究中，分析了该方法与传统随机介质建模方法相比具有的优势，并提出了基于FFT-MA的非平稳随机介质建模方法.建模实验表明，与传统的基于谱分解定理的随机介质建模方法相比，基于FFT-MA的方法在空间域产生随机数序列，使随机数序列与结构参数分离，因此，随机数序列与所建模型存在空间上的对应关系，可以分区域建模和局部修改模型.在非平稳随机介质模型建模时，滑动窗口的大小可以根据自相关长度变化而变化，避免了每个采样点都建立一次完整大小的模型，提高了建模效率.因此，FFT-MA随机介质建模方法能准确构建满足自相关函数要求的平稳及非平稳随机介质模型，具有建模效率高、灵活、实用的优点.  相似文献   

基于Remez迭代算法求解差分系数的双相介质自适应有限差分正演模拟          下载免费PDF全文

王文化  文晓涛  张波  张懿疆  王威 《地球物理学报》2020,63(6):2400-2414

利用传统有限差分方法对基于Biot理论的双相介质波动方程进行数值求解时,由于慢纵波的存在,数值频散效应较为明显,影响模拟精度.相对于声学近似方程及普通弹性波方程,Biot双相介质波动方程在同等数值求解算法和精度要求条件下,其地震波场正演模拟需要更多的计算时间.本文针对Biot一阶速度-应力方程组发展了一种变阶数优化有限差分数值模拟方法,旨在同时提高其正演模拟的精度和效率.首先结合交错网格差分格式推导Biot方程的数值频散关系式.然后基于Remez迭代算法求取一阶空间偏导数的优化差分系数,并用于Biot方程的交错网格有限差分数值模拟.在此基础上把三类波的平均频散误差参数限制在给定的频散误差阈值和频率范围内,此时优化有限差分算子的长度就能自适应非均匀双相介质模型中的不同速度区间.数值频散曲线分析表明:基于Remez迭代算法的优化有限差分方法相较传统泰勒级数展开方法在大波数范围对频散误差的压制效果更明显;可变阶数的优化有限差分方法能取得与固定阶数优化有限差分方法相近的模拟精度.在均匀介质和河道模型的数值模拟实验中将本文变阶数优化有限差分算法与传统泰勒展开算法、最小二乘优化算法进行比较,进一步证明其在复杂地下介质中的有效性和适用性.  相似文献   

Impact of fractional probability distributions on statistics of hydraulic conductivity,dynamics of groundwater flow and solute transport at a low-permeability site     

Chengpeng Lu  Wei Qin  Ching-Sheng Huang  Yong Zhang  Tema Koketso Ealotswe 《水文研究》2020,34(21):4112-4127

The fractional Brownian motion (fBm) and fractional Lévy motion (fLm) can easily describe the geometry and the statistical structure of hydraulic conductivity (K) for real-world. However, the fBm and fLm models have not been systematically evaluated when building the K field for a low-permeability site. In this study, both the fBm and fLm are used to simulate the low-K field at NingCheGu (NCG), Tianjin, China. Groundwater flow and solute transport are then computed using MODFLOW and MT3DMS, respectively, and the influence of the fBm/fLm models for K on groundwater flow and solute transport is discussed. Results show that the fLm fits better the statistics of the low-K medium than fBm, and the random logarithmic K (LnK) field generated by fLm is more stable because the resultant LnK field captures more of the measured properties at the field site than that generated by fBm. In contrast, the LnK generated by fBm is more likely to form both high-K channels and low-K barriers. The fBm therefore predicts more extreme behaviours in flow and transport, including the preferential flow, low-concentration blocks and solute retention. The overall groundwater renewal period and solute travel time for the fLm simulation are slightly shorter than those for fBm. The impacts of the fLm and fBm models on the statistics of the resultant LnK fields and the dynamics of groundwater flow and solute transport revealed by this study shed light on the selection and evaluation of the fractional probability distribution models in capturing the K fields for low-K media.  相似文献   

VARIATION OF THE TRAVEL TIME DIFFERENCE BETWEEN TELESEISMIC PS CONVERTED WAVE AND PP TRANSMITTED WAVE IN FOCAL REGION BEFORE AND AFTER WENCHUAN M_S8.0 EARTHQUAKE          下载免费PDF全文

SONG Chun-yan  SHAO Xue-zhong  MA Jin 《地震地质》2019,41(3):726-742

When P waves from distant earthquakes meet a velocity discontinuity in the earth's crust and upper mantle, they give rise to a series of converted PS waves besides PP refracted waves. It is possible to monitor the variation of the physical properties of the medium in the limited formation space above the transition zone of the seismogenic zone by measuring the time difference between the teleseismic PS converted wave and the first arrival P wave, that is, time-variation Δt_PS=t_PS-t_P. The advantage of this method is that the transition point of the teleseismic source with similar source is relatively stable at the transition interface, and the accuracy of the measured relative time is high, and the change of the medium in a small range of the seismogenic zone above the conversion interface can be monitored. This paper studies the variation of the travel time difference Δt_PS in focal region before and after Wenchuan M_S8.0 earthquake. We select 2001 to 2012 as the research period, use teleseismic waveforms which occurred in the southern region of Sumatra and Hindu Kush area recorded by Sichuan station YZP and JJS. These teleseisms satisfy 5.0 ≤ M ≤ 6.5, and their waveform signal-to-noise ratio is high with clear initial P-wave motion. The epicentral distance of teleseisms is less than 3 degrees. Then we obtain the variation of the travel time difference Δt_PS between teleseismic PS converted wave and PP transmitted wave recorded during the study period in the two stations. The results show that there is a slow increasing trend of Δt_PS before 2006, and an obvious low value process of Δt_PS appeared in the period about 2 years before the Wenchuan earthquake. The maximum decline was about 0.2~0.3s, more than 4~5 times the measurement error. The low value has a certain degree of return about 2~3 months before the earthquake. The change of arrival time difference indicates that the medium is in different states in different periods of seismogenic process. The sharp decrease in Δt_PS from 2006 to January 2008 may be due to the strong disturbance caused by the stress accumulation of the medium. At this stage, the velocity of P wave and S wave increases with the increase of stress, and the increase of S wave velocity will result in the decrease of Δt_PS. The change of Δt_PS is greatly affected by S wave velocity, so Δt_PS appears to decrease rapidly. Regarding the low value that has a certain degree of return about 2~3 months before the Wenchuan earthquake, the possible reason is that the release of stress is much higher than the accumulation of stress in meta-instability stage. At this stage, the velocity of S wave decreases and the decrease of S wave causes Δt_PS to increase. Then, the Wenchuan earthquake of magnitude 8.0 occurred. It is shown that the teleseismic converted wave method in this paper can monitor the variation of medium's wave velocity before large earthquakes, and it has a good prospect in seismic monitoring and worth further experimental study.  相似文献   

冀鲁豫交界区岩石圈磁场异常分布与磁化率结构分析     

张海洋  李博  苏树朋  王朝景  刘德强  云萌  赵慧琴  王慧 《地震地磁观测与研究》2019,40(6):47-52

利用重磁反演软件,对2017年8月冀鲁豫交界区流动地磁总强度数据进行处理并解析延拓,研究该区岩石圈磁场及其深部、浅表磁场异常分布特征与地下介质磁化率结构的性质。结果表明:冀鲁豫交界区岩石圈磁场异常分布是岩石圈浅表与深部磁异常综合叠加的结果,岩石圈地下介质的浅表与深部磁化率结构具有一定差异性,且受断裂控制呈分段性特征,可能与该区域复杂的岩石圈磁化率结构及地质构造有关。  相似文献   

含裂隙介质中的视电阻率各向异性变化   总被引：1，自引：0，他引：1       下载免费PDF全文

解滔  卢军 《地球物理学报》2020,63(4):1675-1694

我国50多年的视电阻率连续观测结果表明,大地震前近震中区域的视电阻率呈现出与主压应力方位有关的各向异性变化,即:垂直于主压应力方向观测的变化幅度最大,平行方向最小或不明显,斜交方向介于二者之间.目前我国定点台站视电阻率观测的探测范围主要在浅层沉积层以内,通常含有较多的含水裂隙.本文将地下岩土介质简化为由固体基质和含流体/气体裂隙组成的固液气三相介质,且基质、流体和气体具有标量形式的电阻率,推导出了包含基质和流体电阻率、裂隙率、饱和度和裂隙面积率因子的电阻率张量表达式.以裂隙的扩展/闭合表示应力作用下裂隙的变化,得到了电阻率随裂隙变化的微分形式,电阻率变化对裂隙体积变化放大系数的表达式和裂隙横向变化对纵向电阻率影响的横向权系数的表达式.在此基础上得到了介质电阻率和视电阻率的各向异性变化特征:对于含水裂隙介质,无论裂隙如何变化,均是最小主轴方向电阻率的变化幅度大于其他方向;对于含水孔隙介质,沿孔隙主要变化方向的主轴电阻率变化幅度大于其他方向.对于各向异性变化,视电阻率和介质电阻率存在π/2的方向差异.相较于含水岩石,无水岩石介质电阻率的各向异性变化不显著.本文提出的电阻率表达式可以对实验室和野外实际观测的许多结果做出合理的解释.  相似文献   

基于Ruger近似式预测裂缝型储层的可行性研究     

吴姗姗  高刚  桂志先  王鹏  李成毅 《煤田地质与勘探》2020,48(1):189-196

叠前纵波方位各向异性检测裂缝是目前应用最为广泛的方法，该方法主要是利用AVAZ地震资料椭圆拟合的长轴与短轴信息来评价裂缝，该方法在具体应用时，认为裂缝型储层AVAZ响应特征的变化仅由各向异性参数决定，而没有考虑非各向异性参数的影响，从而引起裂缝评价的多解性。为提高该技术的裂缝评价精度，从模型分析入手，将模型中各向异性参数及岩性参数设置为概率密度分布函数，然后采用Monte Carlo随机方法进行叠前AVAZ正演模拟。首先，对各向异性参数对椭圆拟合的影响进行简单分析，得出各向异性参数γ对椭圆扁率B/A与各向异性因子B影响最大、δ次之、ε最小。然后，进一步重点分析不同标准差的速度与密度模型的B/A与B响应特征，得出地层的速度与密度的变化在较小范围内时，该技术评价裂缝才具有一定的可行性，其中，地层纵波和横波速度的变化（尤其是纵波速度）对叠前AVAZ响应影响较大，地层密度几乎没有影响，故当横向岩性变化较大时，该方法在评价裂缝型储层时的结果具有一定的不确定性。   相似文献   

一种变形监测中中长基线的单历元模糊度解算算法     

温亚鑫  戴吾蛟 《大地测量与地球动力学》2019,39(9):942-946

针对已有单历元模糊度解算算法只适用于短基线的情形，根据误差传播率确定模糊度精度，依据误差确定模糊度搜索空间，提出一种利用GPS双频观测信息实现中长基线下单历元模糊度固定的算法。首先，利用电离层约束法得到双差模糊度之间的线性约束条件；其次，利用变形监测中监测站点坐标近似已知的条件进行约束，反算获得近似模糊度；最后，根据监测站点已有的观测数据精度给出搜索空间，搜索固定模糊度。结果表明，该方法可以实现中长基线下单历元的模糊度固定，对实时变形监测有一定的参考价值。  相似文献   

A stress and displacement discontinuity element method for elastic transversely anisotropic rock     

A. Bobet  V. García Marín 《国际地质力学数值与分析法杂志》2014,38(18):1898-1922

The paper presents closed‐form solutions for stress and displacement influence functions for stress discontinuity (SD) and displacement discontinuity (DD) elements, for a two‐dimensional plane‐strain elastic, transversely anisotropic medium. The solutions for SD elements are based on Kelvin's problem and for DD elements on the concept of dipoles. Stress and displacement influence functions are derived for the following elements: constant SD, linear SD, constant DD, linear DD, square root DD, parabolic DD, constant DD surface, and linear DD surface elements. The formulations are incorporated into FROCK, a hybridized boundary element method code, and are validated by providing comparisons between the results from FROCK and the finite element code ABAQUS. A limited parametric analysis shows the effects of slight anisotropy on the stress field around the tip of a crack and of the orientation of the crack with respect to the axes of elastic symmetry. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

黏声介质平面波有限差分叠前深度偏移及在莺歌海盆地的应用     

李添才  谢玉洪  李列  王华忠  王新领  但志伟  胡江涛 《地球物理学报》2014,57(5):1612-1622

莺琼盆地诸探区中存在底劈现象，深层气源产生的气体沿底劈产生的裂隙通道向上漫溢.漫溢过程中，一是充填在遇到的砂体中形成气藏；二是弥漫在上溢通道中，使得通道中的纵波速度发生变化，进而纵波波阻抗差异变小，反射变弱.另外，通道中气体的存在，会加强地层的吸收衰减，使得地震波振幅变弱、高频成分损失导致同相轴分辨率降低.利用OBC数据进行多波地震勘探和利用黏声介质的叠前深度偏移都是改善模糊区成像质量的重要方法技术.为此，本文提出用黏声介质平面波有限差分法叠前深度偏移成像方法改善气体充填区域的成像质量.黏声介质成像目的是补偿地震波的吸收衰减；平面波偏移成像目的是适应海上单炮数量巨大，提高波动方程叠前深度偏移成像的效率；有限差分法叠前偏移的目的是适应该区浅层气分布局域性极强、Q值的空间变化大的情况.在莺歌海某探区的实际数据上的黏声介质平面波有限差分叠前深度偏移试验证明，本方法是改善模糊区成像质量的较为有效的途径.  相似文献