简介：本文计算了地面磁电阻率法在监测地下水体中垃圾渗漏液运移的动态磁电阻率异常。基于对垃圾场地概念模型的水文地质仿真,建立了数值地电模型,并采用改进的有限差分软件MMR2DFD,分别计算了点源,MMR-TE方式,MMR-TM方式下半空间模型以及断层模型的磁电阻率异常响应。半空间模型的正演结果表明,垂直于主要电流方向的磁场分量包含更多地下异常信息,MMR-TE模式能有效进行地下水污染监测。断层模型的计算结果更进一步验证了MMR在地下水污染物运移监测中的有效性。

简介：基于岩石物理模型和广义各向异性Zoeppritz方法在频率域计算裂缝型反射体模型反射波方位地震AVO响应。反射体模型为粘弹各向异性、有限厚度的地质体，其地震反射波形序列包含如下动力学信息，即分界面处介质的波阻抗和非弹性差异、反射体内部波的各向异性传播、在传播路径上的频散与衰减，以及来自顶底界面的反射波的调谐与干涉等。计算表明，速度频散和衰减增顶界面反射波大入射角反射时的振幅，而减弱底界面反射振幅。对于固定人射角的方位地震方位地震响应，PP波反射特征表现为随方位角的增加反射波形序列延续时间变长，而PSV和PSH转换类型反射波的方位各向异性变化特征稳定且受储层厚度影响较小，表现为PSV波反射振幅随方位角增加而增加，PSH波在0°和90°方位无反射能量，在45°方位反射振幅最强。

简介：致密储层与围岩波阻抗差异小，常规储层反演方法难以有效预测出有利区。中国西部致密储层非均质性强，常规储层反演预测难度大，本文提出了基于Xu-White模型的叠前储层预测方法。首先修正Xu-White模型的试验参数，直至该模型计算出来的纵、横波速度及密度与实际测井得到的数值基本吻合，然后用修正的模型计算不同储层物性、岩性、含油性下的纵、横波速度、密度及泊松比，进而建立不同的砂泥岩互层地质模型。其次，采用Zoeppritz方程近似式进行叠前正演，分析不同储层物性、岩性、含油性时的叠前道集响应特征，在对中国西部Z区实际叠前道集优化处理的基础上，选择其中1种对岩性、含油性、物性最敏感的弹性参数，对其有利储层分布区进行预测，其预测结果具有较高的精度。

简介：常规地震勘探大多基于反射波理论,对尺度小、倾角陡、构造复杂的灾害性地质体,很难得到精确的成像结果。本文建立2类典型矿井地质模型,采用有限差分正演模拟算法,研究了矿井地震散射波的特征,将等效偏移距偏移（EOM）和基于等效偏移距干涉叠加偏移的散射成像方法用于矿井地震勘探中,重点分析了散射波成像相比常规反射波地震成像的技术优势。研究表明：1）散射成像方法能够提高覆盖次数,充分利用有效的弱散射信号,对不均匀复杂地质体的散射波的成像效果明显优于反射波叠后偏移;2）对于巷道超前探测,它弥补了基于射线理论的叠加偏移的不足,为提高矿井地震勘探分辨率提供了有效的成像方法。

简介：针对祁连山冻土区DK-4井孔含水合物岩层构建岩石物理模型，分别采用K-T方程模型方法和区分填充模式的等效介质模型方法（模式I和模式II）。K-T方程主要模拟地震波在两相介质中传播，基于弹性模量计算速度；而等效介质模型主要依据对水合物地层介质的两种假设：模式I是将水合物作为孔隙填充物的一部分，模式II是将水合物作为岩石骨架的一部分。首先根据粉砂岩层段的测井数据提取了水合物地层骨架的物性参数，包括纵波速度、横波速度、密度、体积模量和剪切模量。然后依据水合物地层各主要成分的物性参数，建立了基于K-T方程的岩石物理模型和区分填充模式的等效介质模型。将两类模型的速度曲线分别与实际地层数据进行了对比：由K-T方程建立的岩石物理模型，其理论计算的速度偏离实际值；而区分填充模式的等效介质模型，其理论计算的速度符合实际值，并且填充模式II模型的速度曲线比填充模式I模型更接近实际地层情况。采用区分填充模式的等效介质模型，能够更好地实现对祁连山冻土区水合物粉砂岩地层的模拟。

简介：本文利用x射线CT获取反映岩心微观结构的三维数字岩心，利用数学形态学中的开运算模拟了岩石的油驱水排驱过程中，不同含水饱和度下油和地层水在孔隙空间中的分布。利用有限元方法计算了岩石电阻率，进而得到岩石地层因素和电阻率指数，并考查了岩石润湿性对岩石电阻率指数的影响。数值模拟结果表明：基于数字岩心的水湿岩石地层因素和电阻率指数数值模拟结果与实验结果一致，拓展了岩石电阻率实验的能力；岩石润湿性对岩石电性有重要影响，在相同含水饱和度下，油湿岩石电阻率高于水湿岩石电阻率，油湿岩石饱和度指数远大于水湿岩石饱和度指数。

