天然气工艺流程图_天然气动态建模过程描述怎么写的过程

1.Jason反演技术在天然气水合物速度分析中的应用

2.Simpleware数字岩心建模与数值分析解决方案，求大神详细介绍

3.我想问一下关于石油工程的课程设置，既整个专业包括的所有课程？尽可能详细一些，帮帮忙

4.三维地质建模方法

能源终端利用系统集成建模与优化

工业能量系统集成建模和优化，节能技术的系统集成研究开发；

建筑物能量系统集成优化；

交通系统节能的集成建模和优化研究；

能源--经济结构--生态环境三者协调互动的全生命周期（LCA）和全方位评价（GSA）模型和优化理论及方法

传递过程强化与控制理论和能源新材料研究

宏观与微观传热传质过程的强化与控制的场协同理论；

多种介质、多种相态、及多种过程耦合的传热传质强化；

新一代强化型的管壳式换热器管、壳程各种强化技术；

节能、储能和能量高效转换新材料的制备及关键性技术；

天然气能源利用技术和策略

天然气能源优化利用技术研究；

分布式能源系统集成优化研究；

液化天然气（LNG）的冷能及相关低温设备研究；

城市天然气管网规划及平衡匹配优化研究；

亚热带及沿海可再生能源利用技术

生物质能、沿海近海风力及变化规律调查研究；

太阳能功能材料、太阳能风力互补发电系统、太阳能吸附蓄冷空调以及地热应用的研究；

固体废弃物处理过程机理研究。

Jason反演技术在天然气水合物速度分析中的应用

赵忠泉

（广州海洋地质调查局 广州 510760）

作者简介：赵忠泉，男，（1983—），硕士，主要从事海洋油气调查研究工作，E-mail：zzqhello＠163。

摘要 利用地震数值模拟技术结合实际资料，可以建立各种地质体的地震识别模型，有效地避免地震现象的多解性，从而可以提高解释的精度。本文介绍了二维地质建模的方法流程及两种模拟方法-褶积法和PSPI波动方程法，前者无边界条件约束和频率域中的信号损失，简洁易行，计算稳定，应用广泛，是最早的地震波场模拟方法；后者通过求解波动方程，包含丰富的波场信息，能够充分反映地震波的动力学和运动学特征。实际应用中利用褶积法对三维潮道模型及简化的碳酸盐岩多旋回倾斜薄互层沉积模型进行了模拟；利用零炮检距的频率波数域的波动方程法模拟了生物礁的地震响应，结果对于碳酸盐岩生物礁识别有一定指导意义。

关键词 地质建模 数值模拟 褶积法PSPI法

不同地质体由于其岩性、物性、含油气性、内部结构和岩石组合等的差异，在地震上具有不同的反射特征，包括内部结构、外部形态、振幅、频率等参数。由于地震波在地下地质体中传播的复杂性，加上各种干扰，造成了地震剖面中的各种反射现象存在多解性，大大增加了地震解释的难度。利用地震数值模拟技术结合实际资料，在建立不同地质体的地震识别模型的同时也有效地避免地震现象的多解性，从而可以提高解释的精度。

1 地质建模

地震数值模拟技术的基础是地质地球物理模型的建立，可归结为对地质及地球物理模型结构的数学描述。

二维封闭结构模型用于建立复杂地质模型。二维封闭结构模型就是定义相同地质属性为一独立封闭的地质单元，按照地质属性将地质模型划分成多个独立封闭的地质单元，把所有独立封闭地质单元按照空间分别有序地排列起来，这样组成的集合体就构建了一个二维地质模型。封闭结构模型是以积木方式定义地下地质结构，可以描述非常复杂的地质体。二维封闭结构模型被描述为具有相同地质属性（速度、密度等），并被地层界面、断层界面或模型边界所围成的地质单元的有机组合。对封闭结构模型的描述，实际上就是描述封闭地质单元和封闭地质单元之间的关系，前者包括对封闭地质单元属性和封闭地质单元边界的描述；后者是对地质单元空间关系的描述，也就是描述封闭地质单元边界相接关系及地层属性［1］。

在进行数值模拟过程中，为了验证某些复杂地质体的波场特征，需要绘制多种不同的地质模型，通常可借助常规绘图软件（绘图板、Photoshop、CorelDraw，AutoCAD等）绘制好二维封闭结构面，再根据图像处理中的区域填充算法（填充和扫描转换填充），对不同二维封闭结构面进行不同颜色的填充。其中不同颜色代表不同的二维封闭结构面属性（速度、密度等）；合并相同属性的封闭面，形成最终的二维封闭结构模型［1］。为了得到二维封闭结构模型的属性（速度、密度等）模型，需要对二维封闭结构模型的彩色图进行速度像素空间和属性空间转换，根据颜色空间和属性空间的相互映射，就可以得到复杂地质体的属性（速度、密度等）模型，如图1为模型创建流程图。

图1 二维封闭结构模型建立流程图

2 两种数值模拟方法

2.1 褶积模型

在褶积模型中，我们把地震反射信号s（t）看作是地震子波w（t）与地下反射率r（t）的褶积。地震子波w（t），使用实际地震系统记录到的地下一个单独的反射界面反射的波形（如图2，理想的无噪声褶积过程）。反射率r（t）则代表理想的无噪声地震记录。记录到的地震道s（t）可看作是地震信号w（t）* r（t）与可加噪声n（t）之和，因此可以把地震道看作是一种有噪声干扰的，经过了滤波的地下反射率的变形。

在无噪声褶积模型中，我们把地震信号S（t）看作是地震子波w（t）和地下反射系数r（t）的褶积：

南海地质研究.2012

式中：s（t）——合成地震记录；

r（t）——反射系数；

w（t）——地震子波。

图2 褶积过程

2.2 PSPI波动方程法

通过求解波动方程的数值模拟方法，能够充分反映地震波的动力学和运动学特征，波场信息丰富，模拟结果较为准确。这里仅介绍适合横向速度剧烈变化的频率-波数域相移加插值的波场延拓方法［2］。

相位移加插值的波场延拓方法，简称PSPI法，基本思想是在波场向下延拓的每个深度步长Δz之内，将波场的延拓分成两部进行，首先用L个参考速度V1，V2，…VL，将位于深度zi处的波场p（x，zi，ω）延拓到zi＋1＝zi＋Δz处，得到L个参考波场p1（x，zi＋1，ω），p2（x，zi＋1，ω），…，PL（x，zi＋1，ω）。第二步，按实际的偏移速度V（x，z）同参考速度V1，V2…，VL的关系，用波场插值的方法求出zi＋1处的波场p（x，zi＋1，ω），按同样的步骤，可将zi＋1处的波场值p（x，zi＋1，ω）延拓到深度zi＋2，得p（x，zi＋2，ω），直到延拓到最大的深度zmax为止。

对于各向同性介质，取二维标量声波方程作为延拓的基本方程：

