第1章地球物理中的计算问题概述论文资料.ppt

计算地球物理计算地球物理地球物理与信息工程学院地球物理与信息工程学院 物探系物探系周周 辉辉2012年年总目录总目录第一章第一章 地球物理中的计算问题概述地球物理中的计算问题概述 （（3 3））第二章第二章 地球物理中常用数值解法地球物理中常用数值解法 的基本原理的基本原理（（1414））第三章第三章 电磁场数值模拟电磁场数值模拟（（5 5））第四章第四章 弹性波场数值模拟弹性波场数值模拟 （（5)5)第五章第五章 地球物理反问题基本方法地球物理反问题基本方法（（3 3））讨讨 论论 （（2 2））上上 机机 自主选题、编程、上机计算、自主选题、编程、上机计算、 分析问题及解决方法分析问题及解决方法 （作业）（作业）计算地球物理计算地球物理地球物理与信息工程学院地球物理与信息工程学院 物探系物探系周周 辉辉2011年年第一章第一章 地球物理中的计算问题概述地球物理中的计算问题概述 内容提要内容提要第一节第一节 绪论绪论第二节第二节 地球物理中的引力位方程地球物理中的引力位方程第三节第三节 地球物理中的磁位方程地球物理中的磁位方程第四节第四节 地球物理中的电位方程地球物理中的电位方程第五节第五节 地球物理中的电磁场方程地球物理中的电磁场方程第六节第六节 地球物理中的弹性波场方程地球物理中的弹性波场方程第七节第七节 地球物理中的声波方程地球物理中的声波方程  计算地球物理学（计算地球物理学（Computational Geophysics)是是地球物理、数学和计算机科学三者相结合而产生的一地球物理、数学和计算机科学三者相结合而产生的一门门边缘学科边缘学科。

计算地球物理学是地球物理学的一个计算地球物理学是地球物理学的一个分支学科分支学科，，已经成为现代地球物理学理论与应用的各个方面必不已经成为现代地球物理学理论与应用的各个方面必不可少的科学手段和有力支柱可少的科学手段和有力支柱 第一节第一节 绪论绪论1.1 1.1 计算地球物理学的兴起与发展计算地球物理学的兴起与发展 地球物理学地球物理学实质上是一门以观测数据对地球内部实质上是一门以观测数据对地球内部各种尺度的地质体和矿体进行研究的一门科学，因此，各种尺度的地质体和矿体进行研究的一门科学，因此，从它形成独立学科的时候起，就离不开计算问题从它形成独立学科的时候起，就离不开计算问题 反问题反问题 正问题正问题 理论问题研究、提出新的理论或方法时理论问题研究、提出新的理论或方法时第一节第一节 绪论绪论1.1 1.1 计算地球物理学的兴起与发展计算地球物理学的兴起与发展n计算机问世前，地球物理的计算是粗略的在上世计算机问世前，地球物理的计算是粗略的在上世纪纪5050年代以前虽然已经广泛地进行地球物理计算，年代以前虽然已经广泛地进行地球物理计算，但它并未形成一门学科分支。

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

些特殊的解法，例如平均值法，降维法些特殊的解法，例如平均值法，降维法些特殊的解法，例如平均值法，降维法 第一节第一节 绪论绪论1.2 1.2 计算地球物理学的研究内容计算地球物理学的研究内容n n正问题正问题正问题正问题 数值计算数值计算数值计算数值计算————效率高、费用低效率高、费用低效率高、费用低效率高、费用低 数值计算方法多种多样：有限差分法、有限元法、数值计算方法多种多样：有限差分法、有限元法、数值计算方法多种多样：有限差分法、有限元法、数值计算方法多种多样：有限差分法、有限元法、边界元法、快速傅里叶变换法、拟谱法、数值积分边界元法、快速傅里叶变换法、拟谱法、数值积分边界元法、快速傅里叶变换法、拟谱法、数值积分边界元法、快速傅里叶变换法、拟谱法、数值积分法、元胞自动机法和射线追踪法法、元胞自动机法和射线追踪法法、元胞自动机法和射线追踪法法、元胞自动机法和射线追踪法 第一节第一节 绪论绪论1.2 1.2 计算地球物理学的研究内容计算地球物理学的研究内容n n反问题反问题反问题反问题 有关地球内部的知识的来源有关地球内部的知识的来源有关地球内部的知识的来源有关地球内部的知识的来源————地球物理观测。