简介：Webuiltathree-dimensionalmodeltosimulatethedisturbanceofthestressfieldnearthereversefaultinZhaziao,LeyiTownshipowingtohydraulicfracturing.Theporepressure,andshearandnormalstressesduringfracturingareanalyzedindetail.Inputrockmechanicsparametersaretakenfromlaboratorytestdataofshalesamplesfromthestudyarea.Thesimulationresultssuggestthatafter16hoursoffluidinjection,thepore-pressurevariationcanactivatethereversefault,i.e.,weobservereverseslip,andtheshearstressanddisplacementonthefaultplaneincreasewithtime.Thebiggeststress–strainchangeoccursafteronehouroffluidinjectionandtheyieldpointappearsabout0.5hafterinjection.Toobservethestressevolutionineachsection,thenormaldisplacementontheboundaryisconstrainedandthefaultplaneissetasnonpermeable.Thus,theslidingislimitedandthesheardisplacementisonlyinthescaleofmillimeters,andthecalculatedmagnitudeoftheinducedearthquakesisbetweenMw-3.5andMw-0.2.Thesimulationresultssuggestthatfluidwaterinjectionresultsininhomogeneousfracturing.Themainrupturedareasarearoundtheinjectionpositions,whereastheextentofrupturingandcracksinotherareasarerelativelysmall.Nevertheless,nonnegligiblefaultactivationisrecorded.Sensitivityanalysisofthekeyparameterssuggeststhattheporepressureismostsensitivetothemaximumunbalancedforceandtheinternalfrictionanglestronglyaffectsthefaultslip.Finally,thecomparisonbetweentheeffectivenormalstressandthemaximumandminimumprincipalstressesonthefaultplaneexplainsthefaultinstability,i.e.,theMohrcirclemovestowardstheleftwithdecreasingradiusreducesandintersectsthecriticalslipenvelope,andcausesthefaulttoslip.

简介：地震波场数值模拟方法对理解和分析地震波的传播规律具有着重要的意义。弹性波动方程能够模拟地下介质的实际情况，为偏移和成像提供有效的依据。在弹性波波场数值模拟中，旋转交错网格数值模拟（RSM）修改了标准交错网格数值模拟（SSM）方法，将同类的参数定义在同样的节点上，拓宽了稳定性条件的约束，但在低速区会出现较严重的频散。变阶数差分方法是自适应空间算子长度方法的一种变化和推广。它以理论频散误差研究为基础，结合实际波场传播的情况进行误差计算，对不同速度匹配不同的差分阶数。本文研究了变阶数旋转交错网格数值模拟（VRSM），即是籽变阶数方法应用到RSM中，它可以很好地解决RSM在低速区域的数值频散问题，以及减少不必要的时间损耗；同时讨论了旋转交错网格的理论频散特性，并基于波场分离的方法分析了实际波场传播的频散误差，将原方法的应用范围由声波推广到剪切波，由理论值推广到时变值。在数值模拟试验中，VRSM将被应用于水平层模型和Overthrust模型。通过阶数分配以及相应波场传播效果和计算时间的分析，验证了该方法应用于复杂介质波场模拟中的实用性和有效性。实验的结果表明VRSM能够合理分配不同速度所对应的差分阶数，能保证计算的精确性，并合理控制计算的时间。

简介：波动方程有限差分法是地震数值模拟中的一种重要的方法，对理解和分析地震传播规律、分析地震属性和解释地震资料有着非常重要的意义。但是有限差分法由于其离散化的思想，产生了不稳定性。精细积分法在有限差分法的基础上，在时间域采用解析解的表达形式，在空间域保留任意差分格式，发展成为半解析的数值方法。本文结合并发展了以往学者的成果，推导了任意精细积分法的三维弹性波正演模拟计算公式，并对其稳定性进行了数值分析。在计算实例中，实现了精细积分法二维和三维弹性波模型的地震正演模拟，对计算结果的分析表明，精细积分法反射信号走时准确，稳定性好，弹性波场相较于声波波场，弹性波波场成分更为丰富，包含了更多波型成分（PP-和PS-反射波、透射波和绕射波），这对实际地震资料的解释和储层分析有重要的意义。实践证明，该方法可直接应用到弹性波的地质模型的数值模拟中。