南海地质研究.2012

式中，p＝p（x，z，t）为二维地震波场值；x，z分别为水平方向和垂直方向坐标轴；t为时间轴；v（x，z）为纵、横向都可变的地震波传播速度。将式（2）分别对x、t作傅氏变换，考虑到并考虑到?2/?x2和与（-ikx）2和（iw）2的对应关系，可得：

南海地质研究.2012

式中， 是p（x，z，t）的二维傅氏变换；v为地震波速度；w为圆频率；kx为水平波数；kz为垂直波数。零炮检距情况下的地震记录模拟只考虑单程波，因此可得到相位移波场延拓公式如下：

南海地质研究.2012

式中， （kx，zi，w）为频率波数域波场值；Δz为深度延拓步长；kx为测线方向波数；kz为深度方向波数。式（4）为二维波场正演公式，其延拓方向为由地下向地面延拓；式（5）为二维波场偏移公式，其延拓方向为由地面向地下延拓。

为了适应地下地震波场速度在纵横向均可变的要求，在同一延拓深度内用几个不同地震波速度分别作相移，再用拉格朗日插值公式进行插值，就可求出所有的以不同速度传播的延拓波场值P（x，zi＋1，t），从而近似地解决了横向变速时的波场延拓问题［3］。

3 模拟实例

3.1 三维潮道数值模拟

运用褶积原理建立了一个简单三维潮道模型，此三维潮道事实上为多个（128）二维剖面排列而成，三维模型的样点为128×128×128，利用MATLAB实现。选用子波为雷克（Ricker）子波，其公式为：

南海地质研究.2012

其中fp为主频。在处理过程中选用主频为fp＝40 Hz、样间隔2 ms，对称样点数为24，子波波形如图3。

图3 雷克子波

图4 潮道平面图

图4为潮道平面图，该图仅反映了潮道的平面形态，作为计算机实现三维建模的边界控制，横坐标代表inline线，纵坐标代表xline（crossline）线，图5为三维地质模型示意图，模型较简单，整体由三个水平层叠置而成，在第二层和第三层之间镶嵌了形如图4的潮道，此潮道没有考虑进水方向，根据此地质模型进行计算机地震正演模拟，可得到相应三维地震数据体，从图中可以看到，**虚线（上）和蓝色虚线（下）位置上，分别横跨了三个潮道分支和两个潮道分支，就是说在相应两条虚线位置上的两条测线应该分别有三个和两个潮道显示，提取相应的两条剖面如下图6和图7：

图5 三维地质模型

图6 xline＝100（**虚线）剖面

图7 xline＝100（蓝色虚线）剖面

再在三维数据体中沿水平方向做切片，即提取时间切片。图8为时间切片在地震剖面上的位置示意图，图中五条标示线从上到下依次为白色实线、**虚线、白色实线、红色虚线和白色实线，与之对应的时间分别为70 ms、85 ms、95 ms、99 ms和110 ms（时间范围是0～128 ms），图9～图13为相应切片，从图中可以看出，随着所做切片时间的增大（深度的增加），潮道的展布范围逐渐减小，由于地层是水平层状的，使得时间切片等同于地层切片和沿层切片，其切片效果非常明显，切片中潮道形态得到了很好的展示，但是在多个切片中发现，从可以见到潮道形态一直到潮道消失的时间范围是在70～110 ms之间，而潮道的真实范围是在80～100 ms之间，显然依据切片所圈定的潮道的范围相比真实的范围扩大了，究其原因是由于不管选取哪一波，子波都有一定的延续长度和有限频宽，这就限制了合成地震记录本身的分辨率并不能达到等时厚度反射系数序列的分辨率。因此在对实际地震资料进行解释的时候，对地质异常体边界的识别应该考虑地震子波并非脉冲波所带来的影响。

图8 剖面示意图

图9 切片t＝70 ms

图10 切片t＝85 ms

图11 切片t＝95 ms

图12 切片t＝99 ms

图13 切片t＝110 ms

3.2 薄互层沉积模型

图14为简化的碳酸盐岩多旋回倾斜薄互层沉积模型（Zeng，2003），模型简化是为了更好地突出由岩相控制的波阻抗结构和地震信号之间的相互关系。该模型所有倾斜的倾角都相同，每层都有相同的垂直时间厚度（5 ms或15 m，速度为6000 m/s），泥岩与低孔隙度颗粒灰岩的波阻抗差，以及低孔隙度颗粒灰岩与高孔隙度颗粒灰岩的波阻抗差都相同，所有高孔隙度颗粒灰岩具有相同的深度范围，综合起来形成了一个水平的岩性地层单元。

其时间域地震响应（图15）中，高频情况下（60 Hz雷克子波），地震反射被建设性地调谐到时间地层单元，因此地震同相轴沿着时间地层单元分布（图15a）。当子波频率减到40 Hz时，地震反射对时间地层单元和岩性地层单元都有响应（图15b）。当用30 Hz雷克子波时（图15c），地震同相轴破坏性地调谐到时间地层单元和建设性地调谐到岩性地层单元，因而时间地层单元的反射进一步变弱，地震同相轴被岩相反射所控制［4］。

这个模拟过程强调了了解地质格架和时间地层单元以及岩性地层相带厚度尺度的重要性。时间地层（图15a）和岩性地层（图15c）成像都是有用的，前者用于对比，后者用于粗略的储层评价。然而，这两种响应不能混淆在一起。图15b中的两组相互矛盾的地震同相轴会造成地震象［4］。

图14 简化的碳酸盐岩多旋回倾斜薄互层沉积模型

3.3 生物礁数值模拟［5～7］

频率—波数域的相移加插值偏移（PSPI）在每一个深度间隔内使用多个参考速度进行偏移，由多个偏移结果插值生成最终的偏移剖面，所用插值的速度越多，越能反映实际介质的速度变化情况，此方法在成像精度及横向变速适应性上具有很大的优越性，但处理所需的时间稍长，鉴于本文的二维叠后建模对处理时间没有过高要求，因此应用PSPI方法做正演、偏移。

图16为某区块过生物礁的原始地震剖面，图17为根据此剖面建立的生物礁速度模型：模型速度变化范围是5600 m/s到5980 m/s，从图16中可以看出生物礁的底界面清晰可辨，围岩有披覆现象，内部呈杂乱反射。为了检验该地质建模的正确性，先用PSPI方法对该模型进行了波场正演模拟计算，其模拟剖面如图18所示。由于生物礁埋藏深，生物礁顶底反射的弧度较大，不规则点的绕射波杂乱，因此用图15的速度模型对其进行叠后时间偏移，得到了偏移剖面（图19），横向表示256个地震道，纵向表示零偏移距反射时间，礁体最大时间厚度约40 ms。从图19可以看出，模拟记录中的礁体顶界与原始剖面有一定差距，但是生物礁底界反射和内幕反射以及侧翼反射与原始剖面基本一致，其他的地层界面形态与原始剖面也吻合较好，在一定程度上验证了地质模型的正确性，说明当生物礁与围岩之间存在一定波阻抗差异时，在地震剖面上必然出现异常反映，经过有效的构造和参数反演，能够将其分辨出来。相信通过模型改进以及算法中参数的调整，能够与原始剖面更好地吻合，从而为生物礁的地震解释提供一种有力的验证工具。

