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

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

简介：油层注入聚合物溶液后,电阻率测井响应发生与水驱不同的变化特征,严重影响水淹层测井评价的效果。本文首先分析了聚合物溶液的导电性质,其次按照聚合物驱油层不同施工流程开展岩电实验研究,分析不同驱替阶段的岩石电阻率的变化特征。并对比分析聚合物驱与水驱油层电阻率测井响应的差异。结果表明,不同聚合物驱替饱含油岩心过程,岩心电阻率随饱和度变化呈现"w"型或"S"型和单调下降的形态;清水聚合物驱油层电阻率测井响应呈现弱、未水淹的特征;污水聚合物驱油层电阻率测井响应呈现与污水水驱相似的变化特征。

简介：根据Ecker的水合物沉积物的三种微观模式,计算含水合物沉积层和含游离气沉积物的弹性模量,分析对比了水合物的不同微观模式、不同水合物饱和度以及不同游离气饱和度对沉积物弹性模量的影响;从纵横波分离的弹性波动方程出发,采用交错网格空间有限差分方法模拟地震波在海底天然气水合物沉积地层的传播,得到纵、横波的海底地震（OBS）共接收点道集。数值算例表明,当水合物作为流体的一部分或胶结颗粒骨架时,仅纵波记录上存在BSR;当水合物胶结颗粒接触,纵、横波记录上均存在BSR。并且,OBS会接收到上行纵波和上行横波在海底界面形成的转换波,干扰横波记录上BSR的识别。

简介：本文利用取自南海神狐海域钻井岩芯进行了核磁共振测量,研究了岩石中所含水合物对地层渗透率的影响,该岩芯样品为一＂孔隙充填型＂水合物样品。核磁共振测量结果表明岩石中所含水合物对地层渗透率影响很大,随着含水合物饱和度的增加,地层渗透率迅速降低。通过对测量结果的进一步分析可知＂孔隙充填型＂水合物对地层渗透率的影响-存在多方面的因素,包括水合物对地层孔隙度、孔隙迂曲度的影响以及对水合物对孔隙喉道的堵塞等。最后本文利用Masuda下降模型建立了适用于南海神狐海域的相对渗透率与含水合物饱和度的经验公式,下降指数N=7.9718。

简介：天然气水合物沉积层在常规叠后地震剖面上，一般表现为似海底反射（BSR）特征，常作为水合物识别的重要标志，但由于地震数据分辨率的限制与多解性的存在，类似的反射特征不一定是水合物的表现。本文先从南中国海神狐海域地质条件分析人手，综述研究区域天然气水合物发育的构造、沉积与运移环境；其次以加权方程与Zoeppritz公式的正演模拟为基础，讨论天然气水合物与BSR特征的关系；然后选取经过三口已钻获天然气水合物站位的地震测线，提取角度域共成像道集用于BSR的AVA特征研究；最后依据叠前地震弹性参数同步反演思路，对水合物重点赋存区域开展精细预测研究，获取了研究区域含水合物沉积层稳定带的分布特征与饱和度数。该方法可作为采用多道地震数据估算天然气水合物含量的一种有效技术。

简介：传统上，有限差分的差分系数一般可以通过泰勒级数展开法或优化方法来极小化频散误差得到。基于泰勒级数展开的差分法在有限的波数范围内精度较高，但在这个范围之外会产生较强的数值频散；基于最小二乘的优化有限差分法能在更大的波数范围内达到较高的精度，并可以在较小的计算需求内获得全局最优解。本文将基于最小二乘的优化有限差分法从二维正演模拟推广到三维，形成了计算效率高、高精度范围宽、适合并行计算的三维声波优化有限差分方法。频散分析及正演模拟表明本文发展的有限差分方法可以很好地压制数值频散。最后，将本文发展的有限差分方法应用到三维逆时偏移的震源波场延拓和检波点波场延拓中，并结合有效边界存储策略与checkpointing技术在GPU集群上实现三维逆时偏移以提高计算效率、减少存储量。三维逆时偏移试算结果表明本文三维优化有限差分方法与传统的有限差分法相比可以获得更高精度的偏移成像结果。

简介：本文提出了一种将高分辨率阵列侧向和方位电极系综合在一起的三维侧向测井电极系3D-LS,该电极系具有径向、纵向和周向探测能力。通过有限元数值模拟计算,考察了井眼尺寸、冲洗带电阻率、侵入深度、层厚及围岩电阻率对六种不同探测模式的影响,确定了电极系尺寸和探测特性。分析伪几何因子,低侵时电极系的探测深度最深可达1.5m,其值接近斯伦贝谢双侧向电极系深探测深度,而大于高分辨率方位侧向成像仪深探测深度,并且三维侧向测井电极系可提供多条径向不同深度曲线,可更好地描述地层侵入剖面。无限厚地层条件下,方位电极可识别出厚度0.1m的异常体,利用方位侧向曲线半幅点对应异常体厚度判断,对异常体纵向分层能力可达0.5m。高阻背景下,异常体的电阻率越低,越靠近井眼,方位越大于15度,越易被方位电极探测。数值模拟结果为后续三维侧向测井电极系的研究奠定了基础,对低阻异常评价具有一定的指导意义。

