简介:本文提出的储层物性参数同步反演是一种高分辨率的非线性反演方法,该方法综合利用岩石物理和地质统计先验信息,在贝叶斯理论框架下,首先通过变差结构分析得到合理的变差函数,进而利用快速傅里叶滑动平均模拟算法(FastFourierTransformMovingAverage,FFT-MA)和逐渐变形算法(GradualDeformationMethod,GDM)得到基于地质统计学的储层物性参数先验信息,然后根据统计岩石物理模型建立弹性参数与储层物性参数之间的关系,构建似然函数,最终利用Metropolis算法实现后验概率密度的抽样,得到物性参数反演结果。并将此方法处理了中国陆上探区的一块实际资料,本方法的反演结果具有较高的分辨率,与测井数据吻合度较高;由于可以直接反演储层物性参数,避免了误差的累积,大大减少了不确定性的传递,且计算效率较高。
简介:常规多波联合反演采用Zoeppritz方程的近似式构建正演方程,反演过程中需要假定背景纵横波速度比为常数,其反演精度不高,稳定性不好。本文提出了一种基于精确Zoeppritz方程的多波联合反演方法,结合贝叶斯方法进行广义线性反演。本方法基于精确Zoeppritz方程构建正演方程,避免了近似式反演在大角度时引起的误差;利用贝叶斯方法引入模型参数的先验分布信息,作为反演的正则化项,降低了反演的不适定性;反演目标函数中引入低频软约束,稳定了反演低频结果,提高了反演的鲁棒性;在求解反演目标函数时,利用快速算法,降低了反演的运算量。经过模型试算,证明了该方法的优越性和抗噪性;并在实际资料的应用中证明了该方法的实用性和有效性。
简介:本文以中观孔隙结构的White模型为基础,构建了部分饱和孔隙介质模型,利用Biot方程的建立思路和Johnson推导的体变模量,推导了部分饱和孔隙介质中的纵波方程,并以平面波为例,求取了方程的衰减系数,分析了地震勘探频带范围内地震波的衰减特性。结果表明:在部分饱和孔隙介质中,地震波在低频段也会发生明显的衰减和频散现象,频率越大,衰减越大;且第二纵波的衰减比第一纵波更为明显;这一结论弥补了Biot理论在描述地震勘探频带范围内波的衰减现象的不足。文中还研究了孔隙度、饱和度和模型内径尺寸对第纵波衰减特性的影响机理,主要表现在在地震勘探频带范围内,波的衰减随孔隙度的增大而增大,随含油气饱和度的增大而减小,当孔隙内径尺寸小于二分之一外径尺寸时,波的衰减随内径尺寸的增大而增大,当内径尺寸大于二分之一外径尺寸时,波的衰减随内径尺寸增大而减小。
简介:用基于张量格林函数的体积分方程法对三维异常体进行瞬变电磁响应的正演模拟,首先在频率域内计算电磁场分量的频率域响应,然后利用快速数字滤波技术将计算结果转换到时间域。设计和计算了水平电偶极子源激发下层状水平地层模型背景下的常见地形如山谷、山峰地形的模型,并考察分别把源和接收器放于这些地形中的瞬变电磁场响应,详细分析了这些地形对长偏移距瞬变电磁测深(LOTEM)的影响。结果表明,山谷和山峰地形对LOTEM的结果均有不同程度的影响。当电偶极子源放在山谷谷底时,地形对观测异常场的畸变非常严重;当接收器放在山谷中时,接收器处地形的影响强烈但该影响在空间和时间上只是局部的。总体来讲,不论山峰地形位于何处,其对LOTEM的影响相对较小。当地形处于发射源与接收器之间时,地形对LOTEM的影响非常小,表明在进行LOTEM勘探时,选择发射源的放置比接收器的位置更加重要,野外勘探是尽量把发射源选择在开阔的平坦位置。
简介:频率-波数域单程波算子能高效地模拟地震波在复杂介质中的传播,但是在描述波的大角度传播和速度横向扰动变化较大介质中传播的问题时仍然存在一定误差。这类误差是由于对单平方根算子使用Taylor展开式的近似程度不足所造成。为了进一步提高泰勒展开式的精确性,本文提出一种利用粒子群智能算法优化级数展开系数的高阶广义屏算子对单平方根算子的展开级数进行优化处理。新的偏移算法能在保持单程波偏移算法高效的前提下进一步提高偏移算子在大角度的成像精度和对强横向速度变化介质的适应性。通过脉冲响应实验,验证了基于粒子群算法优化级数的高阶广义屏算子能够提高常规的高阶广义屏算子的成像精度和成像角度。根据对二维SEG/EAGE盐丘模型的成像处理,基于粒子群算法优化级数的高阶广义屏算子对盐丘下面的断层取得了更高质量的成像,说明粒子群优化级数的高阶广义屏算子比常规的高阶广义屏算子具有更好的横向速度适应性。为了检验本文所提算法对实际资料的处理能力,我们利用常规的偏移处理技术和本文所提算法对一条海上二维数据进行了偏移成像处理,对比分析成像剖面发现本文所提算法描述了更加清晰的层位信息和更高质量的偏移剖面。本文所提算法能有效提高高阶广义屏偏移在广角度成像的能力,具有一定实际应用价值。
简介:提出一种对整幅高光谱图像的稀疏表示结果进行直接显示的方法,图中不仅包含了稀疏表示中保留的光谱信息,还可显示整幅图像的空间信息。稀疏表示后,将字典中的各有效原子根据光谱特性选择颜色标签,之后根据稀疏系数进行混合颜色显示,此时的图像能够同时满足可分性及距离保持特性。针对局部地物时,提出的单像素混合阵列表示法及改进的裂片纹理技术能够直观且完整的显示出每个像元的具体组成情况,还能够根据所生成图像中的信息对原始HSI进行重建,进而提高数据的利用率。该模型不仅能够良好地显示地物的空间特性,同时能够显示稀疏系数的组成,同时单像素混合阵列表示法及裂片纹理技术弥补了混合像素彩色显示中颜色表达混乱的弊端。对真实地物数据进行实验,结果证明该模型产生的彩色图像具有良好的视觉效果及可分性,满足距离保持特性。
简介:针对转换波各向异性叠前时间偏移实践中的四参数速度模型估计中参数间相互影响、难于准确的确定,本文将简化的转换波两参数动校正方程从叠加速度分析拓展到叠前时间偏移速度模型修正,形成了基于两参数动校正方程的新的四参数法转换波叠前时间偏移速度修正新方法和流程。在该方法和流程中,先进行转换波两参数叠加速度分析,获得转换波叠前时间偏移初始的速度和各向异性参数,然后通过分析偏移后的共成像点道集中的剩余时差修正速度和各向异性参数。其中叠前时间偏移速度模型的垂直速度比仍是需要利用纵波和转换波叠加剖面,层位标定后,利用相关的方法获得。有效速度比需结合纵波速度分析,利用Thomsen(1999)公式计算得到初始模型,最终依靠百分比扫描偏移处理获得最终的有效速度比模型。该方法简化了转换波高质量成像的速度估计方法,减小了多个参数估计的不确定性,在实际应用中也取得了较好的成像效果。
简介:为了解决海底起伏变化对地震波场的影响问题,本文提出将(x-z)域中的曲网格映射成(ξ,η)域中的矩形网格,推导出(ξ,η)域中的二维标量声波方程,根据推导出来的波动方程采用逆时有限差分法将海面上采集到的地震波场在(ξ,η)域中向下延拓至海底面,延拓时采用海水的速度,然后采用顺时有限差分法将延拓后的地震波场再反延拓到海面上,延拓时采用海底面以下地层的速度,从而消除了海底起伏带来的负面影响。模型及实际地震资料的计算分析表明该方法不但能够校正由于海底起伏所引起的海底面下地层反射波场的不连续性还能够校正由于海底起伏所引起的地震波的动力学特征的变化。对延拓前后的地震波场进行速度反演,延拓后反演的地层速度比延拓前反演的地层速度的精度提高很多,延拓前后地震波场的叠加剖面对比表明该延拓方法能够明显提高地震波场的成像质量。
简介:波动方程有限差分法是地震数值模拟中的一种重要的方法,对理解和分析地震传播规律、分析地震属性和解释地震资料有着非常重要的意义。但是有限差分法由于其离散化的思想,产生了不稳定性。精细积分法在有限差分法的基础上,在时间域采用解析解的表达形式,在空间域保留任意差分格式,发展成为半解析的数值方法。本文结合并发展了以往学者的成果,推导了任意精细积分法的三维弹性波正演模拟计算公式,并对其稳定性进行了数值分析。在计算实例中,实现了精细积分法二维和三维弹性波模型的地震正演模拟,对计算结果的分析表明,精细积分法反射信号走时准确,稳定性好,弹性波场相较于声波波场,弹性波波场成分更为丰富,包含了更多波型成分(PP-和PS-反射波、透射波和绕射波),这对实际地震资料的解释和储层分析有重要的意义。实践证明,该方法可直接应用到弹性波的地质模型的数值模拟中。
简介:本文提出利用在不同深度处倾斜放置气枪子阵,对不同深度处单枪或相干枪延时激发的方法压制震源端鬼波,从而改善震源子波分辨率。采用延时激发的方法可以构造同相叠加的首波以及被相对压制掉的鬼波。为了合理评价倾斜震源对于鬼波压制的效果,提出了以实际子波和不受鬼波影响的期望子波频谱之间的归一化平方误差为标准,以深度间隔、气枪组合次序、子阵数目等为变量的最优化问题。针对680cu.in经典枪阵研究表明,在综合考虑去鬼波效果和野外施工可操作性,相邻单枪之间的深度间隔为1m-1.5为较为合理的选择。多个子阵在进行组合时,需要对组合效果进行前期模拟确定最佳深度组合,过多的子阵数目由于会引人大量陷波频率,反而会影响鬼波压制效果,在本次实例研究中,2-3个子阵是较为合理的选择。气枪组合次序会一定程度上影响鬼波压制效果,但不同次序结果差异可以通过简单匹配滤波消除。方向性对比表明多深度倾斜气枪震源能够有效地补偿消除陷波效应的影响,同时改善能量传播的稳定性。
简介:随着煤矿的深部开采,多层积水采空区探测成为煤田水文物探工作的新内容。中心回线瞬变电磁(简称TEM)法因场结构与地层的耦合关系,对低阻层探测有利。但低阻层的屏蔽作用不仅使得探测同样的深度需要更长的观测时间,而且还会减弱下伏地层的异常响应。本文通过直接时域数值模拟和水平分层大地的模型正演,估算了探测目的层所需要的时间长度,根据噪声对观测数据造成的影响,给出了多层积水采空区可分辨的标准。山西大同达子沟煤矿水文勘探中获得的TEM实测曲线,表现了对多层积水采空区的探测能力。在实测曲线的定量反演解释中,利用电测井资料作为初始参数进行定量反演,约束了煤系薄层的等值性。所推断的三层和两层积水采空区,为钻孔所证实。研究结果表明,当观测时间有足够的长度、下伏地层的异常显示大于观测误差,中心回线TEM法探测多层采空区积水是可行性的。
简介:瑞利波勘探主要是建立在弹性介质分层半空间模型基础上的。当实际地层中包含孔隙介质层时,需要将孔隙介质简化为弹性介质来进行分析和处理。这种简化处理究竟对瑞利波勘探造成怎样的影响是本文所研究的问题。本文基于弹性介质与孔隙介质共同构成的分层半空间模型,首先推导了分层半空间中两种介质处于不同相对位置时的瑞利波频散方程,解决了不同阶数矩阵之间的变量传递问题然后,针对传递矩阵法在求解频散函数时可能出现的溢出问题,给出了一种可以有效提高计算范围的解决方案;同时,提出了一套新的数值算法用于复频散方程的快速求解。数值计算结果表明:当孔隙介质位于半空间表面时对低频瑞利波频散特性的影响最为显著,而在其它情况下对瑞利波频散特性的影响相对较小。
简介:地震反演是当今最广泛应用于含油气储层预测的技术之一,取得了很多很好的预测效果,但也有失败的例子,达不到区分岩性和识别流体的目的。而本文介绍的建立在岩石物理建模和分析基础上的地震弹性反演,可将含油气储层预测由定性向(半)定量推进一步。根据岩石物理建模和正演扰动分析,可深刻理解岩石物性参数与地震弹性参数之间的内在关系,进而寻找储层岩性、物性和含油气性的敏感地震弹性参数,建立起理论岩石物理解释图版。岩石物理分析结果和所建立的岩石物理解释图版分别用以指导地震反演和反演结果的解释,实现油气储层分布预测和流体检测的目的。文中的含气砂岩分布预测实例研究应用结果表明,这种方法较叠后地震反演储层预测技术具有无可比拟的优越性,效果更佳,效率更高。