简介：通过引入岩石骨架参数表达式到Gassmann近似流体方程中，实现了基于该方程的地震孔隙度反演。但是岩石骨架参数表达式有很多种，使得地震孔隙度反演公式也各不相同，实际应用不方便，反演结果也难以评价。作者在研究现有常用岩石骨架参数表达式如Esheby－Walsh，Pride，Geertsma，Nur，Keys－Xu以及Krief等的基础上，提出了含有调节参数的岩石骨架模型统一表示式，发展了基于Gassmann方程的地震孔隙度反演方法。该方法应用范围宽，参数调节方便、灵活。为验证所提公式和方法的实用性，作者们结合ZJ地区钻井岩心的岩石物理样品测试数据和测井数据，进行了地震孔隙度反演。由于该地区油气储层的存在与中等程度孔隙度的分布有关，因此地震孔隙度反演对预测和识别油气层、干层和含水层有重要作用。反演结果与工区内孔隙度测井瞌线吻合度高，说明本文所提公式与方法是有效的，可靠的。

简介：非裂化的完全匹配层(NPML)吸收边界条件首先用于有限差分方法模拟固体中弹性波的传播。本文将NPML加以推广用来模拟孔隙介质中弹性波传播。采用空间四阶精度和时间二阶精度中心差分将Biot方程离散化到交错网格。我们选取柱状双层模型来验证NPML有效性。并将NPML应用到储层环境下裸眼井模型,数值结果和离散波数法得到结果进行了对比,符合较好。

简介：本文在Oristaglio等(1984)和Adhidjaja等(1985)工作基础上,给出线源二维时间域瞬变电磁二次场的DuFort-Frankel有限差分数值解,有效避免了在总场求解法中场源附近的奇异问题,并对地-空边界电导率的处理、归一化感应电动势偏导数的计算、推进时间步的确定,提出了改进方法;吸取前人成就中二次场地-空边界向上延拓和零值边界处理技术,从而简化了计算方法;通过对均匀大地、水平层状大地模型的计算,二次场求解法与解析法的最大相对误差小于0.01%,计算速度比总场求解法提高了约3倍;模拟计算不同时刻瞬变电磁场在地下的分布形态,描绘出感应涡流向下向外的传播特征,以及与地下异常体相互作用的物理过程。

简介：目前观测系统设计中接收道数的选择基本上是基于水平层状介质假设，这显然不能满足复杂地区采集的需求。本文在朱金平研究基础上，通过野外试验法获得最优接收道数问题，进行数值模拟实现，在固定道间距情况下，建立接收总能量与平均效率能量之间的数学模型，获得最优接收道数，并利用一个复杂工区的地震数据对方法的正确性进行验证，模型数据计算的结果与实际数据处理的结果表明，本文方法与实际处理结论一致，说明该方法具有较高的科学性，可以为复杂地区观测系统设计提供新的思路。

简介：频率空间域地震波数值模拟具有独特的优势:可以同时模拟多源的波传播、每个频率之间独立并行地计算、计算频带选择灵活、不存在累计误差、容易模拟粘弹性介质中地震波传播。但是该方法的最大瓶颈是对于计算机内存的巨大需求。我们使用压缩存储系数矩阵的方法,极大地减少了计算机内存的需求量。同时为了减少短差分算子的数值频散,引用了频率空间域25点弹性波波动方程的差分格式,并使用了最小二乘意义下求出的优化差分系数。为了克服边界反射,采用了最佳匹配层吸收边界条件。数值模拟试验证明:用压缩存储系数矩阵及优化差分系数的频率空间域25点差分格式进行弹性波正演模拟,可以减少数值频散,提高计算精度。使用较大的网格间距,降低计算机内存需求,并保持较高的计算效率。该正演方法为后续弹性波偏移和弹性参数反演提供较好的基础。

简介：许多地震属性,比如地震的纵横波速度,泊松比,声波阻抗,都能用来区分不同的岩石类型,利用拉梅阻抗作为进一步的约束来代替声波阻抗能更好的降低结果的非唯一性。我们用这种方法构建了渤海湾地区一条剖面的岩石成分模型。我们简要介绍了区分岩石类型的地震参数,重点在于不同的参数组合对于分辨岩石类型的敏感性。为了与室温常压条件下测得的有代表性的岩石的地震波速度和密度进行对比,我们首先需要做温度和压力的校正。接着,我们需要用已知资料和实验测得的各种岩石样品做对比,来确定研究区的岩石类型。我们收集了渤中凹陷一条剖面的纵波速度资料。该剖面10km深范围内的属于同种岩石类型的不同岩石最终通过地震波速度、泊松比、密度、声波阻抗、以及拉梅阻抗综合约束。更多还原

