用于起伏地表条件下地球物理场数值模拟的正交曲网格生成技术 被引量：15

1

作者 孙建国 蒋丽丽 《石油地球物理勘探》 EI CSCD 北大核心 2009年第4期494-500,共7页

起伏地表条件下地球物理场数值模拟所面临的主要困难在于如何求解自由边界(地—空界面)上的边界条件。如今虽有多种解决途径,如映射法、插值法和不规则网格法,但均存在一些问题,本文针对上述方法的问题提出了一种新的解决途径。首先引... 起伏地表条件下地球物理场数值模拟所面临的主要困难在于如何求解自由边界(地—空界面)上的边界条件。如今虽有多种解决途径,如映射法、插值法和不规则网格法,但均存在一些问题,本文针对上述方法的问题提出了一种新的解决途径。首先引入了在计算流体力学中发展起来的贴体正交网格技术,通过对弹性波场和直流电场在地—空界面上所应满足的微分方程(边界条件)的分析,得到了正交曲网格优于非正交曲网格的结论;然后利用基于泊松方程的网格生成方法,采用Hilgenstock法控制源项,并结合代数辅助方法生成了正交贴体曲网格。结果表明,本文所用的网格生成方法适用于解决起伏地表条件下地球物理场的数值模拟问题。 展开更多

关键词 网格生成 数值模拟 泊松方程 正交贴体曲网格

在线阅读 下载PDF 职称材料

P-SV波在提取SS波零偏移距反射率中的应用 被引量：3

2

作者 孙鹏远 孙建国 卢秀丽 《石油物探》 EI CSCD 2004年第2期103-105,共3页

利用Gardner's经验关系式对P-SV波反射系数近似公式进行了简化,给出了一个描述P-SV波反射系数和SS波零偏移距反射率之间关系的统一近似公式。4个不同的含油气砂岩模型的定量计算表明,该统一近似公式不受岩性弱反差条件的限制,入射角... 利用Gardner's经验关系式对P-SV波反射系数近似公式进行了简化,给出了一个描述P-SV波反射系数和SS波零偏移距反射率之间关系的统一近似公式。4个不同的含油气砂岩模型的定量计算表明,该统一近似公式不受岩性弱反差条件的限制,入射角在30°以内能够准确地反映P-SV波反射系数和SS波零偏移距反射率之间的近似关系,而且能够适用于不同形式的P-SV波反射系数近似公式。因此,利用该关系式可以由P-SV波直接构建SS波零偏移距剖面,而且可以实现SS波阻抗的定量反演和多波AVO分析,还可以通过SS波资料提取P-SV波振幅随入射角变化的信息。 展开更多

关键词 P-SV波 SS波 零偏移距反射率 砂岩 入射角 反射系数 波阻抗 地震资料 地震勘探技术

在线阅读 下载PDF 职称材料

应用波前构建法实现零炮检距等时线的构建 被引量：2

3

作者 韩复兴 孙建国 王坤 《石油地球物理勘探》 EI CSCD 北大核心 2012年第2期207-211,224,共6页

本文针对基尔霍夫偏移及反偏移成像过程中所需的等时线构建问题,从波前构建的角度着手,分析了波前构建法在实现零炮检距等时线构建过程当中的理论基础和实现方案,并针对均匀介质模型,应用波前构建的思想构建不同时间步长上所对应的等时... 本文针对基尔霍夫偏移及反偏移成像过程中所需的等时线构建问题,从波前构建的角度着手,分析了波前构建法在实现零炮检距等时线构建过程当中的理论基础和实现方案,并针对均匀介质模型,应用波前构建的思想构建不同时间步长上所对应的等时线,同均匀介质模型理论计算的等时线进行误差对比,验证该方法的准确性。通过应用波前构建的思想实现梯度介质、渐变介质、层状介质以及Marmousi模型不同时间步长上等时线和射线路径的构建实例,验证该方法的可行性。 展开更多

关键词 波前构建 等时线 零炮检距 射线

在线阅读 下载PDF 职称材料

Kirchhoff型偏移理论的研究历史、研究现状与发展趋势展望——与光学绕射理论的类比、若干新结果、新认识以及若干有待于解决的问题 被引量：8