地球物理观测地球物理观测地球物理观测 地球物理学的基本问题用地面或地表附近的各种观地球物理学的基本问题用地面或地表附近的各种观地球物理学的基本问题用地面或地表附近的各种观地球物理学的基本问题用地面或地表附近的各种观测资料定性和定量地对地下的地质结构和矿产资源测资料定性和定量地对地下的地质结构和矿产资源测资料定性和定量地对地下的地质结构和矿产资源测资料定性和定量地对地下的地质结构和矿产资源做出判断，这就是解地球物理反问题做出判断，这就是解地球物理反问题做出判断，这就是解地球物理反问题做出判断，这就是解地球物理反问题 第一节第一节 绪论绪论1.2 1.2 计算地球物理学的研究内容计算地球物理学的研究内容n n反问题反问题反问题反问题 地球物理反问题：地球物理反问题：地球物理反问题：地球物理反问题： （（（（1 1 1 1））））根据现在的物理状态去确定物理过程的过去根据现在的物理状态去确定物理过程的过去根据现在的物理状态去确定物理过程的过去根据现在的物理状态去确定物理过程的过去状态。

状态温度分布、地震偏移成像、电磁波成像和（温度分布、地震偏移成像、电磁波成像和（温度分布、地震偏移成像、电磁波成像和（温度分布、地震偏移成像、电磁波成像和地震层析成像）地震层析成像）地震层析成像）地震层析成像） （（（（2 2 2 2））））从微分方程的解的某种泛函来求方程的系数从微分方程的解的某种泛函来求方程的系数从微分方程的解的某种泛函来求方程的系数从微分方程的解的某种泛函来求方程的系数或右端项或右端项或右端项或右端项（速度、密度、电阻率等参数的反演）速度、密度、电阻率等参数的反演）速度、密度、电阻率等参数的反演）速度、密度、电阻率等参数的反演）不适定问题、非线性不适定问题、非线性不适定问题、非线性不适定问题、非线性 第一节第一节 绪论绪论1.2 1.2 计算地球物理学的研究内容计算地球物理学的研究内容n n信号处理、信息提取与可视化信号处理、信息提取与可视化信号处理、信息提取与可视化信号处理、信息提取与可视化 反问题求解之前，要对实际地球物理数据进行旨在反问题求解之前，要对实际地球物理数据进行旨在反问题求解之前，要对实际地球物理数据进行旨在反问题求解之前，要对实际地球物理数据进行旨在加强有用信息的信号处理。

加强有用信息的信号处理加强有用信息的信号处理加强有用信息的信号处理 第一节第一节 绪论绪论1.3 1.3 计算地球物理学的特点计算地球物理学的特点计算地球物理学的特点计算地球物理学的特点 （（（（1 1 1 1）不适定问题不可避免，特别是在解反问题中）不适定问题不可避免，特别是在解反问题中）不适定问题不可避免，特别是在解反问题中）不适定问题不可避免，特别是在解反问题中 吉洪诺夫假定不适定问题的解具有某种吉洪诺夫假定不适定问题的解具有某种吉洪诺夫假定不适定问题的解具有某种吉洪诺夫假定不适定问题的解具有某种“ “正则正则正则正则” ”性性性性质，可以通过将不适定问题正则化为适定的问题求质，可以通过将不适定问题正则化为适定的问题求质，可以通过将不适定问题正则化为适定的问题求质，可以通过将不适定问题正则化为适定的问题求解解解解————求得稳定近似解求得稳定近似解求得稳定近似解求得稳定近似解 适定性适定性适定性适定性————解的存在性、稳定性、唯一性解的存在性、稳定性、唯一性解的存在性、稳定性、唯一性解的存在性、稳定性、唯一性第一节第一节 绪论绪论1.3 1.3 计算地球物理学的特点计算地球物理学的特点计算地球物理学的特点计算地球物理学的特点 （（（（1 1 1 1）不适定问题不可避免，特别是在解反问题中）不适定问题不可避免，特别是在解反问题中）不适定问题不可避免，特别是在解反问题中）不适定问题不可避免，特别是在解反问题中多解性多解性多解性多解性第一节第一节 绪论绪论由由由由反反反反射射射射系系系系数数数数求求求求波波波波阻阻阻阻抗抗抗抗1.3 1.3 计算地球物理学的特点计算地球物理学的特点计算地球物理学的特点计算地球物理学的特点 （（（（2 2 2 2）要求计算地球物理的求解精度要有适度性，）要求计算地球物理的求解精度要有适度性，）要求计算地球物理的求解精度要有适度性，）要求计算地球物理的求解精度要有适度性，否则会适得其反。

