简介:描述和估计介质含水量、介电常数等属性参数分布是探地雷达探测技术的重要研究内容。雷达波的旅行时间和反射振幅系数与介质含水量、孔隙度与介电常数密切相关。常规通过旅行时间计算波速以估计介质参数的方法,例如透射波法,共中心点速度分析等,对于复杂介质分辨率有限。基于反射振幅的阻抗反演方法可以直接根据反射系数计算雷达波阻抗以估计介质属性参数,从而有效地避开常规方法在计算波速时精度低的问题。本文首先建立了基于高斯型和指数型混合自相关函数的三维多尺度等效随机介质模型刻画地下随机介质参数分布,并在局部加入高斯椭圆方程描述局部随机异常目标体。其次,通过引入锥形函数以降低随机介质模型在离散网格数值计算方法误差。在此基础上,推导了探地雷达递推阻抗反演的基本流程并结合随机介质模型测试了该方法在复杂介质参数估计中的计算精度。最后,对内蒙地区的实测探地雷达数据利用递推阻抗反演方法来估计地下污染物参数,估计介电常数、含水量结果与钻孔实测数据和同期开展的电阻率成像结果有很好的吻合。说明基于递推阻抗反演方法在探地雷达复杂介质属性参数估计中具有很好的应用前景。
简介:本文以中观孔隙结构的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分量地震资料计算各向异性参数.为此提出了两种工作流程:双扫描和单扫描流程.理论数据和实际数据的应用表明这两种流程得出的结果其精度相似,但是单扫描流程比双扫描更有效.
简介:常规声波测井解释方法建立在各向同性地层中的垂直井基础之上,对各向异性地层中的水平井和斜度井并不适用。在水平井和斜度井的钻井过程中,岩屑会溅到井壁并最终在井壁底部形成一层薄岩屑层。同时,高速层和各向异性也可能会影响声波时差的测量。本文采用结合了Hybrid-PML的交错网格时域有限差分方法(SGS-FDTD),研究阵列声波仪器在各向异性介质中的水平井和斜度井条件下,不同因素对测得时差所带来的影响。模拟结果揭示了声波时差在高速层,薄岩屑层,井斜角,层厚,各向异性等不同影响因素下的相应规律。我们发现当波长远小于井壁与高速层的距离时,能够准确测得目的层时差;然而在两者处于同一量级时的近场条件下,几何声学理论不再适用;当井眼底部存在岩屑层时,费马原理依旧适用,对真实时差测量没有影响;在各向异性地层环境下,井斜角的增大使得测得时差由垂直向时差逐渐趋近于水平向时差;对于特定源距的声系,仪器在目的层中所移动的距离超过一定距离时,方可获得薄层时差值。基于不同井斜角和不同各向异性大小均匀TI地层模型的模拟结果,绘制了时差估计图版,可以对任意井斜角和各向异性比的模型进行时差的定量估计。通过对不同声系和不同弹性参数的模拟实验,证明了这种时差估计方法对TI地层中的水平井和斜度井准确有效。
简介:边界识别是重磁数据解释中的常用方法之一,依据其结果可划分出地质体的水平范围。边界识别结果受地质体埋深及导数计算误差的影响所识别边界与真实边界之间存在一定的差距,且边界识别法无法直观地给出地质体的深度信息。为了获得异常体的水平位置和深度信息,本文提出空间归一化边界识别方法,其对不同深度的边界识别函数进行归一化计算,空间归一化边界识别法的最大值对应于异常体的水平位置和深度。常规边界识别结果的误差随理深的减小而减小,而空间归一化边界识别法是通过最大值来判断地质体的位置,最大值是在地质体处获得,因此归一化边界识别方法所获得的结果是准确的。通过理论模型试验证明归一化边界识别方法能有效地完成异常体的水平位置和深度的计算,所获得的水平位置和深度信息与理论值相一致,为下一步的勘探计划提供了更加可靠的依据。将其应用于实际航磁数据的解释,获得了断裂的具体分布形式。更多还原