有限差分法不同边界条件下的数值模拟
有限差分法的原理与计算步骤
有限差分法的原理与计算步骤有限差分法（Finite Difference Method）是一种常用的数值计算方法，用于求解偏微分方程的数值解。
其基本原理是将连续的偏微分方程转化为差分方程，通过逼近导数，使用离散的点代替连续的点，从而将问题转化为代数问题。
下面将详细介绍有限差分法的原理和计算步骤：一、基本原理：有限差分法基于Taylor级数展开，通过利用函数在其中一点附近的导数信息来逼近函数在该点处的值。
该方法将连续的偏微分方程转化为差分方程，使用离散的点代替连续的点，从而将问题转化为代数问题。
在有限差分法中，常用的差分逼近方式有前向差分、后向差分和中心差分。
二、计算步骤：1.网格划分：将求解区域划分为有限个离散点，并定义网格上的节点和网格尺寸。
通常使用等距离网格，即每个网格点之间的间距相等。
2.离散化：将偏微分方程中的各个导数项进行逼近，利用差分近似来替代和求解。
一般采用中心差分逼近方式，即通过函数值在两侧点的差来逼近导数。
3.代数方程系统：利用离散化的差分方程，将偏微分方程转化为代数方程系统。
根据问题的边界条件和初值条件，构建代数方程系统的系数矩阵和常数向量。
4. 求解代数方程：利用求解线性方程组的方法求解代数方程系统，常用的方法有直接法（如高斯消元法、LU分解法）和迭代法（如Jacobi迭代法、Gauss-Seidel迭代法）。
求解得到各个离散点的解。
5.后处理：根据求解结果进行后处理，包括结果的插值和可视化。
将离散点的解通过插值方法进行平滑处理，并进行可视化展示，以得到连续的函数解。
三、优缺点：1.直观：有限差分法基于网格划分，易于理解和实现。
2.精度可控：可通过调整网格大小和差分逼近方式来控制计算的精度。
3.广泛适用性：可用于求解各种偏微分方程，适用于不同的边界条件和初值条件。
然而，有限差分法也存在一些缺点：1.精度依赖网格：计算结果的精度受到网格划分的影响，因此需要谨慎选择网格大小。
2.限制条件：有限差分法适用于边界对应点处导数有定义的问题，不适用于奇异点和非线性问题。
浓度场数值模拟matlab
浓度场数值模拟matlab1. 介绍浓度场数值模拟是一种通过数值方法来计算和预测物质浓度分布的技术。
这种模拟方法通常基于流体力学和传质过程的基本方程，通过数值求解这些方程来得到浓度场的分布。
在很多领域中，如环境科学、化学工程、材料科学等，浓度场数值模拟都是一种重要的工具。
在本文中，我们将介绍如何使用MATLAB进行浓度场数值模拟。
我们将首先介绍数值模拟的基本原理和方法，然后详细讨论如何在MATLAB中实现这些方法，并给出一些实际应用的示例。
2. 数值模拟方法数值模拟浓度场的方法有很多种，其中常用的方法包括有限差分法（finite difference method）、有限元法（finite element method）和有限体积法（finite volume method）等。
这些方法都是将浓度场分割成网格，然后通过求解离散化的方程来得到各个网格点上的浓度值。
在MATLAB中，我们可以使用这些方法的内置函数或者自己编写代码来实现浓度场的数值模拟。
下面我们将以有限差分法为例，介绍如何在MATLAB中实现浓度场的数值模拟。
有限差分法是一种将微分方程离散化的方法，它将连续的空间域划分为离散的网格点，并在每个网格点上计算浓度的近似值。
有限差分法的基本思想是使用差分近似来代替微分运算，从而将微分方程转化为代数方程。
在浓度场数值模拟中，有限差分法通常用于求解扩散方程。
3. 在MATLAB中实现数值模拟在MATLAB中，我们可以使用内置函数pdepe来求解偏微分方程。
这个函数可以用于求解一维、二维和三维的偏微分方程，并支持不同的边界条件和初值条件。
下面是一个使用pdepe函数求解一维扩散方程的示例：function [c,f,s] = diffusion_eqn(x,t,u,DuDx)c = 1;f = DuDx;s = 0;function diffusion_simulation()x = linspace(0,1,100);t = linspace(0,1,100);m = 0;sol = pdepe(m,@diffusion_eqn,@diffusion_ic,@diffusion_bc,x,t);u = sol(:,:,1);surf(x,t,u)在这个示例中，我们定义了一个名为diffusion_eqn的函数来描述扩散方程。
有限差分法
第四章有限差分方法4.1引言有限差分法：数值求解常微分方程或偏微分方程的方法。
物理学和其他学科领域的许多问题在被分析研究之后, 往往可以归结为常微分方程或偏微分方程的求解问题。
一般说来，处理一个特定的物理问题，除了需要知道它满足的数学方程外，还应当同时知道这个问题的定解条件，然后才能设计出行之有效的计算方法来求解。
有限差分法以变量离散取值后对应的函数值来近似微分方程中独立变量的连续取值。
在有限差分方法中，我们放弃了微分方程中独立变量可以取连续值的特征，而关注独立变量离散取值后对应的函数值。
但是从原则上说，这种方法仍然可以达到任意满意的计算精度。
因为方程的连续数值解可以通过减小独立变量离散取值的间格，或者通过离散点上的函数值插值计算来近似得到。
这种方法是随着计算机的诞生和应用而发展起来的。
其计算格式和程序的设计都比较直观和简单，因而，它的实际应用已经构成了计算数学和计算物理的重要组成部分。
有限差分法的具体操作分为两个部分：(1)用差分代替微分方程中的微分，将连续变化的变量离散化，从而得到差分方程组的数学形式； (2)求解差分方程组。
在第一步中，我们通过所谓的网络分割法，将函数定义域分成大量相邻而不重合的子区域。
通常采用的是规则的分割方式。
这样可以便于计算机自动实现和减少计算的复杂性。
网络线划分的交点称为节点。
若与某个节点P 相邻的节点都是定义在场域内的节点，则P 点称为正则节点；反之，若节点P 有处在定义域外的相邻节点，则P 点称为非正则节点。
在第二步中，数值求解的关键就是要应用适当的计算方法，求得特定问题在所有这些节点上的离散近似值。
有限差分法的差分格式：一个函数在x 点上的一阶和二阶微商，可以近似地用它所临近的两点上的函数值的差分来表示。
如对一个单变量函数f(x)，x 为定义在区间[a,b]的连续变量。
以步长h=Δx 将[a,b]区间离散化，我们得到一系列节点x = a , x = x + h , x = x + h = a + 212132Δx , ..., x = x + h = b , 然后求出 f(x)在这些点上的近似值。
有限差分法在低渗油藏数值模拟中的应用
（２ １）
等 ．
（ ３ １）
＋ 三 ＋［
（ 一 ｓ ） ｓ
一
］ 户一 ・ ［＋－ ｉ ］ 。 一 ， ． ・ ’ ・，吉 Ｆ吉 ・ ． － ‘ 。件
＋ 警 一
ｓ ｓ ｏ—
一 警
（ ４ １ ）
初 始条 件
ｐ Ｏ＝ ｐ ０
维普资讯
２０ 年第 ９ ０８ 期
内蒙 古 石 油 化 工
１３ ９
有 限 差 分 法 在 低 渗 油 藏 数 值 模 拟 中 的 应 用
饶 盛 文 贺 成 才 赵 丰 ， ，
（ ．西 南石 油大学研 究生院 ， １ 成都 ６０ ０ ； ．四川郫县红 兴中学 ， 都 １５０２ 成 ６１３） １ ７０
２
专 ｎ 。 ｃ ・ 川 一［
。 ・ ＋ ｎ
以上 各式 中 ， 标 ｏ Ｗ 分 别 为 油 和 水 ；一 渗 流 下 、 口 速度 ， ｓ ｍ／； 渗 透 率 ， ；， 相 对 渗 透 率 ， 数 ； 一 ｍ。 一 小
一
。 。 ＋三 。 。 “ ＝ ｎ 户 ］ ｎ 一
・
） 一 ・ ｔ 一 ａ
譬一 ［］
（） １ １
饱 和度 方 程
如＋ Ｓ 一 １ （） ５
数值 计 算采 用 Ｉ ＭＰＥ Ｓ方 法 ， 即对 方 程 （ ） 用 ８采
隐格式 ， （） 将 ２ 代入 方 程 （） 采用 显格 式 ， 有 ４并 则
毛 管 力方程
Ｐ 一 户 一 Ｐ （ 。 ５ ） （） ６
低渗 透 油 藏 ， 由于 其 储 层 致 密 ， 透 率极 低 ， 渗 其 渗 流不 再 满足 达西 定律 。石 油 中常含 有数 量不 等 的 氧 化物 ， 低速 渗 流往 往 会伴 随一 些物 理 化学 现 象 ， 在
有限差分法初步
• 引言 • 有限差分法的原理 • 有限差分法的应用 • 有限差分法的实现 • 有限差分法的优缺点 • 结论与展望
01
引言
有限差分法的定义
有限差分法是一种数值计算方法，通 过将偏微分方程离散化为差分方程， 从而求解偏微分方程的近似解。
近似表示微 分，从而将微分方程转化为差分方程。
有限差分法。
COMSOL Multiphysics实现
COMSOL Multiphysics是一款基于有限元法的多物理场仿真软件，也支持有限差分法。 COMSOL提供了友好的用户界面和丰富的物理模型库，使得有限差分法的实现更加便
捷。
有限差分法的并行计算实现
MPI实现
MPI（Message Passing Interface）是一种并行计算的标准，支持多个处理 器之间的通信。通过MPI，可以实现有限差分法的并行计算，提高计算效率。
自适应网格技术
根据解的特性自适应地调整离散点间距，以 提高计算精度和效率。
并行化与优化
通过并行计算和算法优化等技术提高有限差 分法的计算效率。
与其他方法的结合
将有限差分法与其他数值方法或物理模型相 结合，以处理更复杂的问题。
06
结论与展望
结论
01
有限差分法是一种数值计算方 法，通过离散化连续问题为差 分方程，进而求解数值近似解 。
有限差分法原理简单，易于理解和实现，不需要复杂的数学工 具。
有限差分法可以方便地进行并行计算，提高计算效率。
有限差分法可以应用于各种不同类型的偏微分方程，具有广泛 的适用性。
有限差分法的缺点
精度问题
由于有限差分法是一种离散化方法，其精度受到离散点间距的限制， 可能导致计算结果不够精确。
有限差分法及热传导数值计算
其中式(4-20)是绝热边界的一种离散方式，在确定t1(i+1)之值时需要用到t-1(i) 。根据对称性该值等于t2(i)。这样，从已知的初始分布t0出发，利用式(4-17)及(4-19)可以依次求得第二时层、第三时层直到 i 时层上的温度值(见图4-8)。至于空间步长Δx及时间步长Δτ的选取，原则上步长越小，计算结果越接近于精确解，但所需的计算机内存及计算时间则大大增加。此外，Δx及Δτ之间的关系还受到显式格式稳定性的影响。下面我们从离散方程的结构来分析，说明稳定性限制的物理意义。
理论解 在规定的边界条件下积分，有很大局限性； 数值解 借助计算机，前景广阔。
2.1 导热问题数值解法的基本思想——离散化
以有限差分 无限微分 无限划分 实质 达到精度 以差分代数方程 微分方程 计算机帮助 （当离散点足够多时可以满足要求）
据导热问题的控制方程 ( 导热微分方程 ) 若 △x=△y 则有 得
x
y
如图所示 边界节点 (m,n) 只能代表半个元体，若边界上有向该元体传递的热流密度为q ，据能量守恒定律对该元体有： 边界节点离散方程的建立： 平直边界上的节点 2.3 边界节点离散方程的建立及代数方程的求解 傅里叶定律
02
设有一三元方程组：
其中 （ i=1,2,3 ； j=1,2,3 ）及 是已知的系数（均不为零）及常数。
采用高斯——赛德尔迭代法的步骤： （1）将三元方程变形为迭式方程：
假设一组解（迭代初场），记为：
并代入迭代方程求得第一 次解 每次计算均用最新值代入。 以新的初场 重复计算，直到相邻两次迭代值之差小于允许值，则称迭代收敛，计算终止。
小结： 导热问题数值求解的基本思想
基本概念 ： 有限差分法
有限差分法原理
有限差分法原理有限差分法（Finite Difference Method）是一种常见的数值计算方法，广泛应用于工程、物理、地质等领域的数值模拟和求解偏微分方程。
它的原理是将连续的微分方程转化为离散的差分方程，通过对网格节点上的数值进行逼近，从而求解微分方程的数值解。
在本文中，我们将介绍有限差分法的基本原理及其在实际问题中的应用。
首先，我们来看一维热传导方程的数值求解。
假设我们要求解一个长为L的均匀材料棒上的温度分布，其热传导方程可以写为：\[ \frac{\partial u}{\partial t} = \alpha\frac{\partial^2 u}{\partial x^2} \]其中，u(x, t)表示位置x上的温度分布，t表示时间，α为热扩散系数。
为了使用有限差分法求解这个方程，我们需要将空间和时间进行离散化。
假设我们在空间上取N个网格点，将材料棒分为N个小区间，每个小区间的长度为Δx。
在时间上也进行离散化，取时间步长为Δt。
这样，我们可以用u_i^n来表示位置为x_i的温度在时间t_n的值。
将热传导方程在离散点上进行近似，我们可以得到如下的差分格式：\[ \frac{u_i^{n+1} u_i^n}{\Delta t} = \alpha\frac{u_{i+1}^n 2u_i^n + u_{i-1}^n}{(\Delta x)^2} \]通过对时间和空间上的离散点进行迭代计算，我们可以逐步求解出温度在空间上的分布随时间的演化。
这就是有限差分法的基本原理。
除了一维热传导方程，有限差分法还可以应用于更加复杂的偏微分方程，比如二维热传导方程、波动方程、扩散方程等。
在这些情况下，我们需要在空间上取二维甚至三维的网格点，并相应地修改差分格式。
有限差分法的优点在于它简单易实现，而且可以直接应用于一般的偏微分方程，因此在实际工程和科学计算中得到了广泛的应用。
需要指出的是，有限差分法也有一些局限性。
第9讲-有限差分法(二)
i, j i, j
4
1 n 1 2 ( in1, j in, j 1 in h Fi , j ) 1, j i , j 1
ω =1，就退回到高斯-赛德尔迭代法， ω>2，迭代过程变得极 其不稳定，通常1<ω<2，能提高收敛速度。 ω的最佳选择与具 体问题和离散化的情况有关。 对第一类边值问题： 若正方形场域由正方形网格剖分（每 边节点数为p+1），则最佳收敛因子 若长方形场域由正方形网格剖分（两 边节点数分别为p，q，且都大于 15），则最佳收敛因子
离散方程组的系数一般都是有规律的，且各个方 程都很简单，包含的项数不多（取决于差分计算格 式，每个方程待求量的项数不超过5项）。
3
各离散节点上的方程组形式 （节点顺序按坐标先从y轴增 加，再x轴增加，即从下到 上，从左到右，先列后行） 排列。
2 4 h Fi , j 1 i 1, j 1 i, j i 1, j 1 i , j 2 i , j 1 2 i 1, j i , j 1 i 1, j i , j 1 4 i , j h Fi , j 2 i 1, j 1 i , j 2 i 1, j 1 i , j 4 i , j 1 h Fi , j 1
1 4 i 2,... M, j 2,...q 1,q 1,... N 1 4 j 2,...q 1,q 1,... N
对称轴
1, j (2 2, j 1, j 1 1, j 1 )
介质分界面 i ,q 1 (i 1,q 2 i ,q1 i 1,q 2 r i ,q1 )
15
有限差分法边界条件
有限差分法边界条件
《有限差分法边界条件》
嘿呀，今天咱来聊聊这个有限差分法边界条件。
就好像我那次去海边玩一样，这有限差分法边界条件就像是海边的边界。
你想啊，海水有它流动的范围吧，那岸边就是它的一个边界。
而有限差分法边界条件呢，也是给计算划定了一个范围。
就像在海边，我们知道不能超出那个岸边，不然就掉海里啦。
比如说在计算的时候，我们要确定一些边界上的数值，这就好比在海边我们要知道哪里是安全的，哪里不能再往前走了。
如果边界条件没搞清楚，那计算就可能出大乱子，就像在海边不小心走到深水区，那可就危险咯。
我记得那次在海边，我沿着沙滩走啊走，一直走到一个礁石那里，那礁石就像是一个明确的边界。
我不能再往前了，得绕过去或者停下来。
这和有限差分法边界条件是一个道理呀，到了特定的边界，就得按照规则来处理。
而且这边界条件还分不同类型呢，就像海边有不同的地形一样。
有的边界可能是固定的，就像坚固的悬崖；有的边界可能是变化的，就像海浪不断拍打的地方。
咱在研究有限差分法边界条件的时候，也得像在海边探索一样小心翼翼。
要把每个边界都搞清楚，不然就可能在计算的海洋里迷失方向啦。
总之呢，有限差分法边界条件虽然听起来有点专业有点复杂，但其实就和我们生活中的很多边界一样。
就像那次海边之旅让我明白边界的重要性，研究有限差分法边界条件也是为了让我们的计算更准确、更可靠呀。
下次再看到有限差分法边界条件，就想想那片广阔的大海和海边的边界吧，哈哈！。
有限差分法
h
饱和带
非饱和带
M z
隔水底板
30
MODFLOW模块及其作用
**蒸散(EVT)模块
潜水蒸发
可暂且处理为水面 蒸发
蒸发强度随 潜水埋深的 加大而减弱
ET
z
ET = 0
31
MODFLOW模块及其作用
**河流边界(RIV)模块 通过河床与地表水交换 河流水位必须已知 简化处理
地下水排泄：q= C(Hcell Hriv)
2
Wi , j 0
物理意义是什么？
Ti 1 , j
2
H i 1, j H i , j x
2
y Ti 1 , j
2
H i 1, j H i , j x
2
y x Wi , j xy 0
Ti , j 1
H i , j 1 H i , j y
2 2 2
2
如果均质，无补给，无开采
Hi, j
H i 1, j H i 1, j H i , j 1 H i , j 1 4
T的计算
算术平均：
Ti 1 , j
2
Ti , j Ti 1, j 2
; Ti 1 , j
2
Ti, j Ti 1, j 2
几何平均：
H i 1, j H i , j 1 均质各向同性，四边都为定水头边界 求H1、H2、H3和H4
4H1 H 2 H 3 1 0 4H 2 H1 H 4 1 0 4H 3 H1 H 4 1 1 4H 4 H 2 H 3 1 1
K Aside
q=C(H0Hcell); C=
28
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
有限差分法不同边界条件下的数值模拟
文章介绍了地震数据处理中所使用的数值模拟法，对采用有限差分法所使用不同边界条件处理方式进行了数值模拟，通过波场快照直观的得出了不同的边界吸收条件的吸收效果，对结果进行了对比，分析总结了各种方法的优缺点。
标签：数值模拟；有限差分；边界条件
随着近年来国家宏观经济调控，经济增长的速度逐步减缓，能源行业受此影响最为严重，许多煤矿是在亏损的情况下生产，直接导致了地质行业投入的减少。
物探行业压力也越来越大，物探行业应该抓紧发展先进技术，提高能源勘探的效率。
在物探行业中，地震勘探作为一个重要的手段，发挥着巨大的作用。
数据处理作为地震勘探的一大重要环节，所采用的各种方法和技术手段也一直在更新和进步。
在地震勘探处理方法研究中，地震数值模拟技术可以在室内完成地震数据模型的建立，并对其地震数据进行各种方法的处理，查看处理方法的效果和数据的好坏，另一方面，地震数值模拟进行正演获得的数据也可以作为反演的基础进行比对。
在地震数据处理的过程中，如何模拟地震波的传播便是需要解决的问题。
在二十世纪70年代开始采用显示差分格式来模拟地震波的传播。
由于有限差分法适用条件广，计算速度比较快，占用计算机内存少，编程比较容易实现，模拟精度相对较高而得到广泛应用。
但是有限差分法模拟地震波场时，由于计算机运算核心的限制，有限差分方法只能得到有限的数据点，地震波动方程只能是在有限差分方程中求得近似解，这时就考虑到人工边界问题，如果不对边界进行处理，波在通过边界时会产生反射，因此我们希望对添加的边界进行处理来消除这些反射。
在20世纪70年代，地球物理学界陆续采用了不同的边界条件来实现削弱地震波在通过边界时的反射，比如reynold边界、clayton边界、cerjan边界，以及后来提出的PML层边界条件，每种边界条件都在不同程度上实现了地震波通过边界时的衰减。
为了验证以上边界条件在数值模型的效果，在文章中，我们设计了一些简单的数值模型，给出了不同的边界条件，通过波形在通过不同边界条件时反射进行比较，观察每种方法衰减反射的效果。
1 reynold边界条件
文章首先采用的边界条件为reynold边界条件。
这是ALBERT C.REYNOLDS 于1978年提出的一种边界条件的处理方式。
地下介质中地震波的传播规律可以
近似地用声波方程描述，对于二维模型，首先我们给出二维情况下声波方程，
在此方程中S（t）为震源函数，使用了雷克子波，其中时间域采样间隔为0.001s，采样点数为100点，总的时间长度为0.1s，则频率域采样间隔为10hz。
v（x，z）为介质在点（x，z）处的纵波速度，u为描述速度位或者压力的波场。
震源函数波形如图1：
建立的数值模型是一个均匀速度模型，其中网格数为300x300，速度vp=3000，其中上边界采用自由反射边界，其余三个边界采用reynold边界，可以由图2得到的波场快照看出波在通过边界时反射得到了有效衰减。
2 clayton and engquist边界条件
采用同样的模型参数，当下边界使用自由反射界面，其余三面采用clayton and engquist边界条件，得到的波场快照如图3所示。
两种边界条件得到的单炮记录如图4所示。
由所得到的波场快照可以直观的看出，两种边界处理方法都有明显效果，当地震波通過未经过处理的边界时，地震波的反射现象十分严重，产生一系列反射波并多次反射，对数据形成干扰，但是经过处理的边界时反射得到了明显的削弱。
3 Cerjan吸收边界条件
下面所使用的数值模型是时间二阶、空间二阶精度二维弹性波模拟，采用的震源仍然是雷克子波，其时间域采样间隔为0.001s，主频为30。
模型采用300x300的网格均匀常速度介质，其中P波速度Vp=2000，S波速度Vs=1400，边界条件为上边界采用自由反射边界，其余边界采用Cerjan边界条件。
图5为全波震源得到的波场快照，上图为x分量，下图为z分量。
由波场快照可以看出，地震波在传播时，由于纵波速度比横波速度快，所产生的波形为两个同心圆，外层为直达纵波，内层为直达横波。
在经过边界时，由于上边界采用了自由边界条件，可以看出地震波的反射现象也是比较明显，其他边界由于添加了衰减吸收函数，所产生的反射比较微弱。
同时由于该吸收采用的衰减函数比较简单，计算时所需要的时候也比较少，计算起来速度比较快。
4 结束语
由以上图可以看出，在进行地震模型数值模拟时，所采用的边界条件都达到了一定了效果，但是也有一些不足，采用reynold边界条件和Clayton条件时，其边界吸收对比自由边界效果比较明显，但是仍然具有一定的反射，采用cerjan 边界时，波在边界反射后产生的频散也比较严重，除了以上方法之外，PML边界条件是近年来提出的吸收效果比较好的一种边界处理条件，其原理是在构建的地质模型外围添加一层吸收边界，这种边界能迅速吸收边界中的反射，效果比较
好，所需要的计算量也比较少，使用也比较普遍。
由于采用这种边界的模型文献较多，在此不再进行叙述。
参考文献
[1]贾跃玮.基于MATLAB的地震正演模型实现[J].石油工业计算机应用，2009，2：10-13+65.
[2]左莹.基于高阶交错网格的有限差分地震波场数值模拟[D].长安大学，2009.
[3]王开燕，周妍，刘丹，等.有限差分波动方程正演模拟中的吸收边界条件[J].当代化工，2014，5：752-755.
[4]李桂花，朱光明，马德堂.井间地震数据处理方法——以大庆某采油厂井间地震勘探为例[J].石油地球物理勘探，2011，1：76-82+164+171.
[5]朱光明，李桂花，程建远.煤矿巷道内地震勘探的数值模拟[J].煤炭学报，2008，11：1263-1267.
作者简介：韩浩（1992-），男，山东济宁人，在读研究生，专业为地球探测与信息技术，主要从事地球物理勘探数据处理方法的研究。