图15 图14模型时间域地震响应

4 结论

地震数值模拟（正演）技术基于地球物理模型的建立，运用概念二维封闭结构地质模型的建立方法，得到复杂地质体的数学模型，结合各种算法对其进行模拟从而可以验证相应地质体的地震波场特征；结合实际资料建立不同地质体的地震识别模型，可以有效地减少地震现象的多解性，从而提高解释的精度；褶积法无边界条件约束和频率域中的信号损失，简洁易行，计算稳定，应用广泛，本文用此方法模拟的伪三维潮道模型及倾斜薄互层模型取得了较好的效果；通过求解波动方程的数值模拟方法，包含丰富的波场信息，能够充分反映地震波的动力学和运动学特征，PSPI波场沿拓方法为其中之一，利用正演与偏移相结合的流程模拟了生物礁的地震响应特征，检验解释成果的正确性，为生物礁的地震解释提供了一种有力的检验工具。

图16 原始剖面

图17 生物礁地质速度模型（256×256）

图18 正演记录（子波主频30Hz）

图19 偏移剖面（子波主频30Hz）

参考文献

［1］刘远志.碳酸盐岩地震相分析与数值模拟［D］.成都：成都理工大学，2009.

［2］韩建彦.复杂地质体地震正演与偏移［D］.成都：成都理工大学，2008.

［3］贺振华，王才经等.反射地震资料偏移处理与反演方法［M］.重庆：重庆大学出版社，1989.

［4］Zeng Hongliu ＆Kerans，C.Seismic frequency control on carbonate seismic stratigraphy；a case study of the Kingdom Abosequence，West Texas，American Association of Petroleum Geologists Bulletin，2003.87，273～293.

［5］贺振华，黄德济，文晓涛，等.碳酸盐岩礁滩储层多尺度高精度地震识别技术［R］.成都：成都理工大学地球探测与信息技术教育部重点实验室，2009.

［6］熊晓军，贺振华，黄德济.生物礁地震响应特征的数值模拟［J］.石油学报，2009，30（1）：7～65.

［7］熊忠，贺振华，黄德济.生物礁储层的地震数值模拟与响应特征分析［J］.石油天然气学报，2008，30（1）：75～78

The lication and analysis of two kinds of seismic numerical simulation method based on the2D-geological modeling

Zhao Zhongquan

（Guangzhou Marine Geological Survey，Guangzhou，510760）

Abstract：Pick to using seismic numerical simulation technology combined with the actual seismicdata，we can build all kinds of seismic recognition model of geologic body and effectively oidthe multiple solutions of seismic phenomenon，which can improve the precision of the explana-tion.This paper describes the method of the process of 2D geological modeling and two simulationmethods，seismic convolution method and PSPI we equation method，the former has no bounda-ry condition and the signal loss in frequency domain，is concise and easy，it can be calculatedsteadily and be lied widely，is the earliest simulation method in seismic we field，the latterbased on the we equation，it contains the rich information in we field，can fully reflect thedynamics and kinematics characteristics of seismic we.In the practical lication，we use theconvolution model in 3D-tidal channel model and the multi-cyclic simplified deposition model oftilt thin interbed layer of carbonate；We simulate the seismic response of reefs using the method ofzero-offset we equation in frequency and we number domain，it is confirmed that the resulthas definite significance for the identify of the reef.

Key words：Geological modeling Numerical simulation Convolution PSPI method

Simpleware数字岩心建模与数值分析解决方案，求大神详细介绍

梁劲1　王宏斌1，2　梁金强1

（1.广州海洋地质调查局 广州 510760；2.中国地质大学（北京）北京 100083）

第一作者简介：梁劲，男，11年生，高级工程师，1995年毕业于成都理工学院信息工程与地球物理系应用地球物理专业，主要从事天然气水合物调查与研究工作。

摘要 本文用Jason 反演技术对南海北部陆坡A 测线纵波速度进行计算，结合BSR、振幅空白带以及波形极性反转等多种水合物赋存信息的分析，对水合物成矿带的速度特征进行了综合研究，结果表明：低速背景中的高速异常，是天然气水合物赋存的重要特征；高速异常体一般呈平行于海底的带状分布；在高速异常的内部，速度也是不断变化的。一般在异常体的中心速度最高，由中心到边缘速度逐渐降低，反映在水合物矿带内部，水合物饱和度由矿体中心向边缘逐渐降低的特征。本文的研究成果进一步表明高精度速度分析不仅可以帮助寻找水合物矿点，还可以进一步判定水合物的富集层位。

关键词 Jason 反演技术 天然气水合物 速度分析

1 前言

天然气水合物是在低温、高压环境下，由水的冰晶格架及其间吸附的天然气分子组成的笼状结构化合物，广泛分布于海底和永久冻土带。温度和压力是天然气水合物形成和保存最重要的因素（王宏斌等，2004）。针对天然气水合物的野外调查及研究表明：高分辨率的地震勘探方法是天然气水合物调查评价中行之有效的方法。地震反演技术一直是地震勘探中的一项核心技术，其目的是用地震反射资料反推地下的波阻抗、速度、孔隙度等参数的分布，从而估算含天然气水合物层参数，预测天然气水合物分布状况，为天然气水合物勘探提供可靠的基础资料。常用的地震反演技术有Jason、Strata、Seislog和ISIS等，其中Jason反演技术在含天然气水合物层预测中因其分辨率高而得到广泛推崇，它主要由有井约束和无井约束两种方法组成（廖曦等，2002）。

速度异常是判断天然气水合物是否赋存的重要条件之一。结合BSR（Bottom Simulating Reflector）特征、波形极性特征、振幅特征以及AVO特征等目前已成为判断是否存在天然气水合物层主要手段（史斗等，1999）。大量的测试数据显示：水合物的速度与冰的速度较为接近，而比水高。与含水或含游离气沉积层相比，含水合物沉积层的密度降低，声波速率增大，含水合物层的地层速度往往比一般的地层速度高，含水合物沉积层的下部由于充填了水或气，而使水合物底界面出现速度负异常。因此，地层中速度反转是水合物赋存的一个地球物理标志。含水合物地层的声波速度与水合物的含量有关，水合物含量越高，其声波速度越高。从速度方面看，BSR是上覆高速的含水合物地层与下伏较低速的含水层或含气层之间的分界面。通常，海洋中浅层沉积层的地震纵波速度为1600～1800m/s，如果存在水合物，地震波速度将大幅提高，可达1850～2500m/s，如果水合物层下面为游离气层，则地震波速度可以骤减200～500m/s。因此，在速度剖面上，水合物层的层速度变化趋势呈典型的三段式，即上下小、中间大的异常特征（张光学等，2000）。西伯利亚麦索雅哈气田的资料表明，在原为含水砂层内形成水合物之后，其纵波的传播速度会从1850m/s提高到2700m/s；而在胶结砂岩层，这种速度会从3000m/s提高到3500m/s。深海钻探的570站位的测井结果表明，由含水砂岩层进入含水合物砂岩层时，密度由1.79g/cm3降低到1.19g/cm3，声波传播速度从1700m/s提高到3600m/s，且电导率剧烈下降。

