第1章地球物理中的计算问题概述

1、计算地球物理地球物理与信息工程学院 物探系周 辉2014年总目录第一章 地球物理中的计算问题概述（3）第二章 地球物理中常用数值解法 的基本原理（14）第三章 电磁场数值模拟（5）第四章 弹性波场数值模拟 （5)第五章 地球物理反问题基本方法（3）讨 论 （2）上 机 自主选题、编程、上机计算、 分析问题及解决方法 （作业）计算地球物理地球物理与信息工程学院 物探系周 辉2014年第一章 地球物理中的计算问题概述 内容第一节 绪论第二节 地球物理中的引力位方程第三节 地球物理中的磁位方程第四节 地球物理中的电位方程第五节 地球物理中的电磁场方程第六节 地球物理中的弹性波场方程第七节 地球物理中

2、的声波方程 计算地球物理学（Computational Geophysics)是地球物理、数学和计算机科学三者相结合而产生的一门边缘学科。 计算地球物理学是地球物理学的一个分支学科，已经成为现代地球物理学理论与应用的各个方面必不可少的科学手段和有力支柱。 第一节 绪论1.1 计算地球物理学的兴起与发展 地球物理学实质上是一门以观测数据对地球内部各种尺度的地质体和矿体进行研究的一门科学，因此，从它形成独立学科的时候起，就离不开计算问题。 反问题 正问题 理论问题研究、提出新的理论或方法时第一节 绪论1.1 计算地球物理学的兴起与发展n 计算机问世前，地球物理的计算是粗略的。在上世纪50年代以前虽

3、然已经广泛地进行地球物理计算，但它并未形成一门学科分支。n 将计算机引进到地球物理工作中来，特别是数字化采样后，才使地球物理的计算工作向更深入、更精确的方向发展，许多以前无法计算的问题得以解决。这样才逐渐形成了计算地球物理学的分支学科。n 形成的开始时期大约在1960年代初期。 第一节 绪论1.2 计算地球物理学的研究内容n正问题 物理模拟、解析解、数值计算 解析法分离变数法、积分变换方法、Green函数方法、变分法、对于二维和三维Laplace方程的边值问题，也还可以将解表示为特殊的积分公式、保角变换，对于双曲型方程的定解问题，也存在一些特殊的解法，例如平均值法，降维法。 第一节 绪论1.2

4、 计算地球物理学的研究内容n 正问题 数值计算效率高、费用低。 数值计算方法多种多样：有限差分法、有限元法、边界元法、快速傅里叶变换法、拟谱法、数值积分法（克希霍夫积分）、元胞自动机法（格子玻尔兹曼法）和射线追踪法。 第一节 绪论1.2 计算地球物理学的研究内容n反问题 有关地球内部的知识的来源地球物理观测。 地球物理学的基本问题用地面或地表附近的各种观测资料定性和定量地对地下的地质结构和矿产资源做出判断，这就是解地球物理反问题。 第一节 绪论1.2 计算地球物理学的研究内容n 反问题 地球物理反问题： （1）根据现在的物理状态去确定物理过程的过去状态。（温度分布、地震偏移成像、电磁波成像等）

5、。 （2）从微分方程的解的某种泛函来求方程的系数或右端项（速度、密度、电阻率等参数的反演）。不适定问题、非线性。 第一节 绪论1.2 计算地球物理学的研究内容n 信号处理、信息提取与可视化 反问题求解之前，要对实际地球物理数据进行旨在加强有用信息的信号处理。 第一节 绪论1.3 计算地球物理学的特点 （1）不适定问题不可避免，特别是在解反问题中 吉洪诺夫假定不适定问题的解具有某种“正则”性质，可以通过将不适定问题正则化为适定的问题求解求得稳定近似解。 适定性解的存在性、稳定性、唯一性第一节 绪论1.3 计算地球物理学的特点 （1）不适定问题不可避免，特别是在解反问题中多解性第一节 绪论例子1由