简介：海洋天然气水合物是一种潜在的巨大能源,而地球物理探测技术是勘探天然气水合物资源的重要方法。本文介绍了利用可控源电磁法探测海洋天然气水合物的可行性,在参考了大洋钻探计划(ODP)164航次对天然气水合物的采样数据的基础上建立了几种不同水合物含量的海底天然气水合物一维地电模型,并且利用建立的模型讨论了天然气水合物的频率域电磁响应。文章就天然气水合物的电场振幅值和相位随着收发距和频率等参数的变化关系等做了相应研究,同时就水合物厚度不同时其相应的电磁响应特征做了初步探讨,为海洋天然气水合物的勘探和资源评价提供了参考。

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

简介：为了进一步提高大地电磁三维快速松弛反演的计算效率，在深入分析大地电磁三维快速松弛反演算法的基础上，结合MPI自身的优越性，确定了并行计算的思路，实现了三维快速松弛反演的并行计算。通过理论模型合成数据和实测数据对实现的三维快速松弛反演并行程序进行了试算，分析对比了在多种情况下程序的执行效率。测试结果表明，所实现的三维快速松弛反演并行程序运行结果正确，效率提高明显。此成果在普通微机上实现，推动了大地电磁三维反演技术的实用化，可为其它地球物理三维正反演研究所借鉴。

简介：针对矩形网射线追踪存在的模型剖分灵活性差、速度界面描述精度差等问题,研究了复杂结构三角网最小走时射线追踪全局算法。（1）根据剖分区域点、线、面的结构关系,遵循Delaunay三角剖分的优化准则进行三角网格剖分;（2）定义三角单元射线追踪的拓扑关系;（3）波源点及某一时刻波到达的每一个节点点构成波行面,在波行面扩展过程中计算节点的最小走时和次级源位置,实际次级源检索采用双曲线近似算法;（4）利用各节点走时和次级源方向信息,通过最小走时搜索,拾取从接收点到源点的射线路径。数值模拟结果表明,三角网射线追踪方法模型剖分时灵活性强、速度间断面的描述精度高,追踪结果准确。

简介：高分辨率非线性三维整体反演方法是基于非线性理论,在层位控制下,将工区多井(或全部井)的测井数据与井旁地震道数据输入具有多输入多输出的网络,同时进行整体训练,可获得整个工区的自适应权函数,并建立综合非线性映射关系,并根据储层在纵横方向上的地质变化特征更新这种非线性映射关系,这样,就能对反演过程及其反演结果起到约束和控制的作用,从而获得稳定且分辨率高的地震反演剖面(速度反演剖面/波阻抗反演剖面/密度反演剖面),实现整体反演,该方法通过模型试算和实际资料处理,获得较好的地质效果,证明该方法精度高、实用性强,可用于储层的定量分析。

简介：井地电位成像是通过套管向井中供电或将电源放在井中，在地表观测电位异常的一项技术，其供电源有线源和点源两种类型。为了研究这两种电源对地下异常体产生的电位异常特征，本文针对不同激励源，采用有限差分方法进行数值模拟研究，在线性方程组求解电位时引入不完全Cholesky共轭梯度（ICCG）迭代方法，分别实现了点源和线源井地电位成像技术的三维正演。最后，基于阻尼最小二乘法实现了井地电位成像技术的电阻率三维反演。设计不同地电模型分别进行正演和反演试算，正演结果表明，供电电源的类型不同，异常体在地表的电位异常特征也不同；反演结果表明，低阻体的反演结果要好于高阻体，点源置于异常体下方时反演的电阻率对异常体边界的识别比线源更加准确。

简介：本文是三分量地震在陆相薄互层地震勘探中应用成功的首例。大庆长垣喇嘛甸油田处于高含水开发后期,储层预测的主要难题是密井网条件下厚度2m以上砂体的边界识别问题,常规单一纵波地震解释存在多解性,引入横波信息可增加地震解释的可靠性。本文根据纵、横波对气藏的响应特征,利用纵、横波联合振幅属性、分频属性对气藏进行了分析,能够很好地反应气藏的边界,体现了横波受流体影响较小的优势。通过测井曲线岩石物理定量的分析了该地区主要是密度、λρ和μρ三种物性参数能够很好地反映储层岩性变化。在以地震数据分频为核心的高精度地震参考标准层解释基础上利用纵、横波联合反演得出的纵、横波速度比、λρ、μρ等属性与测井解释的储层具有很好的一致性,多波属性λρ沿层切片描述了开发小层的砂体平面展布,与井砂岩图具有较高的一致性,刻画了井间砂体边界变化的细节,并指示了剩余油挖潜的有利区域。