否则会适得其反否则会适得其反否则会适得其反 地球物理的观测数据是间接的物理场响应的综合地球物理的观测数据是间接的物理场响应的综合地球物理的观测数据是间接的物理场响应的综合地球物理的观测数据是间接的物理场响应的综合结果，不是与地下的地层空间点一一对应的关系，结果，不是与地下的地层空间点一一对应的关系，结果，不是与地下的地层空间点一一对应的关系，结果，不是与地下的地层空间点一一对应的关系，不能简单地要求达到想要达到的理想精度，例如不能简单地要求达到想要达到的理想精度，例如不能简单地要求达到想要达到的理想精度，例如不能简单地要求达到想要达到的理想精度，例如地震分辨率、去噪、组合地震分辨率、去噪、组合地震分辨率、去噪、组合地震分辨率、去噪、组合 第一节第一节 绪论绪论1.3 1.3 计算地球物理学的特点计算地球物理学的特点计算地球物理学的特点计算地球物理学的特点 （（（（3 3 3 3）在计算地球物理学中结合实用的创新性是它的发）在计算地球物理学中结合实用的创新性是它的发）在计算地球物理学中结合实用的创新性是它的发）在计算地球物理学中结合实用的创新性是它的发展的基础展的基础展的基础展的基础 计算地球物理学不是直接应用数理方程和计算数计算地球物理学不是直接应用数理方程和计算数计算地球物理学不是直接应用数理方程和计算数计算地球物理学不是直接应用数理方程和计算数学中的已有知识，而是与地球物理问题密切结合，以学中的已有知识，而是与地球物理问题密切结合，以学中的已有知识，而是与地球物理问题密切结合，以学中的已有知识，而是与地球物理问题密切结合，以明确的地球物理问题为出发点和归宿点。

明确的地球物理问题为出发点和归宿点明确的地球物理问题为出发点和归宿点明确的地球物理问题为出发点和归宿点 从地球物理问题提出的方程、公式和求解方法与从地球物理问题提出的方程、公式和求解方法与从地球物理问题提出的方程、公式和求解方法与从地球物理问题提出的方程、公式和求解方法与计算方法都是比较独特的，不但与现成的数学方法不计算方法都是比较独特的，不但与现成的数学方法不计算方法都是比较独特的，不但与现成的数学方法不计算方法都是比较独特的，不但与现成的数学方法不完全一致，而且有时很难对它们的适定性、收敛性和完全一致，而且有时很难对它们的适定性、收敛性和完全一致，而且有时很难对它们的适定性、收敛性和完全一致，而且有时很难对它们的适定性、收敛性和稳定性做出预先的研究常常是先应用，从应用的效稳定性做出预先的研究常常是先应用，从应用的效稳定性做出预先的研究常常是先应用，从应用的效稳定性做出预先的研究常常是先应用，从应用的效果来判断模型、方法和计算过程的可行性果来判断模型、方法和计算过程的可行性果来判断模型、方法和计算过程的可行性果来判断模型、方法和计算过程的可行性。

