1、总目录总目录第一章 地球物理中的计算问题概述（3）第二章 地球物理中常用数值解法 的基本原理（14）第三章 电磁场数值模拟（5）第四章 弹性波场数值模拟（5)第五章 地球物理反问题基本方法（3）讨 论（2）上 机 自主选题、编程、上机计算、分析问题及解决方法（作业）第一章第一章 地球物理中的计算问题概述地球物理中的计算问题概述 内容提要内容提要第一节 绪论第二节 地球物理中的引力位方程第三节 地球物理中的磁位方程第四节 地球物理中的电位方程第五节 地球物理中的电磁场方程第六节 地球物理中的弹性波场方程第七节 地球物理中的声波方程 计算地球物理学（Computational Geophysics
2、)是地球物理、数学和计算机科学三者相结合而产生的一门边缘学科。计算地球物理学是地球物理学的一个分支学科，已经成为现代地球物理学理论与应用的各个方面必不可少的科学手段和有力支柱。第一节第一节 绪论绪论1.1 1.1 计算地球物理学的兴起与发展计算地球物理学的兴起与发展 地球物理学实质上是一门以观测数据对地球内部各种尺度的地质体和矿体进行研究的一门科学，因此，从它形成独立学科的时候起，就离不开计算问题。反问题 正问题 理论问题研究、提出新的理论或方法时第一节第一节 绪论绪论1.1 1.1 计算地球物理学的兴起与发展计算地球物理学的兴起与发展n 计算机问世前，地球物理的计算是粗略的。在上世纪50年代
3、以前虽然已经广泛地进行地球物理计算，但它并未形成一门学科分支。n 将计算机引进到地球物理工作中来，特别是数字化采样后，才使地球物理的计算工作向更深入、更精确的方向发展，许多以前无法计算的问题得以解决。这样才逐渐形成了计算地球物理学的分支学科。n 形成的开始时期大约在1960年代初期。第一节第一节 绪论绪论1.2 1.2 计算地球物理学的研究内容计算地球物理学的研究内容n正问题 物理模拟、解析解、数值计算 解析法分离变数法、积分变换方法、Green函数方法、变分法、对于二维和三维Laplace方程的边值问题，也还可以将解表示为特殊的积分公式、保角变换，对于双曲型方程的定解问题，也存在一些特殊的解
4、法，例如平均值法，降维法。第一节第一节 绪论绪论1.2 1.2 计算地球物理学的研究内容计算地球物理学的研究内容n 正问题 数值计算效率高、费用低。数值计算方法多种多样：有限差分法、有限元法、边界元法、快速傅里叶变换法、拟谱法、数值积分法、元胞自动机法和射线追踪法。第一节第一节 绪论绪论1.2 1.2 计算地球物理学的研究内容计算地球物理学的研究内容n反问题 有关地球内部的知识的来源地球物理观测。地球物理学的基本问题用地面或地表附近的各种观测资料定性和定量地对地下的地质结构和矿产资源做出判断，这就是解地球物理反问题。第一节第一节 绪论绪论1.2 1.2 计算地球物理学的研究内容计算地球物理学的
5、研究内容n 反问题 地球物理反问题：（1）根据现在的物理状态去确定物理过程的过去状态。（温度分布、地震偏移成像、电磁波成像和地震层析成像）。（2）从微分方程的解的某种泛函来求方程的系数或右端项（速度、密度、电阻率等参数的反演）。不适定问题、非线性。第一节第一节 绪论绪论1.2 1.2 计算地球物理学的研究内容计算地球物理学的研究内容n 信号处理、信息提取与可视化 反问题求解之前，要对实际地球物理数据进行旨在加强有用信息的信号处理。第一节第一节 绪论绪论1.3 1.3 （1）不适定问题不可避免，特别是在解反问题中 吉洪诺夫假定不适定问题的解具有某种“正则”性质，可以通过将不适定问题正则化为适定的
6、问题求解求得稳定近似解。适定性解的存在性、稳定性、唯一性第一节第一节 绪论绪论1.3 1.3 （1）不适定问题不可避免，特别是在解反问题中多解性第一节第一节 绪论绪论由反射系数求波阻抗1.3 1.3 （2）要求计算地球物理的求解精度要有适度性，否则会适得其反。地球物理的观测数据是间接的物理场响应的综合结果，不是与地下的地层空间点一一对应的关系，不能简单地要求达到想要达到的理想精度，例如地震分辨率、去噪、组合。第一节第一节 绪论绪论1.3 1.3 （3）在计算地球物理学中结合实用的创新性是它的发展的基础 计算地球物理学不是直接应用数理方程和计算数学中的已有知识，而是与地球物理问题密切结合，以明确
7、的地球物理问题为出发点和归宿点。从地球物理问题提出的方程、公式和求解方法与计算方法都是比较独特的，不但与现成的数学方法不完全一致，而且有时很难对它们的适定性、收敛性和稳定性做出预先的研究。常常是先应用，从应用的效果来判断模型、方法和计算过程的可行性。第一节第一节 绪论绪论 第二节第二节 地球物理中的引力位方程地球物理中的引力位方程 0000,Mx y zm,P x y z2mFGrFr0002PMPMPMmGrrrF0000,Mx y zm,P x y zFr第二节第二节 地球物理中的引力位方程地球物理中的引力位方程 QL第二节第二节 地球物理中的引力位方程地球物理中的引力位方程 0000,M
8、x y zm,P x y zFrQL由引力和功的可加性，引力场强和引力位的概念可推广到任意分布的质量体激发的引力场：第二节第二节 地球物理中的引力位方程地球物理中的引力位方程 从场位确定场强 第二节第二节 地球物理中的引力位方程地球物理中的引力位方程 重力场是保守场，是有位势的场，是无旋场（梯度的旋度=0）。引力位U满足泊松方程（非齐次）或拉普拉斯方程（齐次），这取决于体密度的分布情况 第二节第二节 地球物理中的引力位方程地球物理中的引力位方程 22222224,UUUUGx y zxyz 磁场起源于电荷的运动，因此，磁现象实质上是电流的一种效应。目前在磁法勘探中所研究的静磁场就是稳定电流的磁
9、场，它只与电流的分布有关。第三节第三节 地球物理中的磁位方程地球物理中的磁位方程 稳定电流的磁场与静电场不一样，它的磁力线是闭合的，因此它是个无散、有旋场。22222224,mUUUUx y zxyz 对于磁荷体密度,mx y z不为零的区域，静磁位满足泊松方程 第四节第四节 地球物理中的电位方程地球物理中的电位方程 SdIJS第四节第四节 地球物理中的电位方程地球物理中的电位方程 ,EEJJE第四节第四节 地球物理中的电位方程地球物理中的电位方程 0ldElMUdEl第四节第四节 地球物理中的电位方程地球物理中的电位方程 U E1UU J第四节第四节 地球物理中的电位方程地球物理中的电位方程
10、 0EU E20U,ex y z22222224,eUUUUx y zxyz 第四节第四节 地球物理中的电位方程地球物理中的电位方程 22222224,4,4,meGx y zUUUUx y zxyzx y z第五节第五节 地球物理中的电磁场方程地球物理中的电磁场方程 电磁勘探从理论基础并结合异常的原因分为：(1)传导类电法，观测和利用传导作用而产生的异常，如电阻率法、自然电场法、充电法和激发极化法等，它们满足位场理论。这种方法通常称为电法。（充电法是良导体直接连通电源，观测其电场分布特征和规律进行找矿勘探的一种方法。带电后的良导电体是一个近似等电位的带电体，可在地表、钻孔或坑道中观测电位或电
11、位梯度的变化。根据实测曲线可分析推断矿体形状，产状、埋深、几个矿体是否相连及确定地下水流速、流向等问题。）第五节第五节 地球物理中的电磁场方程地球物理中的电磁场方程 电磁勘探从理论基础并结合异常的原因分为：（2）感应类电法，观测和利用的是由电磁感应作用产生的异常，如瞬变电磁法、大地电磁测深法、频率测深法、无线电波法等，它们满足麦克斯韦方程，且又是一种波场，与地震波场在某些方法有相通之处。这种方法称为电磁法。瞬变电磁法利用不接地回线或接地线源向地下发射一次脉冲磁场，在一次脉冲磁场间歇期间利用线圈或接地电极观测地下介质中引起的二次感应涡流场，从而探测介质电阻率的一种方法。频率测深法指频率在几十到几
12、万Hz的音频范围内，通过改变交变磁场频率的办法探测岩层电阻率随深度的变化以了解地质构造和找矿的一种人工场源电磁法。无线电波频率在几十万赫至几十兆赫的电磁波。第五节第五节 地球物理中的电磁场方程地球物理中的电磁场方程 =+=0ttq DH JBEDBDEBHJE第五节第五节 地球物理中的电磁场方程地球物理中的电磁场方程 对于高频的探地雷达，电磁场满足麦克斯韦方程。对于大地电磁测深，一般f10 Hz，这时在导电介质中的位移电流与传导电流 J 相比可以忽略不计（），在谐变场的情况下，假定介质无磁性，则0=00iEHHEHE第五节第五节 地球物理中的电磁场方程地球物理中的电磁场方程 zy0,y z,y
13、 z二维介质麦克斯韦方程分解为两组独立的方程yzxxyxzEEiHyzHEzHEy yzxxyxzHHEyzEiHzEiHy 第五节第五节 地球物理中的电磁场方程地球物理中的电磁场方程 亥姆霍兹方程式 磁场平行极化 电场平行极化 222222222200 xxxxxxHHEEk Hk Eyzyz2ki第六节第六节 地球物理中的弹性波场方程地球物理中的弹性波场方程 CeTxxyyzzyzzxxy6 6ijcCTxxyyzzyzzxxyeeeeeee广义虎克定律 第六节第六节 地球物理中的弹性波场方程地球物理中的弹性波场方程 应变-位移关系,111,222yxzxxyyzzyyxxzzyzzxxy
14、uuueeexyzuuuuuueeezyzxxy第六节第六节 地球物理中的弹性波场方程地球物理中的弹性波场方程 应力与位移之间满足牛顿运动定律 第六节第六节 地球物理中的弹性波场方程地球物理中的弹性波场方程 各向同性介质 第六节第六节 地球物理中的弹性波场方程地球物理中的弹性波场方程 牛顿运动定律 2jiijijjiuuxx虎克定律 xyxxxxzxyyxyyyzyzyzxzzzzvftxyzvftxyzvftxyzvx，vy，vz为质点振动速度。第六节第六节 地球物理中的弹性波场方程地球物理中的弹性波场方程 牛顿运动定律 2ijijije虎克定律 xyxxxxzxyyxyyyzyzyzxzz
15、zzvftxyzvftxyzvftxyzijyjxizijjivvvvvtxyzxx第六节第六节 地球物理中的弹性波场方程地球物理中的弹性波场方程 应力和质点振动速度的一阶偏微分方程组 xyxxxxzxyyxyyyzyzyzxzzzzvftxyzvftxyzvftxyz,xyyyzyxxzxzzvvvvvvtyxtzxtzy222yxxxxzyyyyxzyxzzzzvvvvtxyzxvvvvtxyzyvvvvtxyzz第六节第六节 地球物理中的弹性波场方程地球物理中的弹性波场方程 虎克定律、应变-位移关系、牛顿运动定律二阶波动方程 第七节第七节 地球物理中的声波方程地球物理中的声波方程 固体介
16、质中的纵波是一种胀缩应变波，与流体中的声波具有同样性质。如果不考虑固体中的转换波问题，地震波的传播问题可以使用声波方程来研究。通常使用的标量波动方程就是这类以标量函数，如压力为未知函数的方程式。111yxzxyzvPvvKtxyzvPtxvPtyvPtz 第七节第七节 地球物理中的声波方程地球物理中的声波方程 二维声波方程和电磁波方程的相似性（一阶方程）11xzxzvvPKtxzvPtxvPtz 111xzxzHEHtzxHEtzHEtx 1,xzzxEPHvHvK第七节第七节 地球物理中的声波方程地球物理中的声波方程 二维声波方程和电磁波方程的相似性（二阶方程）111xzxzHEHtzxHEtzHEtx 2222222222111xzxzHEHtz tx tHEz tzHEx tx 2222221EEEtzx第七节第七节 地球物理中的声波方程地球物理中的声波方程 二维声波方程和电磁波方程的相似性（二阶方程）11xzxzvvPKtxzvPtxvPtz 2222222222211xzxzvpvvtx tz tvpx txvpz tz 2222222pppvtxz第七节第七节 地球物理中的声波方程地球物理中的声波方程 二维声波方程和电磁波方程的相似性2222222pppvtxz2211,cv2222222EEEctzx
