木里地区水合物钻孔井壁构造裂缝特征

木里地区三露天井田断层发育、岩性复杂,且钻探结果揭示水合物主要赋存于地层裂缝中,因此,笔者专门利用声成像测井资料对三露天井田水合物钻孔井壁构造裂缝进行研究,明确了研究区构造裂缝的成像测井响应特征,并从裂缝线密度、倾角、倾向等方面对裂缝的分布特征进行了分析讨论.结果 表明:①研究区水合物钻孔井壁裂缝发育,倾角以高角度为主,低角度次之,水平裂缝与高角度裂缝少见;冻土层及以上裂缝较冻土层以下更为发育.②研究区裂缝主要走向与地层断裂走向基本一致或小角度斜交,裂缝的发育受断层控制明显.③研究区裂缝倾角相对较大,不利于烃类气体的大量聚集,这可能是该区域钻孔少见水合物的原因之一.     

基于电场不连续边界条件的层状介质电磁格林函数计算

给出了适用于全空间层状模型偶极子源格林函数的去奇异性广义描述技术。过源平行于层界面设置一个虚界面,使相邻层界面限制的区域中电磁矢量势函数可以分解为TE模式和TM模式标量势函的叠加,且它们都满足齐次Helmhlotz方程,相应的解可以通过叠加上行波和下行波函数得到。源的奇异性被耦合到切向电磁场所满足的不连续边界条件上,并通过反射系数所描述的振幅递推关系延拓到其他任意层中,实现了全空间任意位置格林函数的计算。该思路被应用到VED源激发的海洋电磁模型数值模拟中,计算结果很好地反映了电磁场在海洋模型中的传播过程,验证了该技术思路的正确性。     

海洋可控源电磁法发射源姿态影响研究

海洋电磁理论是建立在水平发射偶极的基础上.然而,由于海底洋流作用海洋可控源发射偶极可能产生水平摆动、倾斜和水平旋转,造成测量数据与理想状态的水平电偶极子的电磁响应有较大的偏差.为了研究发射源姿态变化对场的影响,本文提出利用欧拉旋转将发射偶极由源坐标系转换到地球坐标系,进而分析发射源姿态变化造成海洋电磁观测数据的误差分布...     

地面可控源频率测深三维非线性共轭梯度反演

讨论了地面可控源电磁勘探三维非线性共轭梯度反演的可行性以及反演过程中考虑场源的必要性.反演采用非线性共轭梯度反演方法.反演过程中,模型响应利用交错网格有限差分技术计算.反演数据采用与发射源平行的电场x分量Ex.利用层状导电模型作为背景,设计了两个理论模型进行数值试验:第一个模型中包含两个电阻率异常,以检验反演的有效性;第二个模型中,在测区外设置了一个低阻异常,以考察源的信息在反演中的作用.两个模型的反演分别从层状背景模型开始,迭代120次后终止.数值试验结果表明,(1)非线性共轭梯度反演所获得的电阻率分布和理论模型吻合较好;(2)非线性共轭梯度算法收敛速度较慢,需要较多的迭代次数完成反演;(3)对于可控源频率电磁勘探,必须考虑源位置信息.因此,本文采用考虑场源信息的地面可控源非线性共轭梯度反演方法能完成真正意义上的可控源频率电磁测深数据的反演.     

基于二级近似离散复镜像法的低频电磁格林函数计算

本文讨论了利用二级近似离散复镜像法实现低频格林函数快速、精确的计算.通过数值分析,给出了在低频电磁场计算中该方法近似参数的选取原则.(1)二级近似的积分区间里谱格林函数采样个数N为十倍的近似多项式的个数M;(2)总的积分区间L2大小约为40/r,r为收发距离;(3)两个积分区间的分界值L1为总的积分区间大小L2与第二个积分区间采样数N1的比值的Q倍;(4)Q的最优值根据谱格林函数拟合相对误差随Q值变化的近似误差曲线与横轴的交点来确定.将该参数选取原则应用到有限长双极源电磁场计算中,计算得到的视电阻率和直接积分结果吻合,且具有更高的精度,表明本文提出的计算方案是可行的.     

地质-地球物理三维可视化建模及其应用——以雄安新区为例