第一节第一节 绪论绪论#### 第二节第二节 地球物理中的引力位方程地球物理中的引力位方程 m矢量形式矢量形式矢量形式矢量形式 Fr对点对点对点对点P P P P 的引力的引力的引力的引力F F大小为大小为大小为大小为 分量形式分量形式分量形式分量形式 对对对对x x，，，，y y，，，，z z的偏导数的偏导数的偏导数的偏导数 mFr引力场对单位质量点所作的功引力场对单位质量点所作的功引力场对单位质量点所作的功引力场对单位质量点所作的功 第二节第二节 地球物理中的引力位方程地球物理中的引力位方程 QL功与运动所经过的路径无关而只决定于运动的起功与运动所经过的路径无关而只决定于运动的起功与运动所经过的路径无关而只决定于运动的起功与运动所经过的路径无关而只决定于运动的起始点和终点的位置这个性质叫力场的始点和终点的位置这个性质叫力场的始点和终点的位置这个性质叫力场的始点和终点的位置这个性质叫力场的保守性保守性保守性保守性对于保守力场可引入一个标量函数来描述固定一点对于保守力场可引入一个标量函数来描述固定一点对于保守力场可引入一个标量函数来描述固定一点对于保守力场可引入一个标量函数来描述。

固定一点Q Q，对空，对空，对空，对空间任意点间任意点间任意点间任意点P P（（（（MM0 0点除外）赋予一个标量点除外）赋予一个标量点除外）赋予一个标量点除外）赋予一个标量V V( (P P)—)—单位质点从单位质点从单位质点从单位质点从Q Q点点点点移动到点移动到点移动到点移动到点P P引力所作的功：引力所作的功：引力所作的功：引力所作的功： 第二节第二节 地球物理中的引力位方程地球物理中的引力位方程 Q Q处于无穷远处处于无穷远处处于无穷远处处于无穷远处mFrQLV V( (P P) )为引力场的位为引力场的位为引力场的位为引力场的位引力场既是引力场既是引力场既是引力场既是保守场又是有位势的场保守场又是有位势的场保守场又是有位势的场保守场又是有位势的场 由引力和功的可加性，引力场强和引力位的概念可推广到任由引力和功的可加性，引力场强和引力位的概念可推广到任由引力和功的可加性，引力场强和引力位的概念可推广到任由引力和功的可加性，引力场强和引力位的概念可推广到任意分布的质量体激发的引力场：意分布的质量体激发的引力场：意分布的质量体激发的引力场：意分布的质量体激发的引力场：第二节第二节 地球物理中的引力位方程地球物理中的引力位方程 从场位确定场强从场位确定场强从场位确定场强从场位确定场强 第二节第二节 地球物理中的引力位方程地球物理中的引力位方程 重力场是保守场，是有位势的场，是无旋场（梯度的旋度重力场是保守场，是有位势的场，是无旋场（梯度的旋度重力场是保守场，是有位势的场，是无旋场（梯度的旋度重力场是保守场，是有位势的场，是无旋场（梯度的旋度=0=0）。

引力位引力位引力位引力位U U满足泊松方程（非齐次）或拉普拉斯方程（齐次）满足泊松方程（非齐次）或拉普拉斯方程（齐次）满足泊松方程（非齐次）或拉普拉斯方程（齐次）满足泊松方程（非齐次）或拉普拉斯方程（齐次），这取决于体密度的分布情况，这取决于体密度的分布情况，这取决于体密度的分布情况，这取决于体密度的分布情况 第二节第二节 地球物理中的引力位方程地球物理中的引力位方程 磁场起源于电荷的运动，因此，磁现象实质上是电流的一种磁场起源于电荷的运动，因此，磁现象实质上是电流的一种磁场起源于电荷的运动，因此，磁现象实质上是电流的一种磁场起源于电荷的运动，因此，磁现象实质上是电流的一种效应目前在磁法勘探中所研究的静磁场就是稳定电流的磁效应目前在磁法勘探中所研究的静磁场就是稳定电流的磁效应目前在磁法勘探中所研究的静磁场就是稳定电流的磁效应目前在磁法勘探中所研究的静磁场就是稳定电流的磁场，它只与电流的分布有关场，它只与电流的分布有关场，它只与电流的分布有关场，它只与电流的分布有关 第三节第三节 地球物理中的磁位方程地球物理中的磁位方程 稳定电流的磁场与静电场不一样，它的磁力线是闭合的，因稳定电流的磁场与静电场不一样，它的磁力线是闭合的，因稳定电流的磁场与静电场不一样，它的磁力线是闭合的，因稳定电流的磁场与静电场不一样，它的磁力线是闭合的，因此它是个无散、有旋场。

