简介:本文提出了一种矿山微地震震源函数反演的新方法。以球坐标系波动方程为基础,对时空域中的微地震观测信号进行层析成像投影到慢度时间域,我们可以获得微地震波成像能量最大值及其对应的震源位置、发震时间和传播速度以及慢度时间信号。通过分析能量最大值处的慢度时间信号与震源函数之间的关系,推导出了计算矿山微地震震源函数的层析成像公式和利用震源函数计算微地震有效辐射能量的公式。进一步利用最小二乘法,将震源函数的振幅谱拟合成模型震源谱,确定了微地震波的零频极限值和拐角频率,最终计算出该微震事件的震源破裂半径。利用这一方法对理论模型数据和实际资料处理,结果表明:本方法具有既能求出任意一个地震事件的震源位置、发震时间和传播速度,又能同时求出震源函数及其频率特征值和震源破裂半径等震源参数的优点,因此,本方法是一种实时快速计算方法。
简介:本文提出利用在不同深度处倾斜放置气枪子阵,对不同深度处单枪或相干枪延时激发的方法压制震源端鬼波,从而改善震源子波分辨率。采用延时激发的方法可以构造同相叠加的首波以及被相对压制掉的鬼波。为了合理评价倾斜震源对于鬼波压制的效果,提出了以实际子波和不受鬼波影响的期望子波频谱之间的归一化平方误差为标准,以深度间隔、气枪组合次序、子阵数目等为变量的最优化问题。针对680cu.in经典枪阵研究表明,在综合考虑去鬼波效果和野外施工可操作性,相邻单枪之间的深度间隔为1m-1.5为较为合理的选择。多个子阵在进行组合时,需要对组合效果进行前期模拟确定最佳深度组合,过多的子阵数目由于会引人大量陷波频率,反而会影响鬼波压制效果,在本次实例研究中,2-3个子阵是较为合理的选择。气枪组合次序会一定程度上影响鬼波压制效果,但不同次序结果差异可以通过简单匹配滤波消除。方向性对比表明多深度倾斜气枪震源能够有效地补偿消除陷波效应的影响,同时改善能量传播的稳定性。
简介:全波形反演是一种高精度的地震成像方法,可以对地下介质物性参数模型进行准确的重构。然而在实际应用中,尤其是在三维复杂介质反演中,计算成本太大是该方法的一个重要缺陷。将混叠震源技术引入到频率域全波形反演中可以大幅度地降低计算成本,提高反演效率。但是使用震源编码技术也带来了两个问题:一方面,参与编码的各个震源之间会产生“串扰噪声”,导致反演结果中出现假象;另一方面,基于震源编码的频率域全波形反演方法周围噪声较为敏感,使该方法对含噪数据反演质量较差。本文引入一种频率组编码方法来压制“串扰噪声”,并基于震源编码技术提出一种频率域自适应全波形反演方法,通过一个与频率相关的自适应选择机制,将常规频率域全波形反演方法和基于震源编码的全波形反演方法联合起来,在保证反演质量的同时也最大程度地提高了反演效率。
简介:描述和估计介质含水量、介电常数等属性参数分布是探地雷达探测技术的重要研究内容。雷达波的旅行时间和反射振幅系数与介质含水量、孔隙度与介电常数密切相关。常规通过旅行时间计算波速以估计介质参数的方法,例如透射波法,共中心点速度分析等,对于复杂介质分辨率有限。基于反射振幅的阻抗反演方法可以直接根据反射系数计算雷达波阻抗以估计介质属性参数,从而有效地避开常规方法在计算波速时精度低的问题。本文首先建立了基于高斯型和指数型混合自相关函数的三维多尺度等效随机介质模型刻画地下随机介质参数分布,并在局部加入高斯椭圆方程描述局部随机异常目标体。其次,通过引入锥形函数以降低随机介质模型在离散网格数值计算方法误差。在此基础上,推导了探地雷达递推阻抗反演的基本流程并结合随机介质模型测试了该方法在复杂介质参数估计中的计算精度。最后,对内蒙地区的实测探地雷达数据利用递推阻抗反演方法来估计地下污染物参数,估计介电常数、含水量结果与钻孔实测数据和同期开展的电阻率成像结果有很好的吻合。说明基于递推阻抗反演方法在探地雷达复杂介质属性参数估计中具有很好的应用前景。
简介:本文以中观孔隙结构的White模型为基础,构建了部分饱和孔隙介质模型,利用Biot方程的建立思路和Johnson推导的体变模量,推导了部分饱和孔隙介质中的纵波方程,并以平面波为例,求取了方程的衰减系数,分析了地震勘探频带范围内地震波的衰减特性。结果表明:在部分饱和孔隙介质中,地震波在低频段也会发生明显的衰减和频散现象,频率越大,衰减越大;且第二纵波的衰减比第一纵波更为明显;这一结论弥补了Biot理论在描述地震勘探频带范围内波的衰减现象的不足。文中还研究了孔隙度、饱和度和模型内径尺寸对第纵波衰减特性的影响机理,主要表现在在地震勘探频带范围内,波的衰减随孔隙度的增大而增大,随含油气饱和度的增大而减小,当孔隙内径尺寸小于二分之一外径尺寸时,波的衰减随内径尺寸的增大而增大,当内径尺寸大于二分之一外径尺寸时,波的衰减随内径尺寸增大而减小。
简介:针对地面场源在地层介质中所产生瞬变电磁场的数值计算方法和响应分析问题,本文展开研究.解法方面,传统离散镜像法多采用复数运算、数字滤波等方法计算量大,针对这些问题,提出一种改进的离散镜像方法:基于Gaver-Stehfest概率变换算法将电磁场解式实数化,选用Prony方法对目标核函数进行指数级数逼近,根据离散镜像原理和近似系数闭合式求解瞬变电磁场.通过试算均质模型中瞬变电磁场并对比汉克尔变换的数字滤波法所得结果,证明该方法有效且具有较好的精度和适用性.继而基于该方法计算地表磁偶源在典型地电模型中产生的瞬变电磁场,对“地面激发一地层中测量”方式所得感应磁场水平分量响应进行分析并得出结论瞬变场水平分量响应与地电结构、观测时间、空间位置等因素有关,感应磁场水平分量响应反映出涡流场分布及其垂向梯度变化情况,在探测异常体的工作中应尽量选取零偏移距、较大偏移距位置钻孔或用较大观测延时以减小背景场对勘测结果的影响。文中所用离散镜像方法与正演计算结论可为相关研究工作提供参考依据。
简介:针对有限差分数值模拟的频散问题,本文将交错网格技术和紧致差分格式相结合,推导了横向各向同性介质一阶速度一应力波动方程的紧致交错网格差分格式;对比分析了紧致交错网格差分格式、交错网格差分格式以及紧致差分格式的截断误差主项,并利用Fourier误差分析方法分析了上述三种差分格式的近似精度;在此基础上,分别采用上述三种差分格式进行了波场数值模拟。结果表明,当差分方程阶数相同时,紧致交错网格差分格式截断误差最小,数值频散最弱,差分精度最高,证实了该方法的有效性。
简介:常规长排列非双曲动校正公式是在VTI介质中得到的,它不能满足任意空间取向TI(ATI)条件下的扩展.本文以VTI介质中非双曲动校正公式为基础,基于我们推导得出的ATI介质中精确四次时差系数解析解和NMO速度解析解,给出ATI介质中长排列优化的非双曲动校正公式.通过与各向异性射线追踪方法计算所得出的"精确走时"结果对比,研究表明优化后的非双曲动校正公式能精确地描述任意强弱、ATI介质中随测线方位变化的走时曲线,可以用来替代耗时、多偏移距、多方位的射线追踪方法正演拟合ATI介质中长偏移距反射走时,为利用非双曲时距的各向异性参数反演提供理论基础性认识。
简介:基于曾新吾和MacBeth提出的时域内双源累积旋转方法,建立了频率相关介质中分析多分量VSP数据的横波分裂参数提取算法(DCTF)。该算法可以在频域中针对单个频率提取横波各向异性参数(快横波的极化方向以及快、慢横波间的时间差),从而避免了目前常用方法中使用窄带通滤波可能带来的误差。通过对地震合成记录的数值分析,确定了该算法的可行性和正确性,并与目前常用方法的应用结果进行了比较。结果显示,频率相关横波分裂参数可以利用DCTF从地震四分量数据中直接提取。在地震频率范围内,含较大尺度裂缝时各向异性参数将表现出频率相关性,这意味着在地震频率范围会出现频散。随着频率的增加,各向异性有降低的趋势。
简介:基于双相介质理论的AVO正演技术是储层性质描述和流体预测的有效技术手段之一,但是输入参数中基质矿物模量的准确性和双相介质模型的的合理性极大地影响双相介质AVO正演效果。因此,本文采用基于流体因子的基质矿物模量反演方法,自适应反演基质矿物体积模量。引入具有岩石物理意义的多约束条件,使得流体替换技术制作的双相介质模型具有岩石物理意义。保证获得的双相介质AVO特征反映实际地层响应,真实可靠。通过不同岩性岩样的对比分析,说明反演方法的优越性和准确性。同时LH地区实际资料应用,获得孔隙度和流体饱和度等重要岩性参数变化时双相介质AVO特征,特别是不同储层孔隙度在同一入射角对应快纵波和横波反射系数幅值的大小差异和突变角差异是分辨储层孔隙度大小的依据。
简介:Inmarineseismicexploration,ocean-bottomcabletechniquesaccuratelyrecordthemulticomponentseismicwavefield;however,theseismicwavepropagationinfluid–solidmediacannotbesimulatedbyasinglewaveequation.Inaddition,whentheseabedinterfaceisirregular,traditionalfinite-differenceschemescannotsimulatetheseismicwavepropagationacrosstheirregularseabedinterface.Therefore,anacoustic–elasticforwardmodelingandvector-basedP-andS-waveseparationmethodisproposed.Inthismethod,wedividethefluid–solidelasticmediawithirregularinterfaceintoorthogonalgridsandmaptheirregularinterfaceintheCartesiancoordinatessystemintoahorizontalinterfaceinthecurvilinearcoordinatessystemofthecomputationaldomainusingcoordinatestransformation.Theacousticandelasticwaveequationsinthecurvilinearcoordinatessystemareappliedtothefluidandsolidmedium,respectively.Attheirregularinterface,thetwoequationsarecombinedintoanacoustic–elasticequationinthecurvilinearcoordinatessystem.Wenextintroduceafullstaggered-gridschemetoimprovethestabilityofthenumericalsimulation.Thus,separateP-andS-waveequationsinthecurvilinearcoordinatessystemarederivedtorealizetheP-andS-waveseparationmethod.
简介:在具有垂直对称轴横向各向同性介质中,利用四种参数来确定中间至远偏移距转换波(C-波)动校正.它们是C-波叠加速度VC2,垂直速度比和有效速度比γ0和γeff,以及各向异性参数χeff.我们将这四种参数作为C波叠加速度模型.C-波速度分析的目的就是确定这种叠加速度模型.C-波叠加速度模型VC2,γ0,γeff,和χeff可以由P-波和C-波反射动校正资料获得.然而错误的传播是C-波反射动校正反演中的严重问题.当前短排列叠加速度由于是从双曲线动校正推算而得,因而其精度不足以为各向异性参数提供有意义的反演值.中间偏移非双曲线动校正不再被人们所勿略,而是可以用一个背景γ加以量化.非双曲线分析通过中间偏移距的γ校正量可以产生VC2,若数据不含燥音,其误差小于1%.方法稳健,允许γ启始假定值的误差达20%.该方法也适用垂直非均匀各向异性介质.精度的提高使能够用4分量地震资料计算各向异性参数.为此提出了两种工作流程:双扫描和单扫描流程.理论数据和实际数据的应用表明这两种流程得出的结果其精度相似,但是单扫描流程比双扫描更有效.