Cascadia海域ODP889站位的VSP测井资料反映水合物底界为强烈的负速度界面，速度从水合物沉积物层的1900m/s陡降到含游离气层的1580m/s，由于VSP测井为地震测井，受钻井因素的影响较少，因此认为VSP测井真实地反映了水合物沉积层底界的速度变化（陈建文等，2004）。

国土部广州海洋地质调查局在2001～2004年在南海北部陆坡进行10000多公里的天然气水合物高分辨地震调查。本研究利用Jason反演技术，通过对南海北部陆坡区的地震速度资料的精细分析，在已圈定BSR分布范围的基础上研究陆坡区各沉积层的速度特征，最后对速度值与水合物的关系进行了分析和探讨。

2 方法原理

纯天然气水合物的密度（0.9g/cm3）和海水密度相近，而游离气的含量又十分有限，这就决定了产生BSR的波阻抗差主要由速度造成。速度反演技术的特点是在无井约束时，以地震解释的层位为控制，对所有的地震同相轴来进行外推内插来完成波阻抗反演，这样就克服了地震分辨率的限制，最佳的逼近了测井分辨率，同时又使反演结果保持了较好的横向连续性。速度反演技术的主要原理是：①通过最大的似然反褶积求得一个具有稀疏特性的反射系数系列；②通过最大的似然反演导出波阻抗；③通过波阻抗计算速度。该方法的主要优点是能获得宽频带的反射系数，是一种基于模型的反演，具有多种建模方法，对所建模型进行比较分析，并使地质模型更趋合理，反演结果更加真实可靠（郝银全等，2004）。

波阻抗反演方法的出发点是认为地下的反射系数是稀疏分布的，即地层反射系数由一系列叠加于高斯背景上的强轴组成。具体反演是从地震道中，根据稀疏的原则抽取反射系数，与子波褶积生成合成地震记录，利用合成地震记录与原始地震道的残差修改反射系数，得到新的反射系数序列，然后再求得波阻抗。其具体步骤是：

设地层的反射系数是较大的反射界面的反射和具有高斯背景的小反射叠加组合而成的，根据这种设导出一个最小的目标函数（安鸿伟等，2002）：

南海地质研究.2006

式中：R（K）为第一个样点的反射系数，M为反射层数，L为样总数，N为噪音变量的平方根，λ为给定反射系数的似然值。

最大的似然反演就是通过转换反射系数导出宽带波阻抗的过程。如果从最大的似然反褶积中求得的反射系数式R（t），则波阻抗：

Z（i）=z（i-1）×（1+R（i））/R（1-i） （2）

利用波阻抗和速度的关系式：

v=Z（i）/ρ （3）

即可得到速度值。其中，ρ为地层密度，可从区域测井资料结合该测线重力资料反演求取。

在上述过程中为了得到可靠的反射系数估算值，可以单独输入波阻抗信息作为约束条件，以求得最合理的速度模型。一方面，速度反演结果是一个宽频带的反射序列和波阻抗及速度数据，同时加入了低频分量，使反演结果更能正确反映速度变化规律；另一方面，它有多种质量控制方法，具体表现为监控子波的选取、同相轴的连续追踪、反演结果准确性的判断和提供多种交汇显示的相关性分析。所以利用速度反演可对地震剖面上任一相位进行速度反演，在每一个CDP点都可得到任一个同相轴速度数据，并利用二维的反射波的速度层析成像反演方法得到高度连续的速度剖面，如果地震测线足够密，还可利用三维速度反演得到速度体图像。

3 实现过程

3.1 初始模型的确立

在地质规律的指导下，利用地震和测井资料开展沉积特征分析和沉积旋回划分；建立岩石-电性关系，进行砂层组和单砂层对比；在地震剖面上提取各含油砂层组反射波属性，建立地震属与矿体的关系，实现地震-测井综合预测矿体平面分布厚度，开展层间矿体组外推预测；建立初始速度场；在地震属性约束下开展地震反演，反演层间小层矿体厚度。细分层反演层位的标定正确与否直接影响反演结果的精度。因此，在反演过程中对子波提取、能谱特点、信噪比、频谱及反射系数的研究至关重要（闫奎邦等，2004）。技术路线流程如图1所示：

3.2 初始速度场的获得

初始速度场的获得首先要对速度谱进行解释，速度谱的解释和取值是否合理，将直接影响均方根速度的计算精度。具体步骤如下：

1）速度谱的解释先从地质条件简单、反射层质量好、能量团强、干扰少的剖面段开始，绘制叠加速度-反射时间曲线，并逐渐向外扩展；

2）结合地震剖面的反射特征，判断速度极值点是否正确，并选择读取能量团最大的极值点。排除干扰波能量团，从而求得有效波的叠加速度；

3）对相邻速度谱进行比较，通过比较速度谱曲线的形状、相同反射层的速度极值等方法予以检查和修改。

4）每隔40个CDP拾取一组数据，利用地震剖面上的反射倾角数据对它们进行校正，便可得到均方根速度（梁劲等，2006）。

图1 速度反演技术线路流程图

Fig.1 The flow chart of the velocity inversion of technical route

3.3 子波的提取

子波提取时，要使能量集中于子波的主瓣，与地震子波形态吻合。如果所提子波近于零相位，则从波峰向两侧能量衰减较快，波峰两侧波形对称；在子波的能谱特征分析，要使能量都集中在地震波的主频范围内；有井资料时，要对井资料都作了子波与地震波自动关联质量控制。保证子波能谱与地震波能谱相吻合，是反演中较为重要的一方面，子波能谱的峰值与地震波主频的能谱峰值相吻合。首先了解合成记录与地震记录之间的偏差。通过合成记录与地震记录之间的偏差分析，对Jason反射系数偏差、能谱偏差进行进一步的校正，使合成记录与地震记录之间的偏差减小。然后通过反射系数与地震资料之间偏差分析，取相应的手段校正，使地层与合成记录反射系数相吻合。再进行信噪析，使反演处理后的信噪比得到最大限度的提高。通过一系列质量控制手段，使各油层合成记录与地震记录的标定精度得到了较大的提高。