此它是个无散、有旋场此它是个无散、有旋场此它是个无散、有旋场对于磁荷体密度对于磁荷体密度对于磁荷体密度对于磁荷体密度不为零的区域，不为零的区域，不为零的区域，不为零的区域，静磁位满足泊松方程静磁位满足泊松方程静磁位满足泊松方程静磁位满足泊松方程 传导类电法勘探，是以地壳中人工电流场、天然电流场在岩传导类电法勘探，是以地壳中人工电流场、天然电流场在岩传导类电法勘探，是以地壳中人工电流场、天然电流场在岩传导类电法勘探，是以地壳中人工电流场、天然电流场在岩层或矿体中的空间分布规律和由于矿体被极化而产生的激发层或矿体中的空间分布规律和由于矿体被极化而产生的激发层或矿体中的空间分布规律和由于矿体被极化而产生的激发层或矿体中的空间分布规律和由于矿体被极化而产生的激发极化电流场的时间特征及空间分布规律为理论基础的极化电流场的时间特征及空间分布规律为理论基础的极化电流场的时间特征及空间分布规律为理论基础的极化电流场的时间特征及空间分布规律为理论基础的 第四节第四节 地球物理中的电位方程地球物理中的电位方程 地中电流分布的连续性地中电流分布的连续性地中电流分布的连续性地中电流分布的连续性 通过不包含电流源的任意闭合面的电流密度通量等于零。

通过不包含电流源的任意闭合面的电流密度通量等于零通过不包含电流源的任意闭合面的电流密度通量等于零通过不包含电流源的任意闭合面的电流密度通量等于零第四节第四节 地球物理中的电位方程地球物理中的电位方程 地中电流密度与电场强度的正比性地中电流密度与电场强度的正比性地中电流密度与电场强度的正比性地中电流密度与电场强度的正比性 ————欧姆定律的微分形式欧姆定律的微分形式欧姆定律的微分形式欧姆定律的微分形式 在解决电法勘探正、反演问题时最为常用在解决电法勘探正、反演问题时最为常用在解决电法勘探正、反演问题时最为常用在解决电法勘探正、反演问题时最为常用 第四节第四节 地球物理中的电位方程地球物理中的电位方程 地中电流场的势场性地中电流场的势场性地中电流场的势场性地中电流场的势场性电场作功与所沿路径无关电场作功与所沿路径无关电场作功与所沿路径无关电场作功与所沿路径无关, , , ,即电场为无旋场即电场为无旋场即电场为无旋场即电场为无旋场 在稳定电流场中，地中任意点在稳定电流场中，地中任意点在稳定电流场中，地中任意点在稳定电流场中，地中任意点M M M M处的电位，等于单位正电处的电位，等于单位正电处的电位，等于单位正电处的电位，等于单位正电荷由观察点荷由观察点荷由观察点荷由观察点M M M M处移到无穷远电场力所作的功，处移到无穷远电场力所作的功，处移到无穷远电场力所作的功，处移到无穷远电场力所作的功，第四节第四节 地球物理中的电位方程地球物理中的电位方程 地中电流场的势场性地中电流场的势场性地中电流场的势场性地中电流场的势场性稳定电流场中电场强度与电位的关系稳定电流场中电场强度与电位的关系稳定电流场中电场强度与电位的关系稳定电流场中电场强度与电位的关系 第四节第四节 地球物理中的电位方程地球物理中的电位方程 地中电流场的势场性地中电流场的势场性地中电流场的势场性地中电流场的势场性在电荷密度在电荷密度在电荷密度在电荷密度 不为零的区域，电位满足不为零的区域，电位满足不为零的区域，电位满足不为零的区域，电位满足第四节第四节 地球物理中的电位方程地球物理中的电位方程 重力位、磁位和电位满足的方程重写如下重力位、磁位和电位满足的方程重写如下重力位、磁位和电位满足的方程重写如下重力位、磁位和电位满足的方程重写如下 第五节第五节 地球物理中的电磁场方程地球物理中的电磁场方程 电磁勘探从理论基础并结合异常的原因分为：电磁勘探从理论基础并结合异常的原因分为：电磁勘探从理论基础并结合异常的原因分为：电磁勘探从理论基础并结合异常的原因分为：(1)(1)传导类电法传导类电法传导类电法传导类电法，观测和利用传导作用而产生的异常，，观测和利用传导作用而产生的异常，，观测和利用传导作用而产生的异常，，观测和利用传导作用而产生的异常，如电阻率法、自然电场法、充电法和激发极化法如电阻率法、自然电场法、充电法和激发极化法如电阻率法、自然电场法、充电法和激发极化法如电阻率法、自然电场法、充电法和激发极化法等，它们满足位场理论。