简介：弹性参数在甜点区预测和页岩气的开发过程中扮演着重要的角色，因此研究等效弹性参数随页岩气储层属性的变化是一项很有意义的工作。研究中我们用x射线CT扫描技术获得了较为精确的页岩样品微观结构图像。从这些图像中，我们可以获得孔隙度和矿物的详细情况，据此，我们构建了三维数字岩心，并应用有限元法对弹性参数进行了数值模拟，其间深入考察了子样选取、网格划分、求解器类型以及边界条件等，该方法易于区别不同的矿物及其百分含量。本文重点研究孔隙度和干酪根含量对弹性参数的影响，计算结果表明，孔隙度和干酪根含量对弹性性质有较大的影响，当孔隙度和干酪根含量增加时，弹性模量降低，且当孔隙度小于0．75％左右、干酪根含量大于3％左右时弹性参数减小速率较缓。因为孔隙度仅仅为4．5％，孔隙中填充油或气对弹性参数的影响甚微。不同岩心样本具有不同的孔隙度和干酪根含量，传统岩石物理实验不仅昂贵而且费时，而数值模拟是基于数字岩心来计算弹性参数，更加经济、方便。本研究证实了将页岩样品的微观结构图像与弹性模量的计算相结合来预测页岩弹性参数的可行性。

简介：标量CSAMT只适合一维及测量方向与构造方向垂直的二维情况,对于复杂的三维地电结构,CSAMT需采用张量测量。本文试图采用矢量有限元法实现三维张量CSAMT的正演模拟。为了验证算法的正确性,本文在层状介质中计算了三维CSAMT远区的电场,磁场及阻抗张量,并且与层状介质中的理论解进行了比较,接着还模拟了均匀半空间中含有三维异常体的模型,并且分析了四个阻抗张量、视电阻率及阻抗相位的响应特征。得出如下结论：采用矢量有限元法来模拟三维张量CSAMT,其电磁场及阻抗张量的实虚部计算精度都比较高,并且该方法本身满足电场法向不连续,不用进行散度校正。

简介：地震波场数值模拟中不可避免地会出现边界反射,一般采用吸收边界条件以压制人工边界反射。目前常用的分裂式完全匹配层（PML）边界条件需要在边界处进行特殊处理,尤其是在三维情况下需要将变量分裂为三个分量,增加了数值模拟的计算时间和内存占用量。与分裂式PML吸收边界条件相比,混合吸收边界条件（HABC）具有易于实现、计算量小和吸收效果好等优点,可以提高三维波动方程数值模拟的计算效率。本文将基于一阶Higdon单程波方程的混合吸收边界条件从二维计算域发展到三维,提出了适用于三维弹性波数值模拟的混合吸收边界条件。均匀模型以及复杂模型的三维数值模拟结果表明,混合吸收边界条件与传统的完全匹配层边界条件相比,具有效率高、吸收效果好的优势。

简介：类中梯装置三维激电成像技术兼顾了激电测深与激电剖面二者的特点,能快捷高效的获得测区三维电阻率及极化率信息,空间分辨率高,探测深度大。本文从类中梯装置的数据采集方式入手,阐述了采用类中梯装置进行三维激电成像的工作方式.建立了一个三维地电模型,采用类中梯装置进行三维正演模拟,并进行反演及成像。数值模拟结果表明,采用类中梯装置进行三维激电成像能较好地刻画实际地电模型的特点。以甘肃省某多金属矿为例,利用阵列式电磁法综合测量系统,进行了三维大功率激电成像技术的应用研究。实际应用结果表明,采用三维激电成像技术,可以多角度、多细节显示测区地下介质电阻率、极化率的分布情况,明确电性突变界面的延展状态,快速圈定成矿有利区。该研究对多金属矿产勘探等领域具有一定的指导性意义。

简介：工区位于在大民屯凹陷荣胜堡-前进地区工业繁华地带,地表水系发育,表层岩性复杂,深层各种油气藏圈闭类型众多,勘探难度大,以往采集资料质量差.本文详细介绍了为改进资料质量所采取的一系列措施,其中包括科学的观测系统设计、认真细致的表层调查、高精度的测量,选择最佳岩性组合激发、低噪声接收、大型障碍区特观施工、定量化质量控制等工作.成功地进行了二次三维采集,实现了二次采集资料质量的飞跃,减少浅层资料缺失,保证了深层资料品质,保证了资料解释的准确性.文中展示了了实际效果.