关于速度反演可信程度，不能完全由反演方法确定，关键在于获取地震记录的质量和反演前处理流程的振幅保真度。另一个影响因素是数值模拟结果应当是比较准确的，这与计算方法有关，也与子波拾取和地质构造模型有关。至于反演结果的灵敏度，主要由拟合误差值和收敛速度来判断。如果给定的初始模型正确，即与实际地质结构一致，则拟合的误差较小且收敛速度快。本文工作由于受实际情况限制，没有实际的测井资料验证，因此反演所得速度的准确性和精度会受到一定程度的影响。

4 速度剖面特征

运用多种特殊地震成像综合分析，是天然气水合物地震资料解释的关键技术。目前一般用识别BSR、振幅空白带、波形极性反转、速度异常、波阻抗面貌和AVO等天然气水合物地震相应特征来综合分析沉积物中是否含有水合物。高精度的层速度分析可帮助判定水合物的富集层位，速度及振幅异常结构是水合物与下伏游离气共同作用形成的特殊影像，剖面上表现为“上隆下坳”结构，多层叠合构成一明显的垂向“亮斑”这一特殊成像结构在未变形的水合物盆地内较适用于寻找水合物矿点，并可据此定量估算水合物盆地内水合物的数量，分析BSR上下的详细速度结构，是水合物地震资料综合解释的重要手段（张光学等，2003）。

图2 南海北部陆坡测线A道积分剖面

Fig.2 Trace integration profile of the line A in north slope of the South China Sea

图2是南海北部陆坡测线A的地震反射道积分剖面，从图中可以看出，该剖面中部及右下角距海底大约350ms处出现一强振幅反射波，大致与海底反射波平行，与地层斜交，BSR特征明显。在波形极性方面，海底反射波和BSR都表现为成对出现的强振幅双峰波形特征，海底反射波表现为蓝红蓝特征，而BSR表现为红蓝红特征，这表明相对于海底，BSR显示出负极性反射同相轴，即所谓的极性反转（与海底反射相反）。反射波的极性是由反射界面的反射系数决定的，而反射系数则与界面两侧的波阻抗差有关。实际上，海底和BSR都是一个强波阻抗面，海底是海水和表层沉积物的分界面，上部为低速层，下部为相对高速层，反射系数为正值；BSR是含水合物层与下部地层（或含气层）的分界面，上部为高速层（水合物成矿带是相对高速体），下部为相对低速层（如含游离气，则速度更低），反射系数为负值，因此造成了BSR和海底反射波的极性相反现象（沙志彬等，2003）。图3是用速度反演法反演出来的纵波速度剖面，该速度剖面明显显示出一近似平行于海底的相对高速地质体，其位置恰好在BSR上方。高速地质体的纵波速度大约在2000～2400m/s，其上面的低速层的纵波速度大约在1500～1800m/s，而下面的低速层的纵波速度大约在1500～1900m/s，没有明显的游离气存在特征，但根据其高速地质体特征、BSR以及波形极性反转分析，可以认为南海北部陆坡测线A的相对高速地质体极可能是水合物成矿带。

图3 用速度反演法计算的南海北部陆坡测线A纵波速度剖面

Fig.3 P velocity profile of the line A in north slope of the South China Sea computed by velocity inversion

由图3可见，水合物成矿带内部速度是变化的，表明水合物分布不均匀，呈平行于海底的带状分布，中心速度最高，由中心到边缘速度逐渐降低。海底以下有3个近似平行海底的低速和高速带：①海底与高速体之间的相对低速带，为水饱和带；②水合物成矿带；③水合物成矿带下的低速带。水合物成矿带下面的低速带在速度剖面上没有明显的低速特征，由此推断水合物成矿带下可能不含游离气，或者是气体的饱和度很低。

5 结论

水合物的生成除了需要一定的温度和压力条件外，还需要大量的碳氢气体和充足的水。这就需要地层具有较高的孔隙度和渗透率。未固结沉积岩的孔隙度很高，渗透率大，具备水合物生成的物理条件。具备这种特征的未固结沉积岩的地震波速度较低，而含水合物地层的地震波速度增大。这就形成了水合物成矿带作为低速背景中的高速地质体特征。另外，水合物的生成受温度和压力控制，一般情况，等温面和等压面近似平行于海底，因此低速背景中近似平行于海底的相对高速地质体是水合物成矿带的特征（刘学伟等，2003）。

通过对南海北部陆坡A测线纵波速度的计算，并且结合BSR和振幅空白带识别以及波形极性反转等多种特殊地震成像进行综合分析，我们可以进一步了解水合物成矿带的速度特征：揭示水合物成矿带的高速异常一般呈平行于海底的带状分布，在高速异常的内部，速度也是不断变化的，一般在异常体的中心速度最高，由中心到边缘速度逐渐降低，该现象反映在水合物矿带内部，水合物分布并不均匀，水合物饱和度由矿体中心向边缘逐渐降低。分析BSR上下的详细速度结构，是水合物地震资料综合解释的重要手段。高精度速度分析可帮助判定水合物的富集层位，较适用于寻找水合物矿点，并可据此估算水合物量。

参考文献

安鸿伟，李正文，李仁甫，等.2002.稀疏脉冲波阻抗反演在XY油田开发中的应用.石油物探，41（1）：56～60

陈建文，闫桂京，吴志强，等.2004.天然气水合物的地球物理识别标志.海洋地质动态，6：9～12

郝银全，潘懋，李忠权.2004.Jason多井约束反演技术在油气储层预测中的应用.成都理工大学学报，31（3）：2～300

梁劲，王宏斌，郭依群.2006.南海北部陆坡天然气水合物的地震速度研究［J］.现代地质，20（1）：123～129

廖曦，马波，沈浩，等.2002.应用Jason软件进行砂体及含气性预测.天然气勘探与开发，25（3）：34～42

刘学伟，李敏锋，张聿文，等.2005.天然气水合物地震响应研究——中国南海HD152测线应用实例.现代地质，19（1）：33～38

沙志彬，杨木壮，梁金强，等.2003.BSR的反射波特征及其对天然气水合物识别的应用.南海地质研究，15（1）：55～61

史斗，郑军卫.1999.世界天然气水合物研究开发现状和前景.地球科学进展，14：330～339

王宏斌，梁劲，龚跃华，等.2005.基于天然气水合物地震数据计算南海北部陆坡海底热流.现代地质，19（1）：67～73

闫奎邦，李冬梅，吴小泉.2004.Jason反演技术在岩性识别中的应用.石油物探，43（1）：54～58

张光学，黄永样，陈邦彦，主编.2003.海域天然气水合物地震学.北京：海洋出版社

张光学，文鹏飞.2000.南海甲烷水合物的地震特征研究，首届广东青年科学家论坛论文集，中国科学技术出版社

The Application of Jason Inversion Technology in Velocity Analysis of Gas hydrate

Liang Jin1　Wang Hongbin1，2　Liang Jinqiang1