这种方法通常称为等，它们满足位场理论这种方法通常称为等，它们满足位场理论这种方法通常称为等，它们满足位场理论这种方法通常称为电法电法电法电法 （充电法是良导体直接连通电源，观测其电场分布特征和规律进行找矿（充电法是良导体直接连通电源，观测其电场分布特征和规律进行找矿（充电法是良导体直接连通电源，观测其电场分布特征和规律进行找矿（充电法是良导体直接连通电源，观测其电场分布特征和规律进行找矿勘探的一种方法带电后的良导电体是一个近似等电位的带电体，可勘探的一种方法带电后的良导电体是一个近似等电位的带电体，可勘探的一种方法带电后的良导电体是一个近似等电位的带电体，可勘探的一种方法带电后的良导电体是一个近似等电位的带电体，可在地表、钻孔或坑道中观测电位或电位梯度的变化根据实测曲线可在地表、钻孔或坑道中观测电位或电位梯度的变化根据实测曲线可在地表、钻孔或坑道中观测电位或电位梯度的变化根据实测曲线可在地表、钻孔或坑道中观测电位或电位梯度的变化根据实测曲线可分析推断矿体形状，产状、埋深、几个矿体是否相连及确定地下水流分析推断矿体形状，产状、埋深、几个矿体是否相连及确定地下水流分析推断矿体形状，产状、埋深、几个矿体是否相连及确定地下水流分析推断矿体形状，产状、埋深、几个矿体是否相连及确定地下水流速、流向等问题。

速、流向等问题速、流向等问题速、流向等问题第五节第五节 地球物理中的电磁场方程地球物理中的电磁场方程 电磁勘探从理论基础并结合异常的原因分为：电磁勘探从理论基础并结合异常的原因分为：电磁勘探从理论基础并结合异常的原因分为：电磁勘探从理论基础并结合异常的原因分为：（（（（2 2 2 2）感应类电法）感应类电法）感应类电法）感应类电法，观测和利用的是由电磁感应作用，观测和利用的是由电磁感应作用，观测和利用的是由电磁感应作用，观测和利用的是由电磁感应作用产生的异常，如瞬变电磁法、大地电磁测深法、产生的异常，如瞬变电磁法、大地电磁测深法、产生的异常，如瞬变电磁法、大地电磁测深法、产生的异常，如瞬变电磁法、大地电磁测深法、频率测深法、无线电波法等，频率测深法、无线电波法等，频率测深法、无线电波法等，频率测深法、无线电波法等，它们满足麦克斯韦它们满足麦克斯韦它们满足麦克斯韦它们满足麦克斯韦方程，且又是一种波场，与地震波场在某些方法方程，且又是一种波场，与地震波场在某些方法方程，且又是一种波场，与地震波场在某些方法方程，且又是一种波场，与地震波场在某些方法有相通之处有相通之处有相通之处有相通之处。