简介：在分析Ziolkowski气泡振动模型局限性的基础上,建立了多种实际因素影响下的海上地震勘探单枪子波模型。此模型考虑了气泡壁的热传导作用、枪口节流作用、气泡上浮、液体粘度和枪体本身等对气枪子波的影响。相对于Ziolkowski模型,改进模型所模拟的气枪子波主峰值减小,气泡振动衰减加快,与实测子波吻合性较好。实验分析表明：（1）枪口节流作用控制着气枪脉冲峰值振幅的大小,（2）上浮过程中气泡周围静水压力值减小,气泡振动的周期随之改变,（3）热传导作用和流体粘度是引起气泡振动衰减的主要因素。

简介：含洞穴的碳酸盐岩地层具有强烈的非均质性及储集空间预测难度大的特点，利用随钻电阻率测井方法对井眼环境含洞穴的储层进行准确识别和划分，是当前研究的一个焦点问题。本文使用一种新型的高效和高精度自适应有限元方法（hp—FEM）模拟和分析了含洞穴地层随钻电阻率测井仪器响应。本文所提的hp—FEM与传统h-FEM相比，其结果具有网格自适应的特点，并且计算能够以指数速率收敛于较高的精度。数值实例使用自适应有限元方法研究地层中洞穴的大小、洞穴距离井眼的远近和仪器发射频率改变对测井响应的影响，并提供了识别含洞穴地层的方法。研究结果可以为实际测井中遇到的各种地层洞穴的准确识别和定量评价提供理论依据。

简介：Finite-difference(FD)methodsarewidelyusedinseismicforwardmodelingowingtotheircomputationalefficiencybutarenotreadilyapplicabletoirregulartopographies.Thus,severalFDmethodsbasedonthetransformationtocurvilinearcoordinatesusingbody-fittedgridshavebeenproposed,e.g.,standstaggeredgrid(SSG)withinterpolation,nonstaggeredgrid,rotatedstaggeredgrid(RSG),andfullystaggered.TheFDbasedontheRSGissomewhatsuperiortoothersbecauseitsatisfiesthespatialdistributionofthewaveequationwithoutadditionalmemoryandcomputationalrequirements;furthermore,itissimplertoimplement.WeusetheRSGFDmethodtotransformthefirstorderstress–velocityequationinthecurvilinearcoordinatessystemandintroducethehighprecisionadaptive,unilateralmimeticfinite-difference(UMFD)methodtoprocessthefreeboundaryconditionsofanirregularsurface.Thenumericalresultssuggestthattheprecisionofthesolutionishigherthanthatofthevacuumformalism.Whentheminimumwavelengthislow,UMFDavoidsthesurfacewavedispersion.WecompareFDmethodsbasedonRSG,SEM,andnonstaggeredgridandinferthatallsimulationresultsareconsistentbutthecomputationalefficiencyoftheRSGFDmethodishigherthantherest.

简介：本文通过人工变换T2分布和建立管-球模型模拟法研究含水合物地层渗透率与水合物含量之间的关系。首先,在渗透率的模拟试验中,我们改变了束缚水与可动水的比例、总孔隙度以及与之关联的T2分布。试验结果表明,相对渗透率与水合物含量之间的关系受到这些因素的制约。随后,我们用管-球模型表示水合物生长的孔隙空间,并把水合物的生长过程看成是向孔隙空间随机扔小球的过程。在此过程中,采用两种方法计算渗透率,一是Schlumberger＇sT2公式（即SDR模型）,二是Darcy定律与Poiseuille流动方程相结合的方法。前人的实验研究表明,在一定的水合物含量范围内,渗透率基本保持不变。以此为参考,我们将计算结果与之进行比较。我们发现,采用SDR模型时,渗透率的数值模拟曲线与Masuda模型N=15时的结果相近。而采用Darcy定律时,渗透率模拟值较高,但与实验结果的趋势相一致,都会出现渗透率的平直阶段。尤其,当水合物晶体在孔隙体内优先生成时,优先的概率越高,渗透率的平直范围越大。

简介：提出一种对整幅高光谱图像的稀疏表示结果进行直接显示的方法,图中不仅包含了稀疏表示中保留的光谱信息,还可显示整幅图像的空间信息。稀疏表示后,将字典中的各有效原子根据光谱特性选择颜色标签,之后根据稀疏系数进行混合颜色显示,此时的图像能够同时满足可分性及距离保持特性。针对局部地物时,提出的单像素混合阵列表示法及改进的裂片纹理技术能够直观且完整的显示出每个像元的具体组成情况,还能够根据所生成图像中的信息对原始HSI进行重建,进而提高数据的利用率。该模型不仅能够良好地显示地物的空间特性,同时能够显示稀疏系数的组成,同时单像素混合阵列表示法及裂片纹理技术弥补了混合像素彩色显示中颜色表达混乱的弊端。对真实地物数据进行实验,结果证明该模型产生的彩色图像具有良好的视觉效果及可分性,满足距离保持特性。