地球物理方法作为探测地质结构的有效手段,可提供丰富的地质资料,利用地球物理数据开展三维可视化建模,能够将其转化为有用的地质信息,对形象、直观、准确地认识地下情况具有重要意义。针对雄安新区三维地质结构探测需求,基于钻孔、地质剖面、地震、测井等数据资料,采用三维构造建模、属性模拟、三维数据体网格分析等技术,建立了雄安新区三维地质结构模型和容东地区浅部三维物性模型,实现了三维模型的可视化显示与分析。地质-地球物理三维可视化建模为查明雄安新区地层结构和构造格架,评价水源地、深层地热资源和工程场地等提供了技术支撑。     

同时考虑退磁和剩磁的有限体积法正演模拟

为分析同时考虑退磁和剩磁对磁测数据解释的影响,探讨了利用有限体积法求静磁场数值解的方法.从静磁场中的麦克斯韦方程出发,导出了有限体积法控制方程的离散表达式,对边界条件近似处理后求解方程组得到磁异常.通过与退磁改正计算结果对比,验证了方法的正确性,并分析得到忽略剩磁的相对误差与科尼斯布格比（Q）相关;利用有限体积法计算长方体模型在无地磁场情况下的磁异常和内部磁化强度,从数值模拟上说明剩磁也需要进行退磁改正,并表明退磁作用对剩磁的影响不仅与磁化率相关,而且与剩磁的方向和磁性体的形态相关;组合模型的计算结果对比表明,退磁作用对剩磁的影响还会因为临近强磁性体的作用发生改变.在青海灶火河西工区的应用说明,开展同时考虑退磁和剩磁的解释方法对准确识别强磁性岩体具有实用价值.

     

矩形大定源层状模型瞬变电磁响应计算

针对矩形大定源介绍了一种层状导电模型瞬变电磁响应的快速计算方法.其基本思想是,将大定源回线源用较少的有限大小方形回线叠加.在频率域,该叠加过程导出矩形大定源回线模型响应可以表示为一个标准的Hankel积分,叠加效应仅仅是相应Bessel函数的空间积分过程.瞬变电磁响应可以利用余弦变换由频率域响应获得.均匀半空间模型计算结果和解析解较为吻合,层状模型模拟结果符合物理规律,不同测点位置感应电动势的视电阻率具有很好的一致性.新方法误差只是来源于圆回线对小方形回线的等效误差.与均匀半空间解析结果比较,只用16个方形回线等效,相对误差就可以小于10-3.     

海洋环境中水平电偶极子电磁场特征分析

海洋环境中人工源电磁信号可以用来进行海洋资源勘查、钻前储层评价以及海洋环境中的通信,具有重要的研究意义.为了研究可控源装置在海洋环境中产生的电磁响应特征,本文针对典型海洋层状导电模型,利用由虚界面法获得的电磁场表达式,结合离散复镜像法实现了长导线源产生的电场和磁场的计算.通过数值计算结果可以看出:1)不同的场分量的空间分布规律不同;2)不同方向的Ex分量相对Hy分量受空气波的影响较小,故海洋电磁法勘探应选用Ex分量测量,而海洋通讯系统应选择接收Hy信号;3)同线方向Ex分量相对于赤道方向对海底高阻层的反映明显.本文的研究内容对海洋可控源电磁法测量方式的选择以及通信技术的研究具有指导意义.     

有限长导线源频率测深有限内存拟牛顿一维反演

本文采用有限内存拟牛顿法实现有限长导线源频率测深阻抗响应数据的一维反演。水平层状介质有限长导线源阻抗频率响应由基于虚界面法获得的地表水平正交电场和磁场计算得到;一维反演优化问题的求解利用有限内存拟牛顿法,结合光滑模型约束,直接对阻抗的频率响应数据进行反演。在反演过程中,正则化参数的调整采用目标函数自适应技术。反演模型剖分为多层,各层厚度自地表按比例增加。反演从均匀半空间开始,终止条件为目标函数相对变化小于10~(-4)。分别对理论模型和实际数据进行了反演模拟。为考察反演的稳定性,还对理论数据添加10%随机噪声后进行了反演。数值计算结果表明:有限内存拟牛顿方法可以用于有限长导线源频率测深阻抗频率响应的反演;该反演方法对初始模型的依赖性弱,从均匀半空间模型出发基本可以恢复到真实模型;反演初期收敛较快,后期收敛速度变慢,反演结束一般需要迭代40次左右。噪声数据反演结果表明,随机噪声对反演结果影响不大,说明有限内存拟牛顿法具有较好的抗干扰能力。本文研究成果给出了可控源电磁数据反演的一种新方法;同时,利用本文的研究成果,可以为二维或三维反演建立合适的初始模型。