各向异性地层中可控源电磁法一维全参数反演及三维有限体积正演算法研究

来源 :吉林大学 | 被引量 : 0次 | 上传用户:yaping3211
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
可控源音频大地电磁法(CSAMT)是在大地电磁法(MT)和音频大地电磁法(AMT)的基础上发展起来的一种地面人工源频率域测深方法,针对MT法场源的随机性和信号微弱,以致观测十分困难这一状况,CSAMT采用可以控制的人工场源,CSAMT广泛的应用于地质普查、勘探石油、天然气、地热、金属矿产、水文、工程、环境保护等各种复杂地电构造的勘探,是一种重要的地球物理勘探方法。一维水平层状各向同性地层中的CSAMT正反演是非常重要的,因为一维模型简单,存在解析解,作为一个基本的地电模型,有助于理解电磁场的基本物理特性,它也是实际数据的初步解释和分析以及多维反演的初始模型选取的基础。实际地层在一定尺度上总是存在电的各向异性,例如,砂岩页岩的薄交互层是一种典型的各向异性构造。对各向异性地层测量数据进行解释,如果忽略各向异性,将很难得到有效的反演结果。因此,研究各向异性地层中的CSAMT的反演理论是一个非常重要的课题。本文研究各向异性地层中CSAMT的一维正反演理论,讨论各向异性地层对CSAMT一维正反演的影响特征和规律,为实际资料的正确解释提供理论依据和指导。与其他地球物理问题的反演一样,各向异性地层中的一维CSAMT反演,涉及正演计算,模型空间的选择,Fréchet导数的计算,以及反演过程中模型的迭代修正来得到有效的反演结果。在第二章,利用二维Fourier变换与矩阵分解技术将层状各向异性地层中Maxwell方程转化成两个独立的关于TM波和TE波的传输线方程;借助传输线理论与叠加原理,根据CSAMT只使用电偶极源的特点,改进传输线算法,只引入电流源传输线Green函数,来求解TM波和TE波,建立各向异性地层中频率-波数域电流源电场和磁场并矢Green函数的新算法和新的表达式,相比传统的传输线理论,提高了CSAMT数值模拟的效率。在此基础上,利用传输线Green函数的基本解以及边界条件,推导出广义反射系数与振幅递推公式,得到各个地层中传输线Green函数的解析解;然后利用Fourier逆变换与Bessel公式将CSAMT响应表示为Sommerfeld形式的积分,针对层状模型的特点,本文提出了三次样条插值结合Lommel积分公式的技术,快速计算其数值解,相比地球物理电磁法中常用的数字滤波算法,精度相当,但计算速度更快。在第三章,选择能够定量反演地下模型参数的参数化模型空间,通过水平和垂直导电率以及水平层界面构成的模型矢量来描述电阻率的空间分布。利用并矢Green函数,结合摄动理论,推导得到了Fréchet导数的解析表达式,表示为Sommerfeld形式的积分,采用第二章提出的三次样条插值结合Lommel积分公式的技术可以有效的求解。进一步的,建立了全参数广义逆反演算法和全参数正则化迭代反演算法来同时反演所有的模型参数,包括横向、纵向电导率和层界面深度。广义逆反演算法采用奇异值分解技术结合阻尼最小二乘法,可以有效的提高反演结果的稳定性和抗噪能力。对于正则化反演,由于CSAMT响应是非线性的,引入偏差原理并结合Cholesky分解自适应选择正则化因子,保证反演过程的稳定性。数值计算结果表明,CSAMT响应对地层纵向电导率的灵敏度要远小于相比对其他模型参数(横向电导率和层界面)的灵敏度。尽管初始模型存在较大误差,反演结果依然能够给出地层模型的主要特征。选择合适的初始层界面可以有效的改善反演结果;全参数的反演相比固定层界面的反演,反演结果更好;即使初始模型的层数与实际地层不符,依然能够给出有效的反演结果。海洋可控源电磁法(MCSEM)是一种用于海底油气勘探的频率域电磁方法。MCSEM具有提供海底地层电阻率空间分布的能力,在油气层识别和海上油气储层定量评价以及降低海上钻探风险方面发挥着重要作用,已发展成为海上油气勘探的一种重要方法。由于海底地形构造复杂以及地层横向电阻率分布不均匀,在海洋电磁勘探的设计以及海洋电磁资料处理和解释过程中,均需要进行大量的数值模拟,一维和二维的数值模拟技术已经较为成熟或者正在趋于成熟,三维数值模拟成为当前MCSEM的一个研究热点,目前模拟三维MCSEM的数值方法主要有有限元法、有限差分法、有限体积法和积分方程法。基于Yee氏交错网格的有限体积法,对Maxwell方程在各个单元上的积分进行离散处理,能有效降低方程的微分阶数,同时也减少了地层电导率不连续对离散结果的影响,在电磁场数值模拟中得到了较广泛应用。本文采用有限体积算法模型各向异性地层中MCSEM的三维响应,分析、总结各向异性对海洋CSEM三维电磁响应的影响特征。在第四章,为了有效模拟各向异性地层中海洋可控源电磁法的三维响应,我们建立了一套基于电场矢势与标势分解的耦合势有限体积法。利用电场的矢势和标势分解,将电场分解为无散场和无旋场之和,Maxwell方程转换为关于矢势与标势的混合Helmholtz方程,有效的克服低感应数问题。在此基础上,借助Yee氏交错网格和有限体积法推导出旋度和散度的离散公式以及非均质单元中等效电导率的计算公式,建立混合Helmholtz方程的离散方程。为保证大范围电磁场分布的稳定精确的求解以及多发射源的快速正演,采用直接法求解器PARDISO求解离散方程。此外,为了在不明显降低计算效率的情况下尽可能提高近场的计算精度,采用差异场方法处理各向异性地层3D模拟过程中的源奇异性问题。对高阻油气藏的各向异性的数值模拟表明,MCSEM沿测线方向的电场,随着油气藏的纵向电阻率的增大而增大,而与油气藏的横向电阻率大小无关。对油气藏上方的覆盖层的各向异性的数值模拟表明,沿测线方向的电场,随着覆盖层横向电阻率的增大而增大,同时随着覆盖层纵向电阻率的增大而增大。因此,在进行海洋CSEM三维数据解释时,要特别注意海底各向异性的影响。
其他文献
由于中西方的历史文化背景不同,作为西方道德学鼻祖的苏格拉底与中国古代思想家孔子的哲学思想呈现出独自的特点,本文试从人本主义、至善理论、教育思想和死亡哲学等方面来剖
<正>近年来,四川省成都市国家税务局从减法入手,减企业成本、减无效供给;以加法破局,加补齐短板、加服务创新;正面倡导诚信纳税,反面打击违法犯罪……这一加一减、一正一反,
目的通过分析南京地区2013-2017年团体献血者的年龄、性别、职业、血型及献血合格情况,为开展团体无偿献血动员工作及制订团体献血计划提供理论依据。方法采用整群抽样的方法
文章提出了一种基于STM32的智能多路无线温度采集系统的设计。该智能温度采集系统分为从机与主机两部分,主机和从机都是以STM32F103RCT6单片机作为中央处理器,并通过NRF24L01
新世纪教育学界的关注对象和关注范围相比过去,出现了很多变化。在变革中引入新的教学理论和教学方法,来提升整体的教学效率,是很多从业者都在认真思考的问题。在诸多的理论
在质疑中坚持自我,在怀疑中不断成长。正是果敢不羁、好奇心和探索精神,成就了这位诺奖得主。“我想要告诉年轻人的重要的一点是,当你有新颖的想法时,要对自己有信心。”2013
<正> 生产硝酸的工厂排放的尾气中,含有一定量残余的氧化氮(简称NOx),共含量主要取决于工厂类型,以及操作方法。通常尾气中NOx的含量为0.2—0.5%,有时还要高些。这样,硝酸生产
<正>2011年10月,上海大学钱伟长学院成为教育部首批"国家试点学院"之一。围绕国家和上海市发展战略,近年来,学院开展了一系列教育教学改革试验,实施"重基础、跨学科、国际化"
<正>目的在无偿献血队伍中,年满十八周岁的高校大学生成为了无偿献血队伍中的生力军,而该群体多为首次接触献血,因此如何高效选择献血者尤为重要。本文旨在分析健康征询在高
会议
<正>学校组织变革像其他任何组织变革一样,都是事关组织发展的重要举措,几乎在任何一所学校都随时随地发生着,但却常常被大部分管理者所忽视,有的缺少系统思考,没有顶层设计;