4

作者 孙建国 《吉林大学学报（地球科学版）》 EI CAS CSCD 北大核心 2012年第5期1521-1552,共32页

对Kirchhoff型偏移理论的研究历史、研究现状以及在今后一段时间内的发展趋势进行了简要的回顾和评述,并对加权函数、偏移孔径、子波拉伸、振幅畸变、垂向与横向分辨率等问题进行了相对详尽的讨论。还对若干有待于解决的问题进行了探讨... 对Kirchhoff型偏移理论的研究历史、研究现状以及在今后一段时间内的发展趋势进行了简要的回顾和评述,并对加权函数、偏移孔径、子波拉伸、振幅畸变、垂向与横向分辨率等问题进行了相对详尽的讨论。还对若干有待于解决的问题进行了探讨,并对笔者在近几年内所取得的一些与Kirchhoff型偏移理论有关的新认识和新结果进行了介绍。为了充分地顾及反射地震成像和光学成像之间的渊源关系,采用了一种新的讨论方式,即首先对Kirchhoff绕射理论的主要内容进行概略的回顾,然后再通过类比的方式引入Kirchhoff型偏移成像的基本概念和基本方程。在过去的几十年中,在Kirchhoff型偏移理论研究方面已经取得了重大进展,一个完整的理论体系已经形成。尽管如此,仍有一些遗留的或新出现的理论问题有待于解决,例如加权等时面分布叠加偏移的输出场分析和模型空间中的最优偏移孔径等等。由于Kirchhoff型偏移理论的研究对象是加权绕射叠加而不是加权等时面分布,也由于GPU/CPU协同计算技术的出现已经使得加权绕射叠加的计算量问题不再形成瓶颈,所以在今后的发展中Kirchhoff型偏移在理论和实践上将达到统一,即将在实现上从加权等时面分布回归到加权绕射叠加。与此同时,在今后一个时期内的研究重点将会与Green函数的数学表述、局部相关反射信号(locally coherent events)的表示和处理、广义加权绕射叠加以及在模型空间中定义的最优偏移孔径等问题有关。 展开更多

关键词 地震波 Kirchhoff型偏移理论 研究历史 研究现状 发展趋势

在线阅读 下载PDF 职称材料

基于近代天线理论的钻孔电磁波法正常场 被引量：2

5

作者 孙建国 《地球学报》 EI CAS CSCD 北大核心 2002年第6期567-571,共5页

钻孔电磁波法的正常场公式是进行各种钻孔电磁波探测理论研究和资料解释的基础,因此有必要将其建立在近代天线理论之上。通过对钻孔电磁波法中的基本电磁系统的分析讨论,证明了关于绝缘天线的广义传输线理论在钻孔电磁探测和电磁波测井... 钻孔电磁波法的正常场公式是进行各种钻孔电磁波探测理论研究和资料解释的基础,因此有必要将其建立在近代天线理论之上。通过对钻孔电磁波法中的基本电磁系统的分析讨论,证明了关于绝缘天线的广义传输线理论在钻孔电磁探测和电磁波测井中的可应用性,并给出了在一般条件下的基于广义传输线理论的钻孔电磁波法正常场公式。通过与经典的正常场公式的定性对比发现,基于广义传输线理论的正常场公式既包含了天线本身绝缘层的作用,又考虑了钻孔的影响。因此,本文提出的公式能提供更多的关于围岩和钻孔流体的电磁性质的信息,应能成为跨孔电磁成像和电磁波测井的现代理论基础。 展开更多

关键词 钻孔电磁波法 正常场 广义传输线理论 磁导率 围岩 钻孔流体

在线阅读 下载PDF 职称材料

Kirchhoff型高保真反偏移理论 被引量：1

6

作者 孙建国 《地球科学与环境学报》 CAS 2011年第2期207-212,216,共7页