6、反射系数求波阻抗例子2 重力异常1.3 计算地球物理学的特点 （2）计算地球物理的求解精度要求有适度性，否则会适得其反。 地球物理的观测数据是间接的物理场响应的综合结果，不是与地下的地层空间点一一对应的关系，不能简单地要求达到想要达到的理想精度，例如地震分辨率、去噪、组合。 第一节 绪论h4 hoo C CD菲涅尓带1.3 计算地球物理学的特点 （3）在计算地球物理学中结合实用的创新性是它的发展的基础 计算地球物理学不是直接应用数理方程和计算数学中的已有知识，而是与地球物理问题密切结合，以明确的地球物理问题为出发点和归宿点。 从地球物理问题提出的方程、公式和求解方法与计算方法都是比较独特的，不

7、但与现成的数学方法不完全一致，而且有时很难对它们的适定性、收敛性和稳定性做出预先的研究。常常是先应用，从应用的效果来判断模型、方法和计算过程的可行性。第一节 绪论 第二节 地球物理中的引力位方程 0000,Mx y zm, ,P x y z2mFGr矢量形式 Fr0002PMPMPMmGrrrF对点P 的引力F大小为 分量形式 对x，y，z的偏导数 0000,Mx y zm, ,P x y zFr引力场对单位质量点所作的功 第二节 地球物理中的引力位方程 QL功与运动所经过的路径无关而只决定于运动的起始点和终点的位置。这个性质叫力场的保守性。对于保守力场可引入一个标量函数来描述。固定一点Q，对

8、空间任意点P（M0点除外）赋予一个标量V(P)单位质点从Q点移动到点P引力所作的功： 第二节 地球物理中的引力位方程 Q处于无穷远处0000,Mx y zm, ,P x y zFrQLV(P)为引力场的位。引力场既是保守场又是有位势的场 。由引力和功的可加性，引力场强和引力位的概念可推广到任意分布的质量体激发的引力场：第二节 地球物理中的引力位方程 从场位确定场强 第二节 地球物理中的引力位方程 重力场是保守场，是有位势的场，是无旋场（梯度的旋度=0）。引力位 U 满足泊松方程（非齐次）或拉普拉斯方程（齐次），这取决于体密度的分布情况 第二节 地球物理中的引力位方程 22222224, ,UU

9、UUGx y zxyz 磁场起源于电荷的运动，因此，磁现象实质上是电流的一种效应。目前在磁法勘探中所研究的静磁场就是稳定电流的磁场，它只与电流的分布有关。 第三节 地球物理中的磁位方程 稳定电流的磁场与静电场不一样，它的磁力线是闭合的，因此它是个无散、有旋场。22222224, ,mUUUUx y zxyz 对于磁荷体密度, ,mx y z不为零的区域，静磁位满足泊松方程 传导类电法勘探，是以地壳中人工电流场、天然电流场在岩层或矿体中的空间分布规律和由于矿体被极化而产生的激发极化电流场的时间特征及空间分布规律为理论基础的。 第四节 地球物理中的电位方程 地中电流分布的连续性 通过不包含电流源的

10、任意闭合面的电流密度通量等于零。SdIJS第四节 地球物理中的电位方程 地中电流密度与电场强度的正比性 ,EEJJE欧姆定律的微分形式 在解决电法勘探正、反演问题时最为常用。 第四节 地球物理中的电位方程 地中电流场的势场性电场作功与所沿路径无关,即电场为无旋场。 在稳定电流场中，地中任意点M处的电位，等于单位正电荷由观察点M处移到无穷远电场力所作的功，0ldElMUdEl第四节 地球物理中的电位方程 地中电流场的势场性稳定电流场中电场强度与电位的关系 U E1UU J第四节 地球物理中的电位方程 地中电流场的势场性0EU E20U, ,ex y z22222224, ,eUUUUx y zx

11、yz 在电荷密度 不为零的区域，电位满足第四节 地球物理中的电位方程 重力位、磁位和电位满足的方程重写如下 22222224, ,4, ,4, ,meGx y zUUUUx y zxyzx y z第五节 地球物理中的电磁场方程 电磁勘探从理论基础并结合异常的原因分为：(1)传导类电法，观测和利用传导作用而产生的异常，如电阻率法、自然电场法、充电法和激发极化法等，它们满足位场理论。这种方法通常称为电法。 （充电法是良导体直接连通电源，观测其电场分布特征和规律进行找矿勘探的一种方法。带电后的良导电体是一个近似等电位的带电体，可在地表、钻孔或坑道中观测电位或电位梯度的变化。根据实测曲线可分析推断矿体