（1.Guangzhou Marine Geological Survey，Guangzhou，5107602.China University of Geosciences（Beijing），Beijing，100083）

Abstract：The P velocity of A seismic profile in the north slope of the South China Sea were calculated by Jason inversion method.The velocity characterostic of the gas hydrate bed was researched in detail based on the calculated result and the information of gas hydrate existing including BSR，amplitude blanking and polarity reversion of the weform.Research shows that：The abnormity of higher velocity in the background of lower velocity is an important characteristic of gas hydrate existing；The abnormity of higher velocity which distribute as a belt usually parallel to the seafloor；The velocity changes gradually at the inner of the abnormity of higher velocity with the highest velocity at the center of the abnormity whereas the lowest velocity at the margin of it，which suggests that the saturation of gas hydrate decreases gradually from the center to the margin.The result that mentioned above suggest that high resolution velocity analysis not only help to search the hydrate spot but also help to estimate the rich layer of gas hydrate.

Key Words：Jason Inversion Technology Gas hydrate Velocity Analysis

我想问一下关于石油工程的课程设置，既整个专业包括的所有课程？尽可能详细一些，帮帮忙

案例展示

1、利用X射线分析土壤团聚体的层析和水的流动模拟

土壤团聚体是土壤的结构单元，它产生复杂的孔隙系统控制土壤水汽储存和通量。在膨胀和收缩或外力作用下，可以与机械压实一样破坏聚集体，物理冲击如何改变骨料结构的认知仍然有限。本研究的目的是量化压实对团聚体的影响，主要对孔径分布及流动。三维模型建模用Simpleware软件完成。本案例来自M. Menon a,*, X. Jia等发表在Soil & Tillage Research上的Analysing the impact of compaction of soil aggregatesusing X-ray microtomography and water flow simulations.

压实前后的三维模型对比

压实前后流速的仿真分析

2、多孔介质流动的孔隙尺度模拟

3、基于碳酸盐岩图像的三维孔隙尺度流动模拟

通过CMT技术对三种岩石做X射线扫描后成像，利用Simpleware软件的ScanIP模块对其做数字图像处理，得到三个不同的岩石三维模型，之后利用有限元软件分别做氮气、丙烷、汽油、水以及发动机机油做三维孔隙流动模拟，得到其不同的流动特性。

4、循环和静态载荷作用下岩石断裂损伤过程区（FPZ）的实验与数值模拟研究

目前,很多国家的地下开釆都已经或即将进入了深部开阶段，并且，地下矿逐渐走向深部开釆是趋势。未来10多年的时间内，我国有1/3的有色金属矿山的开深度将达到或超过1000m。

在深部，矿岩处于高原岩应力、高地温、高岩溶水压的复杂环境下，表现出与浅部截然不同的力学特性，典型的如岩爆和深部分区破坏裂化现象等。因此，为解决深部开中出现的问题，多国矿山企业和以及相应的国际组织都对深部高应力岩体的力学问题开展了研究，中仿科技的岩土重要合作院校如中科院武汉岩土所、东北大学、中国矿业大学、中南大学、重庆大学、北京科技大学等开展了大量关于深部工程的研究和实践，这些工作对深部岩石力学的发展起到了非常重要的作用。

昆士兰大学土木工程学院的M. Ghamgosar等人用巴西试验(CCNBD)对布里斯班凝灰岩标本的断裂损伤过程区的发展（FPZ）进行了研究，探讨各种静态和循环荷载作用下脆性岩石的力学行为。本案例用扩展有限元法（XFEM）计算裂纹尖端断裂过程区对的位移、拉伸应力分布和塑性耗能等。其中一个最重要的步骤就是用Simpleware软件重构数字岩石模型，还原裂纹尖端断裂过程区延伸的形状，之后用数值模拟，结果与实验非常吻合。

5、类似案例

关于Simpleware软件

借助其image-to-mesh技术，Simpleware三维数字图像重构软件已成为图像到数值模型的图像处理先驱者，获得了包括Queen’sAward for Enterprise in Innovation 2012、Institute of Physics’ (IOP) Innovation Award 2013在内的多个国际奖项，为数字图像三维建模的发展做出了重要贡献。目前Simpleware在世界范围内广泛应用于生物医学、材料科学、石油天然气科学、3D打印等众多领域。

关于中仿科技

中仿科技(CnTech)公司成立于2003年，是中国领先的仿真分析软件和系统解决方案的提供者。中仿科技依靠自主创新研发拥有自主知识产权的中仿CAE系列产品，同时与国际上领先的数值仿真技术公司有长期而紧密的合作关系，具备较强的自主研发能力和创新能力，能够为中国企业和科研机构提供世界一流的仿真技术解决方案。公司总部设在上海，目前在北京、武汉设有分公司。

过去的十多年来，中仿科技一直致力于仿真技术领域最专业的系统实施和项目咨询。目前在中国已有超过1500家用户，其中包括中国航天、中国商飞、中石化、中海油、交通部、地震局、国家电网、中广核以及各大高校和中科院所。服务领域涉及高端制造、国防军工、石油化工、水利水电、汽车交通、能源矿、生物医学、教学科研等。

三维地质建模方法

石油工程的主干学科是石油与天然气工程，主要课程有：技术经济学、油气田开发地质、工程力学、计算机程序设计、流体力学、渗流力学、钻井工程、油工程、油藏工程、油田化学、钻新技术等。主要实践性教学环节包括普通地质实习、金工实习、生产实习、毕业实习、毕业设计等（一般安排30周）。

相近专业有矿工程 /石油工/ 矿物加工工程/ 勘察技术与工程/ 勘察工程 /地质工程/ 矿物工程 /油气储运工程/ 煤及煤层气工程 /勘查工程等。

下边我详细解释一下石 油 与 天 然 气 工 程.

石油与天然气工程是研究石油与天然气勘探、评估、开、油气分离、输送理论和技术的工程领域。其工程硕士学位授权单位培养从事石油与天然气生成环境、勘探、油气井工程设计、测井数据集和处理、油气田开、油气储运以及工程管理的高级技术人才。研修的主要课程有：政治理论课、外语课、工程数学、弹塑性力学、计算机应用技术、高等流体力学、高等渗流力学、油藏数值模拟、油田化学、收率原理、现代油气勘探技术、现代油气井工程、现代凿井工程、天然气工程、高等油藏工程、高等油工程、高等输油管道工程、高等输气管工程、油气田输系统、油气管道运行模拟、天然气液化技术、高等管理学基础、能源经济等。