为了使Kirchhoff型反偏移的输出结果与数值模拟结果相等,提出了一种旨在消除Kirchhoff型反偏移中所出现的振幅畸变现象的高保真反偏移方法。与常规的Kirchhoff型反偏移不同,高保真反偏移是一种与反射面有关的反射地震成像方法。因此,在... 为了使Kirchhoff型反偏移的输出结果与数值模拟结果相等,提出了一种旨在消除Kirchhoff型反偏移中所出现的振幅畸变现象的高保真反偏移方法。与常规的Kirchhoff型反偏移不同,高保真反偏移是一种与反射面有关的反射地震成像方法。因此,在高保真反偏移中首先要确定反射面的位置,其次要确定振幅畸变因子的数值,最后要根据在反射点上振幅畸变因子恒等于1这一事实进行振幅畸变校正。在具体实施过程中,为了求出反射点的空间位置和振幅畸变因子的具体数值,采用四重加权叠加法,即在一次反偏移运算中同时采用一个单位加权函数和三个非单位加权函数。其中,非单位加权函数由常规的真振幅加权函数和用于建立Kircohhoff型反偏移算子的全局坐标系中的水平坐标组成。通过在四个反偏移场之间进行比值运算,可以达到提取反射点坐标和振幅畸变因子的目的。一旦得到反射点坐标和振幅畸变因子的具体数值,就可以通过简单的除法消除掉振幅畸变因子对反偏移像场振幅的影响。 展开更多

关键词 Kirchhoff型反偏移 振幅畸变 等时面多重加权叠加 高保真

在线阅读 下载PDF 职称材料

论直流电位场拟解析近似中的电位与电荷反射函数的物理意义

7

作者 孙建国 《吉林大学学报（地球科学版）》 EI CAS CSCD 北大核心 2012年第2期545-553,共9页

为了避开由大型代数方程组数值求解所带来的问题,将电场反射张量的基本概念和拟解析近似的基本思想引入到了直流电位场的数值模拟中,定义了体电位反射函数、面电荷反射函数和面电位反射函数等3个辅助函数,并给出了其在一般条件下的拟解... 为了避开由大型代数方程组数值求解所带来的问题,将电场反射张量的基本概念和拟解析近似的基本思想引入到了直流电位场的数值模拟中,定义了体电位反射函数、面电荷反射函数和面电位反射函数等3个辅助函数,并给出了其在一般条件下的拟解析近似解和在简单条件下(均匀或点电流源背景场、均匀导电柱体或球体)的分析解。简单的对比分析证明:面电荷反射函数是异常体表面上自由面电荷密度的函数;在简单条件下,这种函数关系转化成为比例关系;类似的结论对于简单条件下的面电位反射函数也成立;然而,在一般条件下难以建立面电位反射函数和自由面电荷密度之间的明确关系。 展开更多

关键词 地球物理勘探-电法 直流电位场 拟解析近似 反射函数 物理意义

在线阅读 下载PDF 职称材料

复杂地表条件下基于线性插值和窄带技术的地震波走时计算 被引量：27

8

作者 孙章庆 孙建国 韩复兴 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2009年第11期2846-2853,共8页

为了在复杂地表条件下实现地震波走时计算,提出了一种基于线性插值和窄带技术的走时计算新方法.其中,线性插值用于局部走时计算,窄带技术用于局部波前捕获和追踪.为了逼近起伏地表,采用三角网和矩形网相结合的方法对速度模型进行剖分.... 为了在复杂地表条件下实现地震波走时计算,提出了一种基于线性插值和窄带技术的走时计算新方法.其中,线性插值用于局部走时计算,窄带技术用于局部波前捕获和追踪.为了逼近起伏地表,采用三角网和矩形网相结合的方法对速度模型进行剖分.为了得到局部走时计算公式,利用费马(Fermat)原理和关于入射点位置的限定条件.有关编程实践和数值试验表明:新方法不仅可以有效、灵活地处理地表高程的剧烈变化,而且还具有很好的适应性和稳定性,得到的计算结果满足波前传播规律. 展开更多

关键词 复杂地表 线性插值 窄带技术 混合网格 走时计算

在线阅读 下载PDF 职称材料

P-P波AVO近似对比研究:定性分析 被引量：16

9

作者 孙鹏远 孙建国 卢秀丽 《石油地球物理勘探》 EI CSCD 北大核心 2002年第B09期164-171,222,共8页

