简介:阻抗张量元素的计算是在大地电磁测深数据处理的重要一步。按照常规,阻抗张量被定义为以Zxx,Zxy,Zyx,和Zyy为元素的2×2矩阵。在本次研究中,6个元素的阻抗张量的计算使用了一个含有Zxx,Zxy,Zyx,Zyy,Zxz和Zyz分量的2x3矩阵。对上述两类阻抗张量元素的属性进行了分析。利用由印度古吉拉特邦卡奇沉积盆地采集的5个分量大地电磁数据测试了文中的方法。从视电阻率和相位的计算中我们是观察到在大部分的频带范围内4个元素阻抗和6个元素阻抗Zxy和Zyx两类元素区别不大。然而,较长周期时间的数据,如超过100秒,观察到视电阻率的增加和相位的减少。我们还注意到,倾子幅度在大部分时间几乎是零,但较长周期(超过100秒),逐渐呈增加的趋势。卡奇沉积盆地的地电断面表明在较长的周期内浅层近水平层和深层异常高电导性的不均质层都可能是引起大的Hz分量的原因。这表明,磁场垂直分量Hz对在大的2D/3D结构区域内的电场参数估计发挥的重要作用。
简介:标量CSAMT只适合一维及测量方向与构造方向垂直的二维情况,对于复杂的三维地电结构,CSAMT需采用张量测量。本文试图采用矢量有限元法实现三维张量CSAMT的正演模拟。为了验证算法的正确性,本文在层状介质中计算了三维CSAMT远区的电场,磁场及阻抗张量,并且与层状介质中的理论解进行了比较,接着还模拟了均匀半空间中含有三维异常体的模型,并且分析了四个阻抗张量、视电阻率及阻抗相位的响应特征。得出如下结论:采用矢量有限元法来模拟三维张量CSAMT,其电磁场及阻抗张量的实虚部计算精度都比较高,并且该方法本身满足电场法向不连续,不用进行散度校正。
简介:标量CSAMT只适合一维及测量方向与构造方向垂直的二维情况,对于复杂的三维地电结构,CSAMT需采用张量测量。本文试图采用矢量有限元法实现三维张量CSAMT的正演模拟。为了验证算法的正确性,本文在层状介质中计算了三维CSAMT远区的电场,磁场及阻抗张量,并且与层状介质中的理论解进行了比较,接着还模拟了均匀半空间中含有三维异常体的模型,并且分析了四个阻抗张量、视电阻率及阻抗相位的响应特征。得出如下结论:采用矢量有限元法来模拟三维张量CSAMT,其电磁场及阻抗张量的实虚部计算精度都比较高,并且该方法本身满足电场法向不连续,不用进行散度校正。
简介:Anewmulti-moderesistivityimagingsonde,withtoroidalcoilsassource,canconductthreeresistivitymeasurements:azimuthalresistivity,lateralresistivity,andbitresistivitymeasurements.Thus,theloggingtimeandcostaregreatlysaved.Thetoroidalcoilsaresimplifiedasanextendedvoltagedipoleandtheresponseequationsarederivedforahomogenousformation.Basedon3DFEM,thedepthofinvestigation(DOI),verticalresolution,circumferentialazimuthalcapacity,boreholediameter,mudresistivity,thicknessoftargetformation,andtheresistivityofthesurroundingformationandmudinvasionaresimulated.TheresultssuggestthatthethreemeasurementmodesofthenewsondearedifferentinverticalresolutionsandDOIs.Thecircumferentialdetectionabilityoftheazimuthbuttondependsonthecontrastbetweentheanomalyandformationresistivityandtheopenangleoftheanomaly.Whethertheboreholeistruncatedatthebitornothasagreatinfluenceonthesimulationresults.Theboreholeandmudinvasionaffecttheapparentresistivityinallmodes,buttheeffectsofresistivityofsurroundingformationandthicknessofthetargetformationareonlycorrectedforlateralresistivitymeasurement.
简介:在直流电法有限元数值模拟中,针对常规结构化网格源点附近网格节点数值精度低的问题,并考虑计算效率和反演成像的需求,本文在二维点源问题中提出一种新的网格加密一收缩方法。其核心思想是在结构化网格计算区域中先后引入网格节点加密和网格节点收缩两个环节。通过在计算区域水平方向上均匀加密网格节点密度,降低源点奇异性的影响范围,并提升对地形的模拟能力;通过在计算区域垂向方向上将多个网格收缩到一个网格中,降低网格节点的规模,进而提高数值计算的效率。理论模型检验表明,网格加密-收缩方法能有效地解决源的奇异性问题,与均匀加密网格相比,网格节点总数降低了约80%。
简介:本文基于弹性波动方程,从其弱形式出发,利用Galerkin变分原理,通过对方程进行空间和时间上的离散,在空间域中引入预条件共轭梯度的逐元算法,在时间域中引入时间积分的交错网格预处理/多次校正算法,发展了弹性波模拟的Chebyshev谱元算法。针对均匀固体介质和具有倾斜分层的分区均匀固体介质模型,通过与有限差分算法结果相比较验证其精度的可信性,同时利用该算法模拟了弹性波在具有水平分层的任意起伏自由表面模型中的传播,并分析了其传播特点。研究表明,我们提出的交错网格预处理/多次校正算法的Chebyshev谱元算法,保留了有限元法的优势,并且采用了具有最优张量乘积技术的元到元的算法,能够处理带有起伏自由表面的复杂介质模型,它具有比有限元法收敛快,计算效率较高等优点,特别适合于复杂结构和复杂介质中的弹性波传播的数值模拟。
简介:本文通过人工变换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模型更接近实际地层情况。采用区分填充模式的等效介质模型,能够更好地实现对祁连山冻土区水合物粉砂岩地层的模拟。
简介:含洞穴的碳酸盐岩地层具有强烈的非均质性及储集空间预测难度大的特点,利用随钻电阻率测井方法对井眼环境含洞穴的储层进行准确识别和划分,是当前研究的一个焦点问题。本文使用一种新型的高效和高精度自适应有限元方法(hp—FEM)模拟和分析了含洞穴地层随钻电阻率测井仪器响应。本文所提的hp—FEM与传统h-FEM相比,其结果具有网格自适应的特点,并且计算能够以指数速率收敛于较高的精度。数值实例使用自适应有限元方法研究地层中洞穴的大小、洞穴距离井眼的远近和仪器发射频率改变对测井响应的影响,并提供了识别含洞穴地层的方法。研究结果可以为实际测井中遇到的各种地层洞穴的准确识别和定量评价提供理论依据。
简介:根据Ecker的水合物沉积物的三种微观模式,计算含水合物沉积层和含游离气沉积物的弹性模量,分析对比了水合物的不同微观模式、不同水合物饱和度以及不同游离气饱和度对沉积物弹性模量的影响;从纵横波分离的弹性波动方程出发,采用交错网格空间有限差分方法模拟地震波在海底天然气水合物沉积地层的传播,得到纵、横波的海底地震(OBS)共接收点道集。数值算例表明,当水合物作为流体的一部分或胶结颗粒骨架时,仅纵波记录上存在BSR;当水合物胶结颗粒接触,纵、横波记录上均存在BSR。并且,OBS会接收到上行纵波和上行横波在海底界面形成的转换波,干扰横波记录上BSR的识别。
简介:天然气水合物沉积层在常规叠后地震剖面上,一般表现为似海底反射(BSR)特征,常作为水合物识别的重要标志,但由于地震数据分辨率的限制与多解性的存在,类似的反射特征不一定是水合物的表现。本文先从南中国海神狐海域地质条件分析人手,综述研究区域天然气水合物发育的构造、沉积与运移环境;其次以加权方程与Zoeppritz公式的正演模拟为基础,讨论天然气水合物与BSR特征的关系;然后选取经过三口已钻获天然气水合物站位的地震测线,提取角度域共成像道集用于BSR的AVA特征研究;最后依据叠前地震弹性参数同步反演思路,对水合物重点赋存区域开展精细预测研究,获取了研究区域含水合物沉积层稳定带的分布特征与饱和度数。该方法可作为采用多道地震数据估算天然气水合物含量的一种有效技术。
简介:传统上,有限差分的差分系数一般可以通过泰勒级数展开法或优化方法来极小化频散误差得到。基于泰勒级数展开的差分法在有限的波数范围内精度较高,但在这个范围之外会产生较强的数值频散;基于最小二乘的优化有限差分法能在更大的波数范围内达到较高的精度,并可以在较小的计算需求内获得全局最优解。本文将基于最小二乘的优化有限差分法从二维正演模拟推广到三维,形成了计算效率高、高精度范围宽、适合并行计算的三维声波优化有限差分方法。频散分析及正演模拟表明本文发展的有限差分方法可以很好地压制数值频散。最后,将本文发展的有限差分方法应用到三维逆时偏移的震源波场延拓和检波点波场延拓中,并结合有效边界存储策略与checkpointing技术在GPU集群上实现三维逆时偏移以提高计算效率、减少存储量。三维逆时偏移试算结果表明本文三维优化有限差分方法与传统的有限差分法相比可以获得更高精度的偏移成像结果。
简介:本文提出了一种将高分辨率阵列侧向和方位电极系综合在一起的三维侧向测井电极系3D-LS,该电极系具有径向、纵向和周向探测能力。通过有限元数值模拟计算,考察了井眼尺寸、冲洗带电阻率、侵入深度、层厚及围岩电阻率对六种不同探测模式的影响,确定了电极系尺寸和探测特性。分析伪几何因子,低侵时电极系的探测深度最深可达1.5m,其值接近斯伦贝谢双侧向电极系深探测深度,而大于高分辨率方位侧向成像仪深探测深度,并且三维侧向测井电极系可提供多条径向不同深度曲线,可更好地描述地层侵入剖面。无限厚地层条件下,方位电极可识别出厚度0.1m的异常体,利用方位侧向曲线半幅点对应异常体厚度判断,对异常体纵向分层能力可达0.5m。高阻背景下,异常体的电阻率越低,越靠近井眼,方位越大于15度,越易被方位电极探测。数值模拟结果为后续三维侧向测井电极系的研究奠定了基础,对低阻异常评价具有一定的指导意义。
简介:波动方程有限差分法是地震数值模拟中的一种重要的方法,对理解和分析地震传播规律、分析地震属性和解释地震资料有着非常重要的意义。但是有限差分法由于其离散化的思想,产生了不稳定性。精细积分法在有限差分法的基础上,在时间域采用解析解的表达形式,在空间域保留任意差分格式,发展成为半解析的数值方法。本文结合并发展了以往学者的成果,推导了任意精细积分法的三维弹性波正演模拟计算公式,并对其稳定性进行了数值分析。在计算实例中,实现了精细积分法二维和三维弹性波模型的地震正演模拟,对计算结果的分析表明,精细积分法反射信号走时准确,稳定性好,弹性波场相较于声波波场,弹性波波场成分更为丰富,包含了更多波型成分(PP-和PS-反射波、透射波和绕射波),这对实际地震资料的解释和储层分析有重要的意义。实践证明,该方法可直接应用到弹性波的地质模型的数值模拟中。
简介:针对矩形网射线追踪存在的模型剖分灵活性差、速度界面描述精度差等问题,研究了复杂结构三角网最小走时射线追踪全局算法。(1)根据剖分区域点、线、面的结构关系,遵循Delaunay三角剖分的优化准则进行三角网格剖分;(2)定义三角单元射线追踪的拓扑关系;(3)波源点及某一时刻波到达的每一个节点点构成波行面,在波行面扩展过程中计算节点的最小走时和次级源位置,实际次级源检索采用双曲线近似算法;(4)利用各节点走时和次级源方向信息,通过最小走时搜索,拾取从接收点到源点的射线路径。数值模拟结果表明,三角网射线追踪方法模型剖分时灵活性强、速度间断面的描述精度高,追踪结果准确。