这种方法称为这种方法称为这种方法称为这种方法称为电磁法电磁法电磁法电磁法 瞬变电磁法瞬变电磁法瞬变电磁法瞬变电磁法——利用不接地回线或接地线源向地下发射一次脉冲磁场，在一利用不接地回线或接地线源向地下发射一次脉冲磁场，在一利用不接地回线或接地线源向地下发射一次脉冲磁场，在一利用不接地回线或接地线源向地下发射一次脉冲磁场，在一次脉冲磁场间歇期间利用线圈或接地电极观测地下介质中引起的二次感应次脉冲磁场间歇期间利用线圈或接地电极观测地下介质中引起的二次感应次脉冲磁场间歇期间利用线圈或接地电极观测地下介质中引起的二次感应次脉冲磁场间歇期间利用线圈或接地电极观测地下介质中引起的二次感应涡流场，从而探测介质电阻率的一种方法涡流场，从而探测介质电阻率的一种方法涡流场，从而探测介质电阻率的一种方法涡流场，从而探测介质电阻率的一种方法频率测深法频率测深法频率测深法频率测深法——指频率在几十到几万指频率在几十到几万指频率在几十到几万指频率在几十到几万HzHzHzHz的音频范围内，通过改变交变磁场频的音频范围内，通过改变交变磁场频的音频范围内，通过改变交变磁场频的音频范围内，通过改变交变磁场频率的办法探测岩层电阻率随深度的变化以了解地质构造和找矿的一种人工率的办法探测岩层电阻率随深度的变化以了解地质构造和找矿的一种人工率的办法探测岩层电阻率随深度的变化以了解地质构造和找矿的一种人工率的办法探测岩层电阻率随深度的变化以了解地质构造和找矿的一种人工场源电磁法。

场源电磁法场源电磁法场源电磁法无线电波无线电波无线电波无线电波————频率在几十万赫至几十兆赫的电磁波频率在几十万赫至几十兆赫的电磁波频率在几十万赫至几十兆赫的电磁波频率在几十万赫至几十兆赫的电磁波 第五节第五节 地球物理中的电磁场方程地球物理中的电磁场方程 麦克斯韦方程组麦克斯韦方程组麦克斯韦方程组麦克斯韦方程组麦克斯韦方程组是大地电磁场所遵循的基本微分方程式麦克斯韦方程组是大地电磁场所遵循的基本微分方程式麦克斯韦方程组是大地电磁场所遵循的基本微分方程式麦克斯韦方程组是大地电磁场所遵循的基本微分方程式E E 电场强度电场强度电场强度电场强度H H 磁场强度磁场强度磁场强度磁场强度D D 电感应强度电感应强度电感应强度电感应强度B B 磁感应强度磁感应强度磁感应强度磁感应强度J J 电流密度电流密度电流密度电流密度q q 自由电荷自由电荷自由电荷自由电荷 第五节第五节 地球物理中的电磁场方程地球物理中的电磁场方程 对于高频的对于高频的对于高频的对于高频的探地雷达探地雷达探地雷达探地雷达，电磁场满足麦克斯韦方程电磁场满足麦克斯韦方程电磁场满足麦克斯韦方程电磁场满足麦克斯韦方程。

对于对于对于对于大地电磁测深大地电磁测深大地电磁测深大地电磁测深，一般，一般，一般，一般f f＜＜＜＜1010 HzHz，这时在导电介质中的位，这时在导电介质中的位，这时在导电介质中的位，这时在导电介质中的位移电流与传导电流移电流与传导电流移电流与传导电流移电流与传导电流 J J 相比可以忽略不计（相比可以忽略不计（相比可以忽略不计（相比可以忽略不计（ ），在），在），在），在谐变场的情况下，假定介质无磁性，则谐变场的情况下，假定介质无磁性，则谐变场的情况下，假定介质无磁性，则谐变场的情况下，假定介质无磁性，则第五节第五节 地球物理中的电磁场方程地球物理中的电磁场方程 zy0二维介质二维介质二维介质二维介质麦克斯韦方程分解为两组独立的方程麦克斯韦方程分解为两组独立的方程麦克斯韦方程分解为两组独立的方程麦克斯韦方程分解为两组独立的方程磁场平行磁场平行磁场平行磁场平行极化极化极化极化 电场平行电场平行电场平行电场平行极化极化极化极化 第五节第五节 地球物理中的电磁场方程地球物理中的电磁场方程 亥姆霍兹方程式亥姆霍兹方程式亥姆霍兹方程式亥姆霍兹方程式 磁场平行极化磁场平行极化磁场平行极化磁场平行极化 电场平行极化电场平行极化电场平行极化电场平行极化 第六节第六节 地球物理中的弹性波场方程地球物理中的弹性波场方程 广义虎克定律广义虎克定律广义虎克定律广义虎克定律 第六节第六节 地球物理中的弹性波场方程地球物理中的弹性波场方程 应变应变应变应变- - - -位移关系位移关系位移关系位移关系 第六节第六节 地球物理中的弹性波场方程地球物理中的弹性波场方程 应力与位移之间满足牛顿运动定律应力与位移之间满足牛顿运动定律应力与位移之间满足牛顿运动定律应力与位移之间满足牛顿运动定律 第六节第六节 地球物理中的弹性波场方程地球物理中的弹性波场方程 各向同性介质各向同性介质各向同性介质各向同性介质 第六节第六节 地球物理中的弹性波场方程地球物理中的弹性波场方程 牛顿运动定律牛顿运动定律牛顿运动定律牛顿运动定律 虎克定律虎克定律虎克定律虎克定律 vx，vy，vz为质点振动速度。