根据已发表的 AVO(Amplitude versus offset)文献资料 ,对基于 Zoeppritz方程的 13种 P- P波反射系数近似进行了系统归纳和总结 ,将其分成了 8大类 ,并讨论了各种近似的内在联系。同时 ,对各种近似的实质、产生背景、假设条件、特点及... 根据已发表的 AVO(Amplitude versus offset)文献资料 ,对基于 Zoeppritz方程的 13种 P- P波反射系数近似进行了系统归纳和总结 ,将其分成了 8大类 ,并讨论了各种近似的内在联系。同时 ,对各种近似的实质、产生背景、假设条件、特点及其在参数反演中的作用等进行了详尽的对比和评述。通过定性分析对比发现 :基于射线参数、幂级数及基于拉梅常数的近似公式可以描述具有大反差的弹性参数界面上的 P- P反射 ,特别适合于弹性参数的定量反演 ;AVO属性提取中所使用的背景速度信息的精度将影响反演的精度和可靠性。另外 ,根据弹性参数和波速之间的关系 ,可以在各种近似公式之间进行一定的转换。这些结论将对合理地利用 展开更多

关键词 P-P波反射系数近似 AVO公式 对比 ZOEPPRITZ方程 定性分析 地震勘探

在线阅读 下载PDF 职称材料

针对复杂地形的三种地震波走时算法及对比 被引量：19

10

作者 孙章庆 孙建国 韩复兴 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2012年第2期560-568,共9页

复杂地形条件下地震波走时算法对于研究复杂地形地区的成像问题有着重要的意义.为了得到精度高且适应于复杂地形的走时算法,首先提出阶梯网格迎风差分法.然后将该方法与不等距网格有限差分法和混合网格线性插值法进行对比研究,得出如下... 复杂地形条件下地震波走时算法对于研究复杂地形地区的成像问题有着重要的意义.为了得到精度高且适应于复杂地形的走时算法,首先提出阶梯网格迎风差分法.然后将该方法与不等距网格有限差分法和混合网格线性插值法进行对比研究,得出如下结论:混合网格线性插值法的计算精度最高,但其计算效率最低;阶梯网格迎风差分法的计算精度最低,但其计算效率最高;不等距网格有限差分法的计算精度和计算效率均居中;而究竟选取哪种算法作为给定复杂地形模型的地震波走时算法,应该综合考虑地形的特点、所研究问题对计算精度及计算效率的要求等因素.最后通过一个计算实例验证了三种算法在面对复杂地形、近地表及地下复杂介质等复杂地质条件时均有很好的适应性和稳定性. 展开更多

关键词 复杂地形 走时计算 迎风差分 不等距差分 线性插值法

在线阅读 下载PDF 职称材料

三维起伏地表条件下地震波走时计算的不等距迎风差分法 被引量：25

11

作者 孙章庆 孙建国 韩复兴 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2012年第7期2441-2449,共9页

三维起伏地表条件下的地震波走时计算技术是研究三维起伏地表地区很多地震数据处理技术的基础性工具.为了获得适应于任意三维起伏地表且计算精度高的走时算法,提出三维不等距迎风差分法.该方法采用不等距网格剖分三维起伏地表模型,通过... 三维起伏地表条件下的地震波走时计算技术是研究三维起伏地表地区很多地震数据处理技术的基础性工具.为了获得适应于任意三维起伏地表且计算精度高的走时算法,提出三维不等距迎风差分法.该方法采用不等距网格剖分三维起伏地表模型,通过在迎风差分格式中引入不等距差分格式、Huygens原理及Fermat原理来建立地表附近的局部走时计算公式,并通过在窄带技术中设定新的网格节点类型来获得三维起伏地表条件下算法的整体实现步骤.精度及算例分析表明:三维不等距迎风差分法具有很高的计算精度且能够适应于任意三维起伏地表模型. 展开更多

关键词 三维起伏地表 走时计算 不等距迎风差分法 窄带技术

在线阅读 下载PDF 职称材料

P-P波AVO近似对比研究:定量分析 被引量：12

12