石油与天然气工程是一个运用科学的理论、方法、技术与装备高效地钻探地下油气、最大限度并经济有效地将地层中的油气开到地面，安全地将油气分离、计量与输运的工程技术领域。石油与天然气作为人类社会能源的重要组成部分，由于其不可替代性和自身的不可再生性，在世界经济的发展、人类社会生活与文明中占有极其重要的地位。由于石油与天然气存在着储层埋藏深，物性有低渗、超低渗，油品有稠油、超稠油，加之高压高温、地层非均质、井眼形成难等特点，给钻探与开发增加了很大的困难。目前，我国石油与天然气收率还比较低、地质条件复杂，深井与超深井钻探与开成本还比较高，因此是一项高投入、高风险、但效益明显的产业。在我国，2l世纪将是石油与天然气工程得以迅速发展的时代。

石油与天然气工程涉及工程力学、流体力学、油气地质、渗流物理、自控理论、计算机技术等基础和应用学科，需要解决的工程问题有钻井、完井、测试、油气藏开发地质、油气渗流规律、油气田开发方案与开技术、提高收率、油气矿场收集处理、长距离输送、储存与联网输配等工程问题。本工程领域与矿产普查与勘探、地球探测与信息技术、矿工程、工程力学、化学工程、机械工程、交通运输工程等学科相关。

培养目标是培养从事石油与天然气工程领域所属油气井工程、油气田开发工程、油气储运工程中科技攻关、技术开发、工程设计与施工、工程规划与管理的高层次人才。

石油与天然气工程领域工程硕士应具有本工程领域坚实的基础理论和宽广的专业知识及管理知识，掌握解决工程问题的先进方法和现代化技术手段，具有独立担负工程技术或工程管理工作的能力以及解决工程实际问题的能力，具有较好的综合素质和较强的创新能力和适应能力。掌握一门外语，能较熟练地使用计算机。

领域范围有以下几个方面：1.油气井工程：油气井工程力学，油气井工作液的化学和力学，油气井工程测量与过程控制，油气井测井数据集、处理与解释。2.油气田开发工程：油气藏描述及开发地质建模的理论与方法，渗流理论和油气藏数值模拟，油气田开发理论与方法，油气工程理论与技术，提高收率理论与技术，油气化学工程与理论。3.油气储运工程：油气长距离管输技术，多相管流及油气田集输和油气处理技术，油气储运及营销系统优化，油气管道和储罐的强度研究，油气储运设施施工及安全、防腐技术。4.石油与天然气工程管理。

课程设置;

基础课：科学社会主义理论、自然辩证法、外语、工程数学、应用弹塑性力学、计算机应用基础、技术经济学等。

技术基础课：高等流体力学、高等渗流力学、油藏数值模拟、油田化学、提高收率原理、渗流物理、油气藏经营管理、运筹学等。

专业课：现代油气井工程、现代完井工程、天然气工程、高等油藏工程、高等油工程、高等输油管道工程、高等输气管道工程、油气田集输系统、油气管道运行模拟、项目管理、能源经济学等。

上述课程可定为学位课程和非学位课程。此外，还可以由培养单位与合作企业根据实际需要确定其他课程。课程学习总学分不少于28学分。

学位论文方面，论文选题应直接来源于生产实际或者具有明确的生产背景和应用价值，或者是一个完整的工程技术项目的设计或研究课题，或者是技术攻关、技术改造专题，或者是新工艺、新设备、新材料、新产品的研制与开发，也可以是工程管理课题。选题要求有难度、有新意、有足够的工作量。

对于技术攻关的成果，应有与国内外同类理论、方法与技术的对析；对于新工具、新工艺设计与开发的技术成果，论文应具有设计方案的比较、评估、参数计算模型与结果、完整的图纸；对于重大工程项目管理的成果，必须给出项目的系统组成、目标分析、风险与效益分析、与管理方案及措施、收益与创新管理方法。

不知道上面这些介绍对这位同学有没有帮助。别人怎么说都只是了解，只有自己才是自己最好的老师！既然这位同学对石油工程这么有兴趣，那么我在这里祝你能心想事成！

自20世纪80年代以来，研究人员提出了许多三维地质模型来模拟地质体，使这方面的研究有了长足的发展。通过对国内外大量的三维地质建模方面的文献和专业软件的研究分析，三维地质建模方法大体可归纳为三类：离散点源法、剖面框架法和多源数据耦合建模法。

1.2.3.1 离散点源法

在地质找矿中，经常需要根据少量的离散点样数据(如地质测绘或钻孔资料)来获取地质体的形状，从而为进一步指导找矿起指导性的作用。因此，研究如何实现空间散乱点数据场可视化的方法具有一定的意义。

Carlson(1987)从地质学的角度提出了地下空间结构的三维概念模型，并提出用单纯复形模型(Simplicial Complex Model)来建立地质模型。Victor(1993)、Pilout(1994)则具体应用Delaunay四面体的三维矢量数据模型研究离散点地质建模问题。Lattuada(1995)对3DDT(3 Dimensional Delaunay Triangulation)在地质领域内的应用进行了研究，表明四面体格网能很好地用于地质体的三维建模，优点包括：四面体单元易于建立索引；模型易于手工编辑；可通过相邻关系导出拓扑结构；约束三角剖分易于实现面约束；四面体非常便于可视化，同时具有较高的表达精度；易于实现搜索和关系查询等。Courrioux et al.(2001)基于Voronoi图实现了地质对象实体的自动重构。Frank et al.(2007)用隐函数法(implicit function)表达三维曲面，对离散点集进行三维重构(reconstruction)，用来模拟断层和盐丘(salt dome)。杨钦(2001，2005)利用离散点源信息构建地层与断层结构面，依此作为约束条件约束Delaunay剖分建立三维地质模型。

钻孔数据也属于一种点源信息。它实质上是将原始的点、线数据进行有效的分层，根据各层面标高应用曲面构造法来生成各个层面或实体。围绕钻孔数据进行三维地质建模已有许多学者进行了研究，其中较早利用钻孔数据进行三维地质模拟的是加拿大学者Houlding(1994，2000)，利用钻孔孔口点位信息进行 Delaunay三角剖分，作为“主 TIN(Primary TIN)”，其他地层面则通过高程映射实现。张煜等(2001)对其建模方法进行了深入研究与发展，在垂直钻孔的理想状态下，用三棱柱(Tri-prism，TP)数据模型建立三维地质模型，并给出了相关的剖切算法。Lemon et al.(2003)用“地层层位法”建立三维地层模型，并用自定义剖面(user-defined cross-sections)的方法对地质模型进行局部交互修正。吴江斌(2003)、朱合华等(2003)提出一种基于钻孔数据的二分拓扑数据结构的建模算法，尝试用基于钻孔数据的四面体体元模型构建地下三维地质模型；四面体结构在表达复杂结构上则较灵活，但是使用四面体表示空间实体会产生大量的冗余，且生成四面体的算法比较复杂。张芳(2005)用Delaunay三角构网技术，利用钻孔数据构建三维地层层面模型，同时引入“界面分片”思想，以适应于海量数据模型的可视化表达，但缺少对地质体属息的表达。在三棱柱模型的基础上，针对钻孔存在偏斜问题，类三棱柱(Analogical Tri-prism，ATP)(齐安文等，2002)、广义三棱柱(Generalized Tri-prism，GTP)(Wu，2004)方法先后被提出，用来进行三维地质建模，已被证明广泛适应于矿山、石油等深部地质问题建模；同时，似三棱柱(Similar Triprism，STP)(Gong et al.，2004)也被提出用于解决钻孔倾斜问题，如郑蔚等(2005)基于钻孔数据用STP建立三维地质模型对地下空间进行虚拟漫游。STP与GTP本质上是相同的。基于钻孔数据建立三维地质建模，这一看似简单的数据模型方法，经历了10多年的发展历程：从初期的TP数据模型，适用于钻孔垂直成层、地层等厚的理想情况，发展到STP、GTP适用于钻孔不垂直且地层不等厚的常见情况。