第六节第六节 地球物理中的弹性波场方程地球物理中的弹性波场方程 牛顿运动定律牛顿运动定律牛顿运动定律牛顿运动定律 虎克定律虎克定律虎克定律虎克定律 第六节第六节 地球物理中的弹性波场方程地球物理中的弹性波场方程 应力和质点振动速度的一阶偏微分方程组应力和质点振动速度的一阶偏微分方程组应力和质点振动速度的一阶偏微分方程组应力和质点振动速度的一阶偏微分方程组 第六节第六节 地球物理中的弹性波场方程地球物理中的弹性波场方程 虎克定律、应变虎克定律、应变虎克定律、应变虎克定律、应变- - - -位移关系、牛顿运动定律位移关系、牛顿运动定律位移关系、牛顿运动定律位移关系、牛顿运动定律二阶波动方程二阶波动方程二阶波动方程二阶波动方程 第七节第七节 地球物理中的声波方程地球物理中的声波方程 固体介质中的纵波是一种胀缩应变波，与流体中的声波具有同固体介质中的纵波是一种胀缩应变波，与流体中的声波具有同固体介质中的纵波是一种胀缩应变波，与流体中的声波具有同固体介质中的纵波是一种胀缩应变波，与流体中的声波具有同样性质如果不考虑固体中的转换波问题，地震波的传播问题如果不考虑固体中的转换波问题，地震波的传播问题如果不考虑固体中的转换波问题，地震波的传播问题如果不考虑固体中的转换波问题，地震波的传播问题可以使用声波方程来研究。

可以使用声波方程来研究可以使用声波方程来研究可以使用声波方程来研究通常使用的标量波动方程就是这类通常使用的标量波动方程就是这类通常使用的标量波动方程就是这类通常使用的标量波动方程就是这类以标量函数，如压力为未知函数的方程式以标量函数，如压力为未知函数的方程式以标量函数，如压力为未知函数的方程式以标量函数，如压力为未知函数的方程式 第七节第七节 地球物理中的声波方程地球物理中的声波方程 二维声波方程和电磁波方程的相似性（一阶方程）二维声波方程和电磁波方程的相似性（一阶方程）二维声波方程和电磁波方程的相似性（一阶方程）二维声波方程和电磁波方程的相似性（一阶方程）第七节第七节 地球物理中的声波方程地球物理中的声波方程 二维声波方程和电磁波方程的相似性（二阶方程）二维声波方程和电磁波方程的相似性（二阶方程）二维声波方程和电磁波方程的相似性（二阶方程）二维声波方程和电磁波方程的相似性（二阶方程）第七节第七节 地球物理中的声波方程地球物理中的声波方程 二维声波方程和电磁波方程的相似性（二阶方程）二维声波方程和电磁波方程的相似性（二阶方程）二维声波方程和电磁波方程的相似性（二阶方程）二维声波方程和电磁波方程的相似性（二阶方程）第七节第七节 地球物理中的声波方程地球物理中的声波方程 二维声波方程和电磁波方程的相似性二维声波方程和电磁波方程的相似性二维声波方程和电磁波方程的相似性二维声波方程和电磁波方程的相似性。