作者 孙鹏远 孙建国 卢秀丽 《石油地球物理勘探》 EI CSCD 北大核心 2002年第B09期172-179,222,共8页

利用两种典型的含气砂岩模型和“P- P波 AVO近似对比研究 :定性分析”一文中所讨论过的各种 AVO近似公式 ,分别计算了在正、负波阻抗界面条件下的 P- P波反射系数随入射角的变化、反射系数的梯度和截距 ,以及与基于 Zoeppritz方程的解... 利用两种典型的含气砂岩模型和“P- P波 AVO近似对比研究 :定性分析”一文中所讨论过的各种 AVO近似公式 ,分别计算了在正、负波阻抗界面条件下的 P- P波反射系数随入射角的变化、反射系数的梯度和截距 ,以及与基于 Zoeppritz方程的解析解相对比的相对误差。结果发现 ,Xu、Mallick、Zheng xd、Wang、Goodway、Fatti等近似在本文选择的两种不同模型下均能与精确解很好地吻合 ,反映出很高的近似程度。另外 ,各种方法对不同模型的近似程度并不相同。只有将多种方法相结合 。 展开更多

关键词 P-P波 反射系数近似 AVO公式 对比 ZOEPPRITZ方程 定量分析 地震勘探

在线阅读 下载PDF 职称材料

P-SV波AVO分析 被引量：12

13

作者 孙鹏远 孙建国 卢秀丽 《石油地球物理勘探》 EI CSCD 北大核心 2003年第2期131-135,共5页

参照 Aki和 Richards的近似方法 ,本文对 Zoeppritz方程所描述的 P- SV波反射系数公式的角度项进行了泰勒展开 ,并保留了密度和速度的二阶项 ,得出一个新的反射系数近似公式。对 4个三层含气砂岩模型的定量计算表明 ,该近似公式所得出... 参照 Aki和 Richards的近似方法 ,本文对 Zoeppritz方程所描述的 P- SV波反射系数公式的角度项进行了泰勒展开 ,并保留了密度和速度的二阶项 ,得出一个新的反射系数近似公式。对 4个三层含气砂岩模型的定量计算表明 ,该近似公式所得出的计算结果在入射角小于 4 0°时与由 Zoeppritz方程给出的精确解吻合得很好 ,甚至优于 Aki和 Richards公式的近似效果 ,很适合大入射角的 AVO参数反演和理论研究。以新的近似公式为基础 ,文中分别从不同岩性参数表示的转换波反射系数、P- SV转换波的 AVO加权叠加方法及 P- P波和 P- 展开更多

关键词 P-SV转换波 反射系数近似 AVO分析 加权叠加 地震勘探 油气勘探

在线阅读 下载PDF 职称材料

基于复程函方程和改进的快速推进法的复旅行时计算方法 被引量：10

14

作者 黄兴国 孙建国 +1 位作者 孙章庆 王乾龙 《石油地球物理勘探》 EI CSCD 北大核心 2016年第6期1109-1118,1147,共11页

直接对复程函方程进行数值求解计算地震波复旅行时存在困难,因为复程函方程融合了复旅行时实、虚两部分。本文结合扰动理论和梯度点积而形成的一种最优化方法,使数值求解复程函方程来计算复旅行时成为可能。为了将这种计算方法应用于弯... 直接对复程函方程进行数值求解计算地震波复旅行时存在困难,因为复程函方程融合了复旅行时实、虚两部分。本文结合扰动理论和梯度点积而形成的一种最优化方法,使数值求解复程函方程来计算复旅行时成为可能。为了将这种计算方法应用于弯曲射线,引入了不等距迎风差分法,并提出了改进的多点起算快速推进法的方案,对以弯曲射线为边界的区域进行处理,从而准确求取复旅行时。数值结果表明,本文的计算方法和处理策略是有效的。 展开更多

关键词 复旅行时 复程函方程 扰动理论 不等距差分 快速推进法

在线阅读 下载PDF 职称材料

2维和2.5维起伏地表直流电法有限差分数值模拟 被引量：17

15

作者 张东良 孙建国 孙章庆 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2011年第1期234-244,共11页