1.2.3.2 剖面框架法

剖面框架法就是在收集整理原始地质勘探资料的基础上，建立分类数据库，人工交互生成大量的二维地质剖面，然后应用曲面构造法生成各层位面表达三维地质模型，或者利用体元表示法直接进行地质体建模(Chae et al.，1999)。

利用地质剖面表达研究区域三维地质现象的初级形式是序列地质剖面法(朱小弟等，2001)。序列地质剖面构模技术实质是传统地质制图方法的计算机实现，即通过平面图或剖面图来描述地质构造，记录地质信息，如图1.2所示。其特点是将3D问题2D化，在空间上用若干平行或近似平行的地质剖面来表达研究区域的地质分布特征，但它在空间表达上是不完整的，它把剖面之间的地层或构造分布情况留给工程设计人员去“想象”。这种构模方法难以完整表达3D矿床及其内部构造。

基于剖面信息建立真三维模型具有很大的发展空间，对于复杂地质构造区域具有很好的适应性，成为当前地质建模的主要方法之一。然而，基于剖面进行三维重构得到完善发展的是在医学领域，后来迅速扩展到其他领域。在医学领域里，通过电脑断层扫描(CAT)或者核磁共振(MRI)等技术，可以获得一系列相互平行的人体切片图像，通过提取对象的边界，基于轮廓线算法，生成三维人体模型。地质剖面信息同医学切片信息一样，都是反映研究对象的某一特定断面上的构造分布，可以借助医学三维人体建模技术来构造三维地质模型。较早将医学领域的切面三维建模引入地学领域的是在考古学方面的应用(Tipper，16，17；Herbert et al.，1995)，主要应用在古生物的三维重构方面，而应用在三维地质建模方面的文献并不很多。

图1.2 序列地质剖面构模实例

公认的剖面三维重构的代表之作是Keppel的文章(Meyres et al.，1992；Herbert et al.，1995，2001；Xu et al.，2003；屈红刚等，2003)。在Keppel的研究基础上，Meyres(1992)将剖面建模方法分为4个子问题：对应问题(correspondence problem)、构网问题(tiling problem)、分支问题(branching problem)和光滑问题(fitting problem)：①对应问题解决相邻剖面之间的轮廓线匹配问题；②构网问题主要解决轮廓线之间的三角形构网问题，考虑满足某个准则，例如最大体积法(Keppel，15)、最小面积法(Fuchsetal.，17)等；③分支问题是解决同一对象在不同剖面上的组成部分的个数不同的问题；④光滑问题主要解决将初始生成的三角网进行插值，从而得到更加光滑的三角网。

屈红刚等(2003)提出基于含拓扑剖面地质建模方法来实现复杂地质的三维建模的对应问题，邓飞等(2007)则对一般意义上的剖面地质建模进行了讨论。

1.2.3.3 多源数据耦合建模法

随着计算机性能的提高，具备了对海量数据的处理能力，人们对建立的地质模型要求也不断提高，希望能够建立高精度和高复杂度的地质模型(Turner，2003，2006；Calcagno et al.，2006；Kaufmann et al.，2008)。提高模型的精度可以通过插值的方法来实现，但更好的方法是通过增加约束信息来对初始地质模型进行细化，这就涉及耦合多源数据来建立地质模型的问题。

早在1993年，Houlding提出三维地学建模概念的时候就强调地质解释信息具备对模型的修正(revision)功能。并且指出矿业工程有地质勘探数据、人工绘制数据及施工数据，还有不确定性的需要通过地质统计学进行估计的数据(Houlding，2000)，最终的地质模型需要综合考虑这些种类不同的数据。

McInerney et al.(2005a，b)认为三维地质建模只能部分上是一个数字地质样过程，更重要的是地质学家的人工解释过程。并且尖锐地指出，不要指望一些计算机软件能够自动并成功地“建模”! 让一个有经验的地质学家输入解释性的信息进行建模，是现实和必要的；而软件只是建模过程中提供便利的一个工具(There is no expectation that some computer software will successfully and automatically“builda model”! The reality is that interpretative input from a skilled geologist is essential to build a model；the software is simply a tool to facilitate the model-building process)。其要求实际上是，地质建模不仅要考虑地质勘探所获取的确定性数据，还应加入地质工程人员对地质构造的解释性数据，这就构成多源地质建模的基本思想。

Mallet(2002)针对地质体建模的特殊性和复杂性，以点、线数据为主要数据源，建立以三角形为基本单元的三维曲面，用离散光滑插值技术(Discrete Smooth Interploate，DSI)使曲面的粗糙度最小，并作为GOCAD的核心技术，得到了许多地球物理公司和石油公司的支持。

相比较国外以石油、矿业工程为主要应用领域的三维地质建模，钟登华等(2006)则从水利水电工程地质领域，研究多源地质数据建立坝区的三维地质模型。Wu et al.(2005)提出一种逐步细分的多源数据集成地质建模方法，考虑到地质数据大多比较稀疏和低样率的特征，用逐步细化的方法对初始地质模型不断修正。

地质构造的复杂性和认识的阶段性，使多源地质建模引起越来越多的研究兴趣。32届国际地质大会(International Geological Conference，IGC)于2004年在意大利佛罗伦萨召开，在“地质的复兴(The Renaissance of Geology)”(Zanchi et al.，2007)议题上，多名国际知名的地学建模专家共同提到了多源地质建模问题。其中，Zanchi et al.(2008)借助商业软件对意大利境内阿尔卑斯山(Alps)利用多源地质建模问题进行研究，并应用于滑坡稳定性分析。西方发达国家主要将地质建模应用于能源与环境领域，这是为数不多的在工程建设领域开辟蹊径的研究。无独有偶，Kaufmann et al.(2008)尝试用多源地质建模，研究在废弃煤矿巷道内进行天然气储存问题。

总体来看，三维地质建模技术是一个从简单地层模拟到复杂地质构造模拟的发展过程。从最初基于单一数据建立简单层状三维地质模型，到综合利用多源数据建立复杂地质模型，能够反映地质构造的空间特征。