12、形状，产状、埋深、几个矿体是否相连及确定地下水流速、流向等问题。）第五节 地球物理中的电磁场方程 电磁勘探从理论基础并结合异常的原因分为：（2）感应类电法，观测和利用的是由电磁感应作用产生的异常，如瞬变电磁法、大地电磁测深法、频率测深法、无线电波法等，它们满足麦克斯韦方程，且又是一种波场，与地震波场在某些方法有相通之处。这种方法称为电磁法。 瞬变电磁法利用不接地回线或接地线源向地下发射一次脉冲磁场，在一次脉冲磁场间歇期间利用线圈或接地电极观测地下介质中引起的二次感应涡流场，从而探测介质电阻率的一种方法。频率测深法指频率在几十到几万Hz的音频范围内，通过改变交变磁场频率的办法探测岩层电阻率随深度

13、的变化以了解地质构造和找矿的一种人工场源电磁法。无线电波频率在几十万赫至几十兆赫的电磁波。 第五节 地球物理中的电磁场方程 麦克斯韦方程组= +=0ttq DH JBEDB麦克斯韦方程组是大地电磁场所遵循的基本微分方程式。E 电场强度H 磁场强度D 电感应强度B 磁感应强度J 电流密度q 自由电荷 DEBHJE第五节 地球物理中的电磁场方程 对于高频的探地雷达，电磁场满足麦克斯韦方程。对于大地电磁测深，一般 f10 Hz，这时在导电介质中的位移电流与传导电流 J 相比可以忽略不计（ ），在谐变场的情况下，假定介质无磁性，则0=00iEHHEHE第五节 地球物理中的电磁场方程 zy0, y z,

14、 y z二维介质麦克斯韦方程分解为两组独立的方程yzxxyxzEEiHyzHEzHEy 磁场平行极化 电场平行极化 yzxxyxzHHEyzEiHzEiHy 第五节 地球物理中的电磁场方程 亥姆霍兹方程式 磁场平行极化 电场平行极化 222222222200 xxxxxxHHEEk Hk Eyzyz2ki第六节 地球物理中的弹性波场方程 CeTxxyyzzyzzxxy6 6ijcCTxxyyzzyzzxxyeeeeeee广义虎克定律 第六节 地球物理中的弹性波场方程 应变-位移关系 ,111,222yxzxxyyzzyyxxzzyzzxxyuuueeexyzuuuuuueeezyzxxy第六节

15、 地球物理中的弹性波场方程 应力与位移之间满足牛顿运动定律 第六节 地球物理中的弹性波场方程 各向同性介质 第六节 地球物理中的弹性波场方程 牛顿运动定律 2jiijijjiuuxx虎克定律 xyxxxxzxyyxyyyzyzyzxzzzzvftxyzvftxyzvftxyzvx，vy，vz为质点振动速度。第六节 地球物理中的弹性波场方程 牛顿运动定律 2ijijije虎克定律 xyxxxxzxyyxyyyzyzyzxzzzzvftxyzvftxyzvftxyzijyjxizijjivvvvvtxyzxx第六节 地球物理中的弹性波场方程 应力和质点振动速度的一阶偏微分方程组 xyxxxxzxy

16、yxyyyzyzyzxzzzzvftxyzvftxyzvftxyz,xyyyzyxxzxzzvvvvvvtyxtzxtzy222yxxxxzyyyyxzyxzzzzvvvvtxyzxvvvvtxyzyvvvvtxyzz第六节 地球物理中的弹性波场方程 虎克定律、应变-位移关系、牛顿运动定律二阶波动方程 第七节 地球物理中的声波方程 固体介质中的纵波是一种胀缩应变波，与流体中的声波具有同样性质。如果不考虑固体中的转换波问题，地震波的传播问题可以使用声波方程来研究。通常使用的标量波动方程就是这类以标量函数，如压力为未知函数的方程式。 111yxzxyzvPvvKtxyzvPtxvPtyvPtz 第七节 地球物理中的声波方程 二维声波方程和电磁波方程的相似性（一阶方程）11xzxzvvPKtxzvPtxvPtz 111xzxzHEHtzxHEtzHEtx 1,xzzxEPHvHvK第七节 地球物理中的声波方程 二维声波方程和电磁波方程的相似性（二阶