起伏地表直流电场数值模拟现多采用有限元法,主要是因为其有灵活的处理曲边界的能力,然而有限元法比有限差分法要复杂,如果让有限差分法也同样具有较好的处理曲边界的能力,那数值模拟将变得更为简单.本文通过在非正则内点处采用不等距差... 起伏地表直流电场数值模拟现多采用有限元法,主要是因为其有灵活的处理曲边界的能力,然而有限元法比有限差分法要复杂,如果让有限差分法也同样具有较好的处理曲边界的能力,那数值模拟将变得更为简单.本文通过在非正则内点处采用不等距差分,在起伏地表点处直接实现边界条件,克服转移法人为改变地表形状的弊端,使得基于笛卡尔网格的有限差分法有了较强的处理曲边界的能力.通过2维和2.5维的水平地表数值解和解析解的对比,以及2维起伏地表情况下数值解和保角变换理论解的对比,均说明该方法的有效性.最后例举了一些数值计算实例. 展开更多

关键词 起伏地表 曲边界 数值模拟 有限差分

在线阅读 下载PDF 职称材料

三维复杂山地多级次群推进迎风混合法多波型走时计算 被引量：7

16

作者 孙章庆 孙建国 +2 位作者 王雪秋 高正辉 江兆南 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2017年第5期1861-1873,共13页

三维复杂山地条件下的各种地震波型的走时计算技术,可以直接用于复杂山地区域地震波运动学特性的分析、地震数据采集观测系统的设计以及直接基于三维复杂地表的地震数据处理技术的研发.为了在三维复杂地表条件下准确、灵活且稳定地计算... 三维复杂山地条件下的各种地震波型的走时计算技术,可以直接用于复杂山地区域地震波运动学特性的分析、地震数据采集观测系统的设计以及直接基于三维复杂地表的地震数据处理技术的研发.为了在三维复杂地表条件下准确、灵活且稳定地计算各种地震波型的走时,提出一种多级次群推进迎风混合法.该算法利用不等距迎风差分法简洁稳定地处理三维复杂地表及附近的局部走时计算问题,利用计算精度不错的迎风双线性插值法处理绝大部分均匀正方体网格中的局部走时计算问题,利用群推进法模拟三维复杂地表条件下地震波前的扩展问题,利用多级次算法处理各种类型的地震波的走时计算问题.算法分析和计算实例表明:新方法具有很好的计算精度与效率,且能灵活稳定地处理三维复杂地表复杂介质条件下的多波型走时计算问题. 展开更多

关键词 三维复杂山地 多级次算法 群推进法 迎风混合法 多波型 走时计算

在线阅读 下载PDF 职称材料

基于完全三叉树堆排序的波前扩展有限差分地震波走时快速算法 被引量：5

17

作者 杨昊 孙建国 +1 位作者 韩复兴 马淑芳 《吉林大学学报（地球科学版）》 EI CAS CSCD 北大核心 2010年第1期188-194,共7页

波前扩展有限差分地震波走时算法具有物理意义明确、因果稳定性强的特点,但每次波前扩展都要寻找波前面上的最小走时点。当计算网格点数较多,特别是涉及到三维走时计算时,寻找波前面上的最小走时点是一项十分耗时的工作。研究发现,波前... 波前扩展有限差分地震波走时算法具有物理意义明确、因果稳定性强的特点,但每次波前扩展都要寻找波前面上的最小走时点。当计算网格点数较多,特别是涉及到三维走时计算时,寻找波前面上的最小走时点是一项十分耗时的工作。研究发现,波前扩展有限差分地震波走时算法的波前点具有两个突出特点:①波前点更新十分频繁,通常每次取出波前最小走时点后都要插入若干新的波前点;②新计算出的波前点的走时通常比较大。数据结构中的二叉树堆排序方法可以提高寻找波前面上最小走时点的效率,根据特点①,在原始二叉树堆排序方法的基础上,优化了插入新波前点和移除波前面上最小走时点的流程,实际计算结果表明,与原始的二叉树堆排序方法相比,改进后的二叉树堆排序方法可以提高大约20%的计算效率。根据特点②,将原始的二叉树堆排序方法推广到多叉树,实际计算结果表明,完全三叉树堆排序方法优于二叉树和四叉树堆排序方法,可以再提高5%的计算效率。 展开更多

关键词 走时 波前扩展 程函方程 堆排序 有限差分法

在线阅读 下载PDF 职称材料

2.5维起伏地表条件下坐标变换法直流电场数值模拟 被引量：7

18

作者 孙章庆 孙建国 张东良 《吉林大学学报（地球科学版）》 EI CAS CSCD 北大核心 2010年第2期425-431,共7页

研究起伏地表对视电阻率分布的影响是进行地形校正的基础。由于很难处理不规则边界,计算简洁且效率高的有限差分方法很少用于解决起伏地表问题。为了解决该问题,引入曲化平思想,实现了一种基于坐标变换法的起伏地表条件下的直流电场数... 研究起伏地表对视电阻率分布的影响是进行地形校正的基础。由于很难处理不规则边界,计算简洁且效率高的有限差分方法很少用于解决起伏地表问题。为了解决该问题,引入曲化平思想,实现了一种基于坐标变换法的起伏地表条件下的直流电场数值模拟方法。方法从传统的2.5维基本方程及边界条件出发,通过坐标变换将起伏地表问题转化为水平地表问题,并利用有限差分法求解水平地表问题,最后再通过坐标映射得到起伏地表条件下的电位和视电阻率分布。精度分析及计算实例表明:本方法实现简洁,计算误差主要集中在震源附近,整个计算区域内的平均相对误差为1.39%,计算结果满足起伏地表条件下的电位和视电阻率的分布规律。 展开更多

关键词 2.5维 起伏地表 直流电场 坐标变换 有限差分法

在线阅读 下载PDF 职称材料

二维起伏地表条件下坐标变换法直流电场数值模拟 被引量：8

19

作者 孙章庆 孙建国 张东良 《吉林大学学报（地球科学版）》 EI CAS CSCD 北大核心 2009年第3期528-534,共7页

针对有限差分方法处理不规则边界难的问题,引入曲化平的思想,实现了一种新的基于坐标变换法的起伏地表条件下直流电场数值模拟方法。通过坐标变换,将物理空间的直流电场满足的基本方程和边界条件变换到计算空间;在计算空间内,利用有限... 针对有限差分方法处理不规则边界难的问题,引入曲化平的思想,实现了一种新的基于坐标变换法的起伏地表条件下直流电场数值模拟方法。通过坐标变换,将物理空间的直流电场满足的基本方程和边界条件变换到计算空间;在计算空间内,利用有限差分方法对变换后的方程和边界条件进行数值分析;最后,将计算结果由计算空间映射回物理空间。精度分析及计算实例表明:本文提出的算法精度高,计算误差主要集中在源附近,且当源离地表较近时平均相对误差在1.0%以下;计算结果满足起伏地表条件下电位及视电阻率的分布规律。 展开更多

关键词 起伏地表 直流电场 坐标变换 有限差分

在线阅读 下载PDF 职称材料

基于线性插值和窄带技术的走时计算方法 被引量：10

20

作者 孙章庆 孙建国 韩复兴 《石油地球物理勘探》 EI CSCD 北大核心 2009年第4期436-441,共6页

为了解决现有走时计算方法中计算精度和计算效率很难兼顾的问题,本文基于线性插值和窄带技术提出了一种新的走时计算方法。该方法在实现过程中利用线性插值技术进行局部走时计算,利用窄带技术进行波前传播过程的模拟。数值分析表明其计... 为了解决现有走时计算方法中计算精度和计算效率很难兼顾的问题,本文基于线性插值和窄带技术提出了一种新的走时计算方法。该方法在实现过程中利用线性插值技术进行局部走时计算,利用窄带技术进行波前传播过程的模拟。数值分析表明其计算精度和效率均高于基于窄带技术的一阶、二阶精度有限差分法;通过对几个典型的复杂速度模型的数值运算证实,该方法切实可行,且对复杂速度模型具有很好的适应性。 展开更多

关键词 线性插值 窄带技术 走时计算

在线阅读 下载PDF 职称材料