为什么“海洋地质十号”下水填补海洋地质调查空白？

地质局下属的广州海洋地质调查局倒是以海洋石油地质调查为主的事业单位，此外青岛海洋地质研究所也是隶属于地质调查局的海洋地质调查单位。

6月28日上午，我国自主设计、建造的“海洋地质十号”综合地质调查船在广东东莞出坞下水。

广州海洋地质调查局 广州海洋地质调查局统一工作平台

广州海洋地质调查局 广州海洋地质调查局统一工作平台

摘要本文通过海洋成果地质资料涉密清理，总结了成果资料涉密清理经验，并对以后的海洋类原始资料等清理工作提出建议。

地质调查局钟自然在当天的下水仪式上说，该船填补了我国小吨位大钻深海洋地质钻探船的空白，提升了海洋地质调查能力，也标志着我国海洋地质综合调查能力跻身世界前列。

“海洋地质十号”是集海洋地质、地球物理、水文环境等多功能调查手段为一体的综合地质调查船。船身总长75.8米、宽15.4米、深7.6米，结构吃水5.2米，排水量约3400吨，续航力8000海里，定员58人。调查船采用电力推进全回转舵浆，二级动力定位等世界先进航行及控制系统，可以实现在全球无限航区开展海洋地质调查工作，配置了我国首套自主研制的举升式海洋钻探系统。

广州海洋地质调查局叶建良如何说？“海洋地质十号”下水后将进入舾装和设备海试阶段，预计在今年年底与“海洋地质八号”“海洋地质九号”两艘新建调查船一起入列，届时“三兄弟”将与外语：英语一般粤语水平：一般其它外语能力：国语水平：一般现有的“海洋六号”等调查船一起，构成覆盖浅水、深水及大洋等全海域海洋地质调查船舶及技术装备探测体系，海洋地质调查能力将有长足发展。

基于遗传算法的叠前波形反演构建虚拟井曲线

模型为： 。

罗文造1杨绍国2王英民3

F值的计算公式和结果为：

（1.广州海洋地质调查局 广州 510760；2.成都理工大学 成都 610059；3.石油大学（） 102200）

作者：罗文造，男，1969年出生，高级工程师，主要从事地球物理技术方法研究。

摘要 叠前波形反演能够提供详细的地下地层特征，但由于其计算量大、数据和模型之间的非线性、目标函数的多极值和反演结果的多解性使叠前反演的实施成为一大难点。本文通过叠前和叠后混合反演技术采用遗传算法实现了深海无井反演。遗传算法思想简单、易于实现和使用、具有隐含并行性和全局搜索能力等优点，基于遗传算法的叠前波形反演得到了与井中实际数据基本吻合的速度和密度数据。

叠前反演 遗传算法 虚拟井

反演分为叠前反演和叠后反演。一般叠前反演方法所得到的物理特性与叠后反演方法相比提供了更详细的地下地层特征。叠前反演分辨率高，但速度慢、稳定性，现处于研究阶段，距大规模生产应用还有一定距离。叠后反演虽然分辨率要低一些，但速度快、稳定性好，可满足大规模生产应用的需要。

叠前波形反演所面临的难题在于：①计算量和数据量非常庞大；②数据和模型之间高度非线性；③目标函数具有多个极小值；④反演结果具有多解性，可能存在多个模型与数据匹配良好（Sen等，19）。非线性、非性和大计算量交织在一起，使叠前波形反演的难度很大。但它对储层岩性和所含流体的高分辨率，对油气勘探开发技术研究人员来说具有较大的吸引力。针对叠前波形反演所面临的难题，近十年来，有许多地球物理工作者进行了大量的探索，取得了重要的研究成果，这些进展基本上解决了叠前波形反演所面临的高度非线性和局部极小值问题；但对非性和大计算量的问题没得到很好解决。为此，有些学者采用了一种折中的办法，即叠前和叠后混合反演的办法，首先在一些控制点进行精细的叠前波形反演构建虚拟井曲线，然后以虚拟井作为控制信息进行叠后反演（Mallick，2000）。叠前和叠后混合反演的办法利用了叠前反演分辨率高，叠后反演速度快、稳定性好的优点，克服了各自的缺点，成为目前的一个研究亮点。特别对于深海无井的情况，具有很好的应用前景。本文采用叠前和叠后混合反演的办法实现深海无井反演，其关键是叠前波形反演构建虚拟井曲线方法的研究。

2 遗传算法基本原理

遗传算法（Genetic Algorithm，简称GA）是美国Michigan大学的John H.Holland在60年代提出的，目的是把自然界的自适应机理应用于计算机系统的设计。目前，遗传算法作为一种有效的全局寻优自适应概率算法，由于其算法思想简单、易于实现、易于使用、具有隐含并行性和全局搜索能力等优点，且对很多优化问题能够较容易地得到令人满意的解，在自适应控制、组合优化、模式识别、机器学习、人工智能、地球物理反演、管理决策等涉及优化计算问题的领域得到了广泛的应用，且影响越来越大。

尽管遗传算法在解决复杂非线性优化问题中具有独特的优势，但它本身也有局限性，其突出的一个弱点就是收敛性能，尤其对于多参数、超大解空间的优化问题，其收敛速度有时让人很难接受，这在一定程度上制约了GA的使用和发展（Sen等，1992）。采用GA求解高维、多约束、多目标的优化问题仍是一个没有很好解决的课题，它的进展将会推动GA在许多工程领域的应用。

遗传算法是一种基于自然选择和基因遗传学原理的随机搜索算法，但又不同于一般随机搜索算法。它是通过将待寻优的模型空间的参数进行编码，并用随机选择作为工具来搜索过程向着更高效的方向发展，其计算简单、功能强大（Mallick，1995）。具有以下特点：

1）GA是对要寻优参数的编码进行作，而不是对参数本身；

2）GA是从“群体”（多个初始点）出发开始的并行作，可以有效地防止搜索过程收敛于局部是优解，是一种全局寻优的方法，且计算速度较快；

3）GA采用目标函数来确定基因遗传的概率，对问题的依赖性较小，适用范围更广；

4）GA的作均采用随机概率的方式，减少人为干预2 回归分析概述对结果的影响；

5）GA以随机选择来搜索过程，是一种启发式的搜索方法，搜索效率更高；

遗传算法运算的原理简单，只涉及参数编码的和部分编码的交换和变异作，标准GA只包括选择、交叉和变异三种基本作：

（1）选择（Selection）/再生（Reproduction）

选择是根据群体成员的适应度值fi，对群体成员进行的过程。其按照一定的概率选择成员并保留下来，体现“适者生存”的自然规律。

选择运算有许多方法，简单的方法是采用法。该方法根据群体中各成员的适应度，计算适应概率fpi，适应度越大，适应概率越大。

（2）基因交换（Crossover）/重组（Recombination）

基因交换或重组由两步构成，一是匹配，即对再生的群体成员做随机匹配；二是交换匹配成员的基因，即对每对匹配成员，按照随机概率Pc选择随机位置，将两成员的编码进行交换或重组。经过再生所得的新群体中并没有新成员，也就是在搜索空间中没有得到新的搜索点。而按随机概率进行基因交换后，新群体中既有上一代的成员，又有由成员交叉后产生的新成员，得到的新搜索点。根据生物遗传的“杂交优势”规律，这种新成员应该优于原成员。

（3）基因变异（mutation）

基因变异的目的是在搜索过程中，不断引入新的信息量，以免再生和交换的遗传过程中丢失潜在的有用遗传物质。基因变异实际上是对群体成员的基因按小概率Pm扰动而使其发生变化，以补充新的信息。

遗传算法求解优化问题的基本思路：

（1）确定目标函数

目标函数（Objects）是刻画解的标准，也是适度计算的依据。一般是以计算值与观测值之间的拟合程度或误大小为标准的。

（2）指定问题参数的搜索范围

给定每一个模型参数的取值范围，对任一参数x给出［Xmin，Xmax）。实际上是问题解空间搜索精度的描述，其决定了搜索空间的大小，精度越高，离散化的搜索空间越大。

（3）模型参数编码

由于GA运算是对模型参数的编码进行的。编码方式有多种，常用的有二进制编码和十进制编码。

（4）初始群体的产生

给定群体成员数n，用随机生成的方法获取n个群体成员构成初始群体；

（5）遗传计算

① 求取适应度值

适应度值是直接由目标函数转换计算出来的。对小化问题可用

指数转换：

i=1，…，n

其中：Objects［i］为第i个成员的目标函数值；σ为群体目标函数值的方；Fitness［i］为第i个成员的适应度值。

② 再生

Ps［0］=0；

i=1，…，n

生成［0，1］随机数r，如果Ps［i-1］≤r＜Ps［i］则第i个成员获得再生。

将群体成员两两配对，组成n/2对母本，并按概率Pc进行交叉；即产生［0，1］随机数r，如果r＜=Pc，则在一随机位置交换两个母本的编码；否则不变。

④ 变异

对每一个新成员，按概率Pm进行基因的随机突变，即产生［0，1］随机数r，如果r＜=Pm，则在一随机位置改变该成员的编码，由1变0，或由0变1。

遗传算法对问题的求解是通过对寻优参数空间进行编码，从多点出发，采用随机选择作为工具来搜索过程向着更高效的方向发展，是一种普适性的搜索方法（Mallick，1999）。由于在搜索过程中使用了其父辈的适应度函数值作为启发知识，因此又是一种启发式搜索方法。这种搜索方法，对于简单的多极值优化问题可能会产生好的效果。尽管如此，GA对于多参数复杂非线性问题收敛速度还是过慢（Xia等，1998）。

就遗传算法本身而言，提高其收敛速度的关键在于合理的适合问题特点的遗传编码方法、适应度函数及变换方法、遗传算子、算法参数的设置和选取。

3 基于遗传算法的叠前波形反演

采用GA叠前波形反演估算弹性参数。叠前反演在角道集上进行，以便减少计算工作量。叠前资料预处理包括：角道集抽取、叠前去噪、压制多次波和高精度速度分析等方。GA叠前波形反演的技术路线如图1。

图1 GA叠前波形反演框图

Fig.1 The flow chart of pre-stack seiic weform inversion by means of GA

其算法如下：

1）准备初始模型、记录Seis等数据

由高精度速度分析构建初始模型，记录Seis为角道集记录。

2）确定地质模型参数及参数搜索范围和搜索间隔

3）对模型参数编码。根据搜索范围的搜索间隔，先确定各参数可能取得的不同值的个数，为节省空间对所有参数进行整数编码。

4）生成拟合模型的初始随机总体P，设生成了n个随机模型

由X=Xmin+Code×Dx，对Vp，Vs和ρ三个参数nt个样点用随机生成的方式生成整数码产生要求的样本量。

5）计算各模型的合成记录Syni。

6）比较Syni与Seis，计算并保存目标函数值

i=1，…，n

其中：n为群体样本数；从理论上说，点越多，搜索效率应该越高。但实际上增加搜索点，也高增加了遗传计算的计算量。因此解决实际问题时，根据问题的性质及解空间的大小，做适当选择。在计算时，由于遗传计算量相对较大，选择了较小的群体。为便于作和增加程序的适应能力，采用人机交互输入的方式选择8到32间的偶整数。

Nt=nt×angles；

nt：道时间取样点数；

angles：角道集所选角度个数。

seis［j］：观测记录角道集

Syn［j］［i］：第i个群体成员的合成记录角道集

7）根据目标函数值对P做再生、交叉、变异作，更新P生成新的随机总体；

8）如果满足结束条件，结束并输出结果；否则重复5）至8）直到结束。

4 叠前波形反演的实施

针对波反演问题这种多参数、复杂非线性问题的特点采用了如下技术措施：

（1）编码方法

对多参数、复杂非线性问题，其编码的优劣直接影响计算效率。为此采用了整数编码方案有效降低码的长度，加快计算速度。参数编码采用整数编码，方法如下：

参数值=参数小值+码值×参数搜索精度；

码值=（参数值-参数小值）/参数搜索精度；

即：Dx=（Xmax-Xmin）/Codemax

X=Xmin+Code×Dx

Code=（X-Xmin）/Dx

（2）适应度函数

适应度函数是由目标函数转换而得的用以刻划个体适应生存能力的函数。对极小值问题一般采用指数变换，但这种变换是一种均匀变换，在计算后期当群体中各样本目标函数值接近时，为增加个体在再生时被选中的可能性，从而加快算法收敛，选择采用了S函数做叠加变换。

开始时使用指数变换：

i=1，…，n

其中：Objects［i］为第i个成员的目标函数值；σ为群体目标函数值的方；Fitness［i］为第i个成员的适应度值。

当群体中样本目标函数值接近时使用S函数变换：

用于在遗传迭代计算后期，当群体各样本适应度很接近时，以指数形式放大平均适应度以上的样本适应度异，缩小平均适应度以下的样本适应度的异，以便更好地选择个体（Sen，2001）。

当然为达到上述目的，也可使用其它函数，选择使用S函数，一是因为S变换的连续性，可使大于平均值的适应度放大，而使小于平均值的适应度缩小。二是基于前人的经验增加个体在再生时被选中的可能性，从而加快算法收敛。

（3）算法过程

一般GA在计算时采用的是上一代的适应度作为启发函数再生后进行的随机启发搜索方法。为提高算法速度在实际处理中除使用上一代的适应度，还充分利用了的隔代遗传的信息作为启发信息，参与遗传过程的计算。采用一种有限深度回溯搜索的方法，避免了迭代计算的反复，从而加快了计算收敛速度。事实上，在超大解空间中，某一代的遗传性能往往很难决定终结果的好坏。另外，在交叉中每对成员交叉变换使用两次概率选择方法，即先选成员对，再选参数，且每个参数分别选择，这样可以有效地增加搜索能力。

5 试算效果实例分析

图2给出了理论记录的基于遗传算法的叠前波形反演的实例。图中展示了反演纵波速度和实际井中速度曲线，反演横波速度和实际井中速度曲线，反演密度和实际井中密度曲线，理论角道集记录和反演结果的合成角道集记录。基于遗传算法叠前反演结果与井中实际数据基本吻合，理论角道集记录和反演结果的合成角道集记录吻合非常好。

图2 基于遗传算法的叠前波形反演实例

Fig.2 A case for GA based pre-stack seiic weform inversion

图中从左到右为反演纵波速度和实际井中速度曲线，反演横波速度和实际井中速度曲线，反演密度和实际井中密度曲线，5°理论角道集记录和反演结果的合成角道集记录，10°理论角道集记录和反演结果的合成角道集记录，15°理论角道集记录和反演结果的合成角道集记录，20°理论角道集记录和反演结果的合成角道集记录

参考文献

Inger L，and Rosen B.1992.Genetic algorithms and very fast simulated annealing：A comparision，Math.Comput.Modelling，16，87～100

Mallick S.1995.Model-based inversion of amplitude-variation-with-offset data using a genetic algorithm，Geophysics，52，1355～1364

Mallick S.2000.Hybrid seiic inversion：A reconnaissance tool for deepwater exploration.The Leading Edge，19，1230～1237

Sen M K，and Stoffa P L.19.Nonlinear one-dimensional seiic weform inversion using simulated annealing，Geophysics，56，1624～1638

Sen M K.2001.Pre-stack weform inversion：Current status and future direction，Institute for Geophysics

Xia G，Sen M K，and Stoffa P L.1998.1-D elastic weform inversion：A divide-and-conquer approach，Geophysics，63，1670～1684

The construction of pseudo-well logs by inversion of pre-stack seiic weform based on genetic algorithm

Luo Wenzao1Yang Shaoguo2Wang Yingmin3

（1.Guangzhou Marine Geological Survey，Guangzhou，510760；2.Chengdu University Sciences and Technologies，Chengdu，610059；3.China Petroleum University，Beijing，102200）

Abstract：The inversion of pre-stack seiic weform is capable of providing the detailed character of subsuce stratigraphy.But its large scale of calculation，the non-linear relationship between data and model，multi-pole of objective functions and its multi-resolution has resulted in the difficulty of its operation.The seiic inversion without well constrained in the deep sea has been realized by hybrid inversion of pre-and t-inversion.The genetic algorithm is ，easy to realizing and using with the character of latent paralleling and global searching capability.The inversion of pre-stack seiic weform based on genetic algorithm gives the velocity and density which consists basically with those from the wells.

Key Words：Pre-stack inversion Genetic algorithm Pseudo-well

石油地质海洋地质局现在是什么单位?

合成记录采用Zoeppritz方程计算。

不明白你说的“石油地质海洋地质6）GA适合于大规模复杂问题的优化，如反演问题。局”是个什么单位

是中4.1 原始资料归档要完整、及时国地质调查局下属的单位

沉积物粒度分析软件程序设计

Fig.3Showing tilted unconformity of single seiic profile D1A

梁蓓雯1王金莲1朱美琼2

（1.广州海洋地质调查局，广州，510760;2.广东省物料实验检测中心，广州，510080）

作者：梁蓓雯，女，1973年生，广东人，计算机应用工程师，1996年毕业于成都理工大学计算机科学与应用专业，现从事计算机程序开发研制工作。

摘要本文介绍了计算机处理沉积物粒度分析数据的程序设计过程和原理，生产和科研中运用该程序，解决了过去人工绘图及计算器处理数据误大，容易出错的问题，使粒度分析数据处理过程效率大为提高，处理结果性提高。

沉积物粒度分析程序设计

1前言

沉积物粒度分析是地质实验中非常重要的测试项目，它为沉积岩石学研究——判别沉积环境提供重要的相关参数（孙永传等，1986）。多年来，该测试项目的大量数据都依靠人工计算与人工绘图，整个数据处理过程耗时长，人为因素控制数据结果影响很大，并且计算很容易出错，结果误也较大，工作效率很低。在这种情况下，为改善这一状况，研制这一测试项目数据处理的先进计算机程序非常必要，通过自动化的数据处理，提高工作效率及结果精度。

本文以计算机程序处理的沉积物粒度分析结果数据为基础，详细介绍了程序设计过程及设计原观测记录与合成记录之间的匹配程度称为模型的拟合度（Inger等，1992），如果随机模型与实际情况相很远，由观测记录计算得到的角道集与相应的合成角道集匹配就会。相反如果所选随机模型接近实际情况，从而使由观测记录计算得到的角道集与相应的合成角道集能很好地匹配。理，以及重要的计算过程和计算原理。

2作平台和编程语言的选择

现在使用的微机多为PⅢ档次以上，作系统多为windows me/XP，为了扩大应用范围，选择了windows me为本程序的编程平台，使本程序既可以在windows98上运行，又可以在windows XP作系统上运行。编程语言种类繁多，应根据具体情况进行选择。Visual Basic是一个集应用开发、测试、查错功能于一体的集成式开发环境，对Active X的支持，使Visual Basic开辟了新的编程领域，Visual Basic以其面向对象的程序设计、灵活的控件运用、良好的编程界面、方便的数据库编程等功能而非常流行（李天启，1999），因此选择Visual Basic6.0作为本程序的编程语言。

3计算模式的选择

沉积物粒度分析结果综合整理方法有计算法即矩值法和作图法两种（成都地质学院陕北队，1976）。作图法是求粒度频率分布的较好方法，常用的粒度图有直方图、频率曲线图和累积曲线图。概率累积曲线图的特征是能够以直线段明显地反映和区分各对数正态总体，分析结果很直观，可以直接从概率累积曲线图上看出沉积岩的沉积相变，并把反映沉积能量条件为敏感的部分——粗、细尾部大大地扩大（孙永传等，1986），十分有利于环境分析，而且作图简便，一个样品即可作出一张概率累积曲线图，这是矩值法无法比拟的。现行标准及很多沉积学家采用了概率累积曲线图进行分析结果的计算，因此，研制以概率累积曲线图为计算模式的计算机程序，继承了粒度分析测试项目的传统和特色，并拓展了该测试项目的程序设计领域。

4程序构成

沉积物粒度分析数据处理程序的构成框图如图1，整个程序分解成上述几个模块，在程序设计中既互相又互相联系。

5功能模块设计

5.1分析数据输入模块

分析数据的输入利用Text Box控件功能，对关键的Text Box控件进行一定的容错设计，利用帮助模块指导用户输入数据，尽量方便作和减少错误，见图2。输入的数据文件保存在一定的目录下，方便用户调用和修改。

根据输入的分析数据，计算沉积物各粒级的累积百分率。因为分析仪器对粒度的筛选存在一定误，本模块具有一定的容错性，允许样品的累积百分率大于或小于。

5.3画正态分布图

式中：σ——标准偏；α——平均值。

当σ＝1、α＝0时，为标准正态分布，在程序设计中，要计算概率0.01～99.99的标准正态分布的分位数，因此，使用了如下正态分布的分位数计算函数：

南海地质研究.2003 图1程序流程图

Fig.1Flow chart of program

图2分析数据输入窗

Fig.2Input intece of data ysis

b0=0.1570796288×10b6=－0.1045274970×10-5

b1=0.3706987906×10-1b7=0.8360937017×10-7

b2=-0.8364353589×10-3b8=-0.3231081277×10-8

b3=-0.2947176×10-3b9=0.3657763036×10-10

b4=0.6841218299×10-5b10=0.6936233982×10-12

b5=0.5824238515×10-5

通过利用上述函数，使用循环设计，概率从0.01递增到99.99，计算出每个概率级（0.01）的分位数：

（保存K值到数据库以备后用）

Next Q

计算好Kβ值，把Kβ换算成对应计算机屏幕上的具体坐标值，把这些一一对应的数据制成一个正态分布概率数据表，供程序画图调用，数据准备好后，使用Picture Boxes控件，把该控件的Auto Redraw属性设置为True，再使用Visual Basic的Line、Circle等作图功能画出正态分布概率图，并利用Scale Height、Scale width属性，方便图形设计。

5.4图形显示模块

利用上一模块中得到的数据，在正态分布概率图上，画出样品各粒级的累积百分率所对应的点，同时显示连线的点号，让用户选择。连线的设计是本程序的重点，也是难点之一，关系到各粒度参数的计算结果。本模块设计了三种连线方式。①自动连线方式：选择此方式，程序会根据一定的误参数值，在屏幕上随机生成M条直线，用户可以用此方式对曲线图先作一个初步的绘制。②抛点自动连线方式：用户可以选择此方式抛弃掉1～2个偏较大的点，然后由程序自动绘出曲线图。③手动绘图方式：用户可根据点的走势，输入点的编号，确定N点和M条线，生成曲线图。见图3。本模块运用Visual Basic的Picture Boxes控件和Line、Point、Scale等绘图功能。

5.5计算各类参数

根据所画的曲线图，得到每条直线的方程，并根据要求的特定的百分比值所对应的Y值，代入所对应直线方程，算出X值即Φ值，然后再按以下公式计算各种粒度参数Mx值，σi值，Ski值，Kg值。

直线方程公式：

式中：（x1,y1），（x2,y2）是直线通过的两点。

粒度参数计算公式：

图3编号1305样品的正态分布概率曲线图

Fig.3Normal probability curve of No.1305sample

式中Φn表示累积曲线上n%对应的Φ值。

5.6打印图表、报告

[object.]Print[{Spc（n）｜Tab（n）｝][expressionlist][{;I,}]

报告为一个包含各种数据的表格和一个正态分布概率曲线图。

5.7输出图表、数据

所得数据和图形与WINDOWS应用软件（如WORD，EXCEL）动态连接，便于地质解释时数据进一步与应用软件（如GRAPHER，SURFER）连接。

Mallick S.1999.Some practical aspects of prestack weform inversion using a genetic algorthm：An example from the east Texas Woodbine gas sand，Geophysics，64，326～3366总结

通过运用本程序，减少了粒度分析数据处理的时间，使整个处理过程得到简化，提高了劳动效率，并提高了样品分析结果的质量。另外，使粒度分析资料结果可以以电子形式得以保存，便于资料的保存和随时查询，为以后大量利用数据资料提供更大的方便。

参考文献

成都地质学院陕北队编.1976.粒度分析在区分沉积环境上的应用.：地质出版社

李天启编著.1999.Visual Bisic 6.0学习捷径.：清华大学出版社

秦建侯，邓勃，王小琴.1990.分析测试数据统计处理中计算机的应用.：化学工业出版社

孙永传，李蕙生编.1986.碎屑岩沉积相和沉积环境.：地质出版社

Programming of the Grain Size Analysis of the Sediment

Liang Beiwen Wang Jinlian Zhu Meiqiong

（Guangzhou Marine Geological Survey,Guang zhou,510760）

Abstract:The introduced the comr programming process and the principle of the grain size ysis of the sediment.Using the program,we can solve many problems in pract,such as great error,fallibility in plotting and calculating.The program will improve the efficiency and precision of the results.

Key Words:SedimentGrain size ysisProgramming

“海洋六号”完成南极科考任务是怎么回事？

使用Visual Basic的Printer对象来把处理结果报告输送到打印机（李天启，1999），其用法为：

海洋六号”综合调查船是由船舶及海洋工程设计研究院自行设计、建造的综合地质地球物理调查船，集2.4.3 变量显著性检验（t检验）：、地质调查等多项调查功能于一体，总造价近4亿元。

2008年10月在武昌造船厂建成下水，该船为广州海洋地质调查局建造，以海底“可燃冰”5.2计算各粒级的累积百分率模块调查为主，能在海上航行60天无需补给。2018年11月11日，“海洋六号”科考船完成深海地质调查第6航次和大洋第51航次科考任务，返回广东海洋地质专用码头

南海东北部晚第四纪地层不整合的发现及其地质意义

读取回归结果如下：

王树民陈泓君钟和贤

③ 交叉

摘要根据D1A段单道剖面反射特征发现了南海东北部晚第四纪地层不整合，并被第四系浅钻ZK3孔岩心的岩性、古生物和测年所证实。上部不整合造成中更新统顶部至上更新统上部间的缺失。该不整合主要是由于末次冰期盛冰期海平面下降形成，其次是受末次冰期前局部构造变动的影响。据测年推算.15.6ka B.P.海平面相对现代海平面下降了146m，当时的古海岸线位于距今海岸约170km处。这就为南海东北部末次冰期盛冰期海平面下降幅度提供了实际证据。

南海东北部晚第四纪地层不整合末次冰期盛冰期海平面变化

本文资料主要来自1999年12月完成的“F50（汕头）1:100万海洋区域地质调查与编图”项目中海洋地质补充调查成果。补充调查在该区除了按规范进行海底地质取样外，还完成了四个第四系浅钻和两条模拟单道剖面测量（图1）。在D1单道剖面的D1A段和ZK3浅钻中发现了中更新世晚期与晚更新世早、中期地层不整合。这个不整合对南海东北部晚第四纪地质研究具较大意义。本文是在对D1A段单道剖面的解释和对ZK3孔岩心样品测试的基础上，进一步论述这个不整合的特征、形成时间及其形成原因。

图1D1A单道剖面和ZK3钴孔位置图

Fig.1Location of single1 前言 seiic profile D1A and core ZK3

1不整合的特征

这个不整合不仅在单道剖面中有非常明显的反映，而且在ZK3孔岩心的岩性、古生物、测年中得到证实。

1.1不整合在反射界面中的反映

D1A段单道剖面位于浅滩以西外陆架至上陆坡部位，水深为45～375m。在获得清晰反射记录中，自上而下可划分出6个主要反射界面（R0～R5）和相应的层组（A～F）.在R2～R3和R3～R4间，还可以划分出 、 、 三个次要界面（图2）。在41点、46点和53点的反射界面中，很明显地发现了R1切割了R2，形成了本区规模较大的地层不整合（上部不整合）；在48～51点间，又明显地发现R2切割了 和 形成下部局部不整合（图2），显示了上、下界面的交角达10°左右（图3）。

1.2不整合在ZK3孔岩心中的岩性、古生物及测年中的反映

1.2.1不整合面为明显的岩性界面

剖面显示的上部不整合被ZK3孔揭露，孔深11.7m为不整合面。该面上覆沉积层为粘土质粉砂，其中粉砂质量分数约54.6%，粘土约33.3%；碎屑中有丰富的贝壳碎片，还有少量的黄铁矿（0.032%）和海绿石（0.005%）。下伏沉积层为粉砂质细砂，其中细砂含量约51.7%，粉砂约29.9%；碎屑中有贝壳碎片和完整贝壳，也有少量黄铁矿（0.022%）和海绿石（0.007%）。下伏沉积层具有明显的微斜层理，与上覆沉积有明显的构造异。可见，这个不整合为明显岩性异界面，而且钻孔揭露的深度与解释界面（R1）较为吻合。

1.2.2不整合面为古生物不连续面

为了满足1:100万编图，对ZK3孔岩心进行较大间距的古生物分析，其中在孔深8.0m和12.5m处做了有孔虫和钙质超微化石鉴定（由广州海洋地质调查局实验测试中心承担）。鉴定结果为：8.0m处的浮游有孔虫含量为42.03%，钙质超微化石为NN21带；12.5m处的浮游有孔虫含量为15.15%，钙质超微化石为NN19带（图4）；据此并结合岩性可推测，上部不整合面（R1）的上覆沉积为水较深的浅海环境；下伏沉积为水较浅的滨海环境。其间钙质超微化石缺失了NN20带，到底是由于不整合导致沉积间断，还是由于取样间距过大造成的？结合岩性特征分析，很可能是不整合造成。可见，不整合面也为古生物不连续面。

1.2.3不整合缺失一段时间的地质记录

为了对本区第四纪地层进行划分，还对ZK3孔岩心用14C和热释光法作年龄测定（由科学院广州地球化学所承担）。其中在孔深4.0m、8.0m、12.5m和18.0m别测得年龄为（14730±330）a、（15.2±1.2）ka、（212.5±17.3）ka和（288.5±23.0）ka。用沉积速率法分别计算了上部不整合面上、下点的年龄值。结果为：不整合面上点年龄为15.6ka，下点年龄为201.0ka。可见，该不整合面以上为晚更新世晚期以来的沉积，以下为中更新世晚期以前的沉积，其间缺失了中更新世晚期顶部至晚更新世早、中期的沉积，即大致缺失了185.4ka的地质记录。

综上所述，由单道发现的南海东北部晚第四纪地层不整合，与ZK3孔岩心岩性、古生物、测年所表现出的特征较为一致。据此，尽管对ZK3岩心取样间距较大，但仍可肯定这个不整合的存在，并能大致确定它的形成时间。

图2D1A单道剖面解释图

Fig.2Single seiic profile D1A interpretation

R0～R5—反射界面； —Qhbj—上更新统南卫组上段至全新统北尖组； —中更新统北卫组上段、中段、下段；Qp1—下更新统；Ⅰ—级陆坡；Ⅱ—第二级陆坡；Ⅲ—第陆坡；Ⅳ—第四级陆坡；50m—ZK3井终井深度；11.7m—不整合并深；－146m—古海岸线距今海面深度

图3D1A单道剖面显示的倾斜不整合

图4ZK3孔综合柱状图

Fig.4Core ZK3 gcneral column

2南海东北部晚第四纪地层不整合的地质意义

众所周知，一般的地层不整合是由于地质构造运动造成，海区第四纪地层不整合主要是由于海平面变化造成。海区第四纪海平面变化与古气候和新构造运动密切相关。下面仅就上述的不整合讨论南海东北部末次冰期盛冰期海平面位置及其形成原因。

2.1南海东北部末次盛冰期古海面高度

某一地质时期的古海面可用沉积、古生物、海底古地貌和古岸线等确定。南海东北部晚更新世晚期的古岸线，据D1A单道剖面（图2）分析，应在上部不整合面（R1）顶超的点（53点）附近，距今岸线约170km。该点现代水深为131m,R0～R1的沉积层厚约15m，即当时古海面位于现代海面—146m；其形成时间，应在据ZK3孔确定的15.6ka以前。关于玉木冰期盛冰期海平面的下降幅度，国内外学者做过较多研究，认为黄海和渤海沿海地带距今20.0～15.0ka的古海平面在－150～－160m左右（赵希涛等，1979），南海北部约在－130m左右（冯志强等，1996）。可见，本区末次冰期盛冰期海平面下降幅度及形成时间与以往推测的较为接近，这就为南海北部玉木盛冰期海面位置提供了实际证据。

2.2南海东北部晚第四纪海平面变化的原因

第四纪全球海平面变化主要与全球性的古气候变化相关，其次还受构造运动和沉积物压实影响。据R.F.Flint（1971）计算，末次冰期盛冰期时，全球冰约71.36×106km3（折合成水），冰后期融化了47.3×106km3，可使海面上升132m，若考虑海底均衡下降1/3，则可使海面上升92m。上述本区末次冰期海平面下降了146m，是相对现代海平面的综合海平面变化，其主要原因，据所推算的形成时间无疑是由于末次冰期盛冰期造成的，但也不可避免地受到构造作用和沉积物压实的影响，若设这些影响也约占1/3，那么实际海平面下降约为97m。据ZK3孔推算，上部不整合缺失了201.0～15.6ka的沉积，这一时期已超出末次冰期下限；而且D1A剖面ZK3孔处B层（R1～R2）厚为2.4m，而37点和42点处B层厚约18m。显然，201.0ka后至末次冰期前是有沉积的，不同地段B层厚度不同是由于断层（F）及局部地壳升降异影响，使下降盘保留了较多的沉积。因此，可以推测本区古海面变化中明显地受到局部构造变动的影响。同样，下部倾斜不整合的形成也可能受到更早期构造影响。

据新构造运动研究，晚第四纪南海北部沿岸以地壳断块异运动和间歇升降为主，而且是西弱东强，据此及D1A剖面所表现的现象，推测本区在末次冰期之前地壳有一次振荡运动，使地层发生轻微褶皱和断裂作用，然后由于末次冰期海平面变化才造成地层被剥蚀而形成不整合。类似这样的过程，在这以前还可能有几次，从图2可以看出，南海东部陆架和陆坡从早更新世以来至少经过四个阶段，由于海平面多次变化和地壳运动，沉积物逐步向海方向进积而形成的（图2、图5）。

3结论

（1）南海东北部晚第四纪以来地层存在两个较明显的不整合，其中以中更新统顶部与上更新统上部不整合明显，且规模较大。这个不整合不仅在单道反射波中有清楚的反映，而且被ZK3孔岩心的岩性、古生物和测年资料所证实。

图5南海东北部晚第四纪地层不整合演化简图

Fig.5A sketch map showing the uncomformable evolution of the Late Quaternary strata in northeastern South China Seas

（2）该不整合主要是由于末次冰期盛冰期海平面下降造成，其次是受末次冰期以前地壳异升降造成地层轻微褶皱和断裂的影响。本区在15.6ka海平面相对现代海平面下降146m，当时的古海岸线位于距今海岸线约170km处。这就为本区末次冰期盛冰期海平面下降幅度提供了实际证据。

参考文献

1.科学院南海海洋研究所海洋地质构造研究室.1988.南海地质构造与陆缘扩张.：科学出版社，307～312.

2.冯文科，薛万俊，杨达源.1988.南海北部晚第四纪地层环境.广州：广东科技出版社，183～1.

3.黄玉昆，邹和平，张珂.1995.新构造学.广州：广东省地图出版社，229～249.

4.冯志强，冯文科，薛万俊等，1996，南海北部地质灾害及海底工程地质条件评价.南京：河海大学出版社，17～49.

5.杨子赓主编.2000，海洋地质学.青岛：青岛出版社，283～301.

THE FOUNDATION OF LA TE QUATERNARY STRATA UNCONFORMITY AND THEIR GEOLOGICAL SIGNIFICANCE IN THE NORTHEASTERN SOUTH CHINA SEAS

Wang Shumin Chen Hongjun Zhong Hexian

（Guangzhou Marine Geological Survey,MLNR,Guangzhou,510760,China）

Abstract

According to the reflection characters of D1A segment in the single seiic profile,we discover the late Quaternary strate unconformity in northeastern South China Seas,which is confirmed by the lithology,palseobiology and age of borehole ZK3.The upper unconformity engenders the sediment absence of top middle Pleistocene to top upper Pleistocene.It formed firstly by the sea ll falling during pleniglacial of kataglacial;secondly,affected by the local tectonic variance before the kataglacial.The sea ll has fallen relatively 146m to modern sea ll in 15.6 ka B.P.,when the paleocoast-line sited off 170km to modern coast-line at that time,and provided the evidence of the falling of sea ll during pleniglacial in northeast area of the South China Seas.

Key words:northeastern South China Seas;late Quaternary;stratal unconformity;kataglacial and pleniglacial;sea ll changes

注释

浅谈海洋地质资料清理的思考和认识

南海地质研究.2003

唐卫寇新琴徐梅影

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

海洋地质资料原始资料涉密清理

通过对保管的成果地质资料进行涉密清理、准确定密，进一步加强资料管理单位对涉密成果地质资料的管理，确保地质资料的安全使用，充分发挥现有地质资料的作用，避免工作重复和资料浪费。

1 海洋地质资料的特点

随着海洋世纪的到来，海洋资源开发、海洋环境保护、海洋权益维护与海洋科学研究等对全方位、高质量的海洋信息服务需求越来越迫切，越来越离不开海洋地质工作。

海洋地质工作与陆地地质工作相比，更为复杂。海洋地质调查一般以项目为主，参加与合作的单位多、周期长、规模大。如1:100万海南岛区调项目自2004年开始已有6个年头，其中调查手段有、重磁、导航、多波束、地质取样、地温梯度测量等，并且不同调查航次开展的调查手段不同，因此海洋地质资料数据具有获取难度大、成本高、种类多、数量大等特点。

五十多年来，广州海洋地质调查局（以下简称“广海局”）在我国黄海、东海、南海、南沙（群岛）以及海域开展了海洋油气、天然气水合物、海洋区域地质调查、近岸海洋地质环境与地质灾害、大洋多金属资源调查等工作，取得大量非常宝贵的海洋地质资料。

2 海洋成果地质资料清理

为了贯彻落实《地质资料管理条例》，更好地发挥海洋地质资料的作用，2007年，广海局资料档案室开展本单位馆藏的成果地质资料涉密清理工作。

2.1 清理的依据

此次海洋类成果地质资料的清理依据主要有：国土资发［2003]147号文、国土资发［2008]69号文、国测办字［2003]17号文、国海密字［1996]450号文等。

过去对成果资料密级界定主要是根据任务书定，没有细分到每一件，密级比较笼统，在此次清理工作过程中，凡是清理范围内的成果资料，不论以往是否定过密级，项目人员一律按照现行的标准，以认真、严谨的态度对每一份成果资料的报告、评审意见、附图、附表、附件等进行逐档、逐件的填表登记，密级和期限进行重新标注。同时如果一份成果地质资料符合几条相关的法律条款，我们本着就高不就低的原则，对成果地质资料进行密级的界定。将采集的所有信息录入到国土资源部统一设定的“地质资料汇交清理登记系统”中，从而建立文件级地质资料涉密清理数据库。

2.3 清理的结果

通过此次清理工作，摸清了家底，对在清理范围之内的多套资料进行了重新分类、编号、数据库录入、涉密定级，对破损的地质资料进行修复、统计，形成了式中：x对应用不同样本的原适应度值；y为变换后的适应度值；θ0为所有样本的平均适应度。a＞0表示用于控制放例参数，越大对平均值以上的部分放大越明显。b＞0表示调节系数，当a=1时，可取b为8到10；b太大达不到对接近值处的适应度的放大，b较小时可用线性变换取代。成果地质资料清理数据库及相关技术文档。

同时依据不同的标准对广海局的资料进行分类的整理，满足不同功能的需求。例如我们按照地调资料与非地调资料分类；按照地域进行分类，如南海资料、东海资料、大洋资料等；按照项目分类如水合物资料，北黄海资料等；按照密级种类的变化等分类。根据这些分类能很快地找到相关的内容，方便以后的资料编研工作。

3 海洋成果地质资料的服务利用

由于海洋地质工作所形成的资料密级程度比较高，为此，我们坚持依法治档，狠抓规章制度建设，对于有密级的成果资料，强化管理，如对已借出的涉密资料进行跟踪管理并与借阅人签订协议，确保秘密和地质资料的安全使用。

3.2 公开资料目录，提供网络查询

成果地质资料涉密清理完毕后，实现了由传统的按档定密管理提供利用模式向按件定密管理提供利用模式的转变，解放出了一大批资料。对能公开的地质资料，全面公开地质资料目录，而对于涉密资料，利用先进的计算机技术将地质资料中涉密要素去除，使其转化为公开类地质资料，提供网络查询。

3.3 加强数字化建设，促进老资料的二次开发

天然气水合物的发现可以说是老资料的二次开发利用的证明。专家学者们主要是通过对以前广州海洋地质调查局在我国南海进行油气调查时形成的剖面上发现了具有水合物特征BSR，从而预言南海可燃冰前景，开启了水合物新能源的探索航程。老资料的实用价值还有待我们进一步挖掘，但由于保存的年代已久，使用频繁，这些资料有的已模糊破损，因此，应充分利用先进的科技力量，加强数字化建设，对老资料进行有效的保护和利用，确保其安全、完整和长期可用，为提供高效、优质的服务。

4 对海洋类原始地质资料清理的建议

此次海洋成果资料清理工作积累了不少经验，为以后的原始资料、实物地质资料清理工作奠定了基础。同此，对海洋类原始地质资料清理提出自己的建议。

一个正常实施的地质项目，归档的原始地质资料必须要齐全成套，从立项开始到审查验收结束，经历的每一个过程，包括直接形成和采集的，以各种载体类型存3.1 加强管理，遵循制度在的原始记录、中间性（过渡性）产品、终地质工作成果等资料均属归档范畴。

4.2 涉密清理的依据要充分

在成果资料清理过程中发现：对于海洋类资料密级界定的依据太过笼统，作性不是很强，应该细化，如果是根据该项目任务书来定密，一个项目一般就对应一个成果报告。但对于种类繁多的原始资料、实物地质资料很容易形成按档定密而无法细分到件。

4.3 “地质资料汇交清理登记系统”要尽早完善

成果资料清理时研发的“地质资料汇交清理登记系统”，对海洋类密级代码设置不是很完善，在进一步设计时可以参考对于海洋类资料密级的划分来设置密级代码，这样能更好更快地完成资料的清理工作。

原始资料汇交清理软件系统能和此次成果“地质资料汇交清理登记系统”共用一个开发平台，这样对档案管理人员查找资料很方便，并且能推进成果资料、原始资料和实物资料一体化。

基于回归分析的海洋地质调查研究及实例应用

2.2 清理的方法

梁广1，2邵长高1，2

计算再生概率：

（1.广州海洋地质调查局 广州 510760；2.国土资源部海底矿产资源重点实验室 广州 510760）

作者：梁广（1972—），男，工程师，主要从事网络管理和数据管理工作，。

摘要 近年来资源勘探已经覆盖大部分陆地区域，越来越多的把目光投向海洋。海洋作为一个巨大的能源和资源宝库在国民经济、军事战略等的重要性也日益显现。各个竞相制定海洋科技开发规划、战略，优先发展海洋新技术［1］。如何有效的从海量海洋地质调查数据中获取有用信息是海洋新技术研究中的重要研究内容。论文针对海洋地质调查数据研究技术应用需求，引入了回归分析模型到海洋地质调查数据库中，详细介绍了回归分析的技术方法和在海洋地质调查数据库研究中的应用优势，为海洋科学研究提供了技术支持。

海洋地质 回归分析 数据库

随着陆地资源的消耗和人类对能源越来越强烈的需求，海洋作为一个尚待大规模开发的能源和资源宝库引起各国越来越多的关注。我国作为世界上的发展家对能源的需求也在大幅增加，近年来我国石油进口数量急剧增长，据估计到2020年我国石油进口依存度将达到60%。和人多次提出“资源、能源、特别是油气资源，已成为我国经济和发展的重要因素，解决后备能源问题是保证经济安全的大事”。随着我国国土资源大调查和海洋地质专项调查的开展，大量的海洋地质数据被收集和积累，并建立了多个满足各自业务需求的信息系统和数据源［2］。如何有效的从海量海洋地质调查数据中获取有用信息是海洋新技术研究中的重要研究内容。论文针对海洋地质调查数据研究技术应用手段的需求，引入了回归分析技术到海洋地质调查数据库中，详细介绍了回归分析的技术方法和在海洋地质调查数据库研究中的应用优势，为海洋科学研究提供了技术支持。

2.1 概述

回归分析是确定两种或两种以上变数间相互依赖的定量关系的一种统计分析方法。回归分析按照涉及的自变量的多少，可分为一元回归分析和多元回归分析；按照自变量和因变量之间的关系类型，可分为线性回归分析和非线性回归分析。如果在回归分析中，只包括一个自变量和一个因变量，且二者的关系可用一条直线近似表示，这种回归分析称为一元线性回归分析。如果回归分析中包括两个或两个以上的自变量，且因变量和自变量之间是线性关系，则称为多元线性回归分粒度分析中运用的图为正态分布概率图，该图的设计是本程序的重点之一，也是难点之一。正态分布为一个定积分。（秦建侯等，1990）：析［3］。回归分析预测法可以从各数据之间的相互关系出发，通过对与预测对象有联系的现象变动趋势的分析，推算预测对象未来状态数量表现的一种预测方法，通过对与预测对象（y）有联系的多个因素X1，X2，……，Xk建立回归模型。求出的回归模型是否合理，是否符合变量之间的客观规律性，引入相关因素是否有效，变量之间是否存在线性相关关系，模型能否付诸应用，这要通过检验决定。本文给出了两方面的检验：一方面为实际意义检验。即利用理论所拟定的期望值与实际结果相比较是否相符。另一方面为统计检验：分别为拟合优度检验（R平方检验）、方程显著性检验（F检验）、变量显著性检验（t检验）［4］。论文主要介绍一元线性回归分析在海洋地质调查数据库中的应用。

2.2 一元线形回归分析模型

线性回归分析可以描述两个要素之间的回归关系。线性回归分析公式为：yi=a+bxi+εi.其中a和b为参数.εi是误.我们定义Q（a，b）a为总误。则：

南海地质研究（2014）

对公式两边的a和b求导得：

南海地质研究（2014） 南海地质研究（2014）

x表示x的平均值.y表示y的平均值.

关系系数R2求值方法为［5］：

南海地质研究（2014）

2.3 多元线形回归分析模型

研究对象y受多个因素x1，x2，x3，…xn的影响，定各个影响因素与y的关系是线性的，则可建立多元线性回归模型：

y=β0+β1x1+β2x2+…+βkxk+ε

式中：x1，x2，……，xk代表影响因子；ε 为随机误；y 代表所研究的对象，即预测目标［3］。

2.4 统计检验

统计检验是运用数理统计的方法，对方程进行检验、对模型参数估计值的可靠性进行检验。这主要包括拟合优度检验、方程显著性检验、变量显著性检验，即常用的R2检验、F检验和t检验。

2.4.1 拟合优度检验（检验）：

拟合优度检验就是检验回归方程对样本观测值的拟合程度。又称为复相关系数检验法，它是通过对总变（总离）的分解得到。

南海地质研究（2014）

其中

南海地质研究（2014）

总变平方和S总是各个观察值与样本均值之的平方和，反映了全部数据之间的异；残平方和S残是总变平方和中未被回归方程解释的部分，由解释变量x1，x2……，xk中未包含的一切因素对被解释变量y的影响而造成的；回归平方和S回是总变平方和中由回归方程解释的部分。对于一个好的回归模型，它应该较好地拟合样本观测值，S总中S残越小越好。于是可以用：

南海地质研究（2014）

2.4.2 方程显著性检验（F 检验）：

对于多元线性回归方程，方程显著性检验就是对总体的线性关系是否显著成立作出推断，即检验被解释变量y与所有解释变量X1，X2，……，Xk之间的线性关系是否显著，

南海地质研究（2014）

即F统计量服从以（k，n-k-1）为自由度的F分布。首先根据样本观测值及回归值计算出统计量F，于是在给定的显著性水平a下，若F＞Fa（k，n-k-1），则拒绝H0，判定被解释变量y与所有解释变量x1，x2，……，xk之间的回归效果显著，即确实存在线性关系；反之，则不显著［4］。

对于多元回归模型，方程的显著性并不意味每个解释变量对被解释变量y的影响都是重要的。如果某个解释变量并不重要，则应该从方程中把它剔除，重新建立更为简单的方程。所以必须对每个解释变量进行显著性检验。

在给定的显著性水平a下，若｜ti｜＞ta/2（n-k-1），则拒绝H0，说明解释变量xi对被解释变量y有显著影响，即xi是影响y的主要因素；反之，接受H0，说明解释变量xi对被解释变量y无显著影响，则应删除该因素［4］。

3 应用实例

论文利用线形回归分析模型对南海海域海洋沉积物温度进行了分析，其中散点图显示如图1所示，回归分析结果见表1。

图1 水深与沉积物温度散点图

Fig.1 Water depth vs.sediment temperature

表1 水深沉积物温度回归分析结果Tab.1 The regression ysis result for Water depth vs.sediment temperature

截距：a=17.56；斜率：b=-0.0014；相关系数：R=0.276；测定系数：R2=0.076；F值：F=89.54。

建立回归模型，并对结果进行检验

南海地质研究（2014）

其中P＜0.0001。回归结果证明，沉积物温度与海水深度有着密切的关系，但是通过散点图显示，并不是温度越深沉积物温度越低。而是受到其他例如海底热流，海洋环流等因素的影响。

4 结语

本文介绍了回归分析在海洋地质调查研究中的应用，同时提供了回归分析的技术原理及实现方法，并通过对南海沉积物与海水深度关系模型进行了应用分析，回归结果显示了两者具有密切但是存在不确定性的关系。实验结果得到有效的应用。

参考文献

［1］单宝强，毛永强.2005.GIS中的坐标系定义与转换［J］.黑龙江国土资源，11，38-39

［2］苏国辉，孙记红，等.2011.海洋地质数据集成中的关键问题和方案［J］.海洋地质前沿，11（27）：51

［3］百度百科.回归分析.：//

［4］沈聪.2009.基于EXCEL的回归分析在足迹分析上的应用［M］.辽宁警官高等专科学校本科

［5］Cottrell A.Regression Analysis：Basic Concepts.：//

The Marine Geological Survey Based on Regression Analysis

Liang Guang1，2，Shao Changgao1，2

（1.Guangzhou Marine Geological Survey，Guangzhou，510760；2.Key Laboratory of Marine Mineral Reasources，MLR，Guangzhou，510760）

Abstract：The new resources survey had covered most of the continental area at present.So，the ocean resources he attracted more and more attention now as it is a huge resource and energy reservoir that had a profound meaning to national economy and military strategy.The energy competition made manly countries dloped new technology project and put the new ocean technology as the primary study area.Howr，how to abstract useful rmation from marine geological survey data is one of the most important study technologies.This focuses on the study of the deficit of marine database technology and introduces regression ysis model and the application aantage of it.The pure of this is to provide the technology support for marine study.Key word：Marine geology；Regression ysis model；Database

电气自动化工程师

求得［4］。

电气自动化工程师范本

Sen M K，and Stoffa P L.1992.Rapid sampling of model space using genetic algorithms：Examples from seiic weform inversion，Geophys.J.Internat.，108，281～292

姓名：网 目前所在：黄埔区年 龄：26 岁户口所在：湖南国 籍：婚姻状况：已婚民 族：汉族培训认证： 未参加 身 高：165 cm 诚信徽章： 未申请 体 重：60 kg 人才测评： 未测评 我的'特长： 求职意向 人才类型：普通求职 应聘职位：自动控制工程师

姓名：网

目前所在：黄埔区年龄：26 岁户口所在：湖南国籍：婚姻状况：已婚民族：汉族培训认证： 未参加 身高：165 cm 诚信徽章： 未申请 体重：60 kg 人才测评： 未测评 我的特长：

求职意向

人才类型：普通求职应聘职位：自动控制工程师/技术员：电气工程师，自动控制工程师/技术员：电气工程师，电气工程师/技术员：电工工作年限：3 职称：中级 求职类型： 可到职日期：随时 月薪要求：5000--8000 希望工作地区：广州,花都区,深圳

工作经历

广州海控自动化工程有限公司 起止年月：2012-10-10 ～ 2012-10-16公司性质：私营企业所属行业：其他行业担任职位：电气自动化工程师工作描述：自动化项目设计 安装 调试 熟悉西门子S7-200 S7-300PLC及PPI MPI DP 以太网通讯 熟悉应用三菱 施耐德 欧姆龙 富士等PLC 交流伺服 变频器 步进电机 软启动器应用 调试

工程例案:广州海洋地质调查局海洋四号万米绞车电气控制改造离职原因：上海凯泉泵业()有限公司 起止年月：2008-04-01 ～ 2012-09-01公司性质：民营企业所属行业：机械/机电/设备/重工担任职位：维修电工 电气工程师工作描述：恒压供水 熟悉西门子.施耐德.三菱.欧姆南海地质研究.2006龙PLC编程及变频器. 触摸屏应用

工程例案:武广高铁金沙洲隧道排污系统改造离职原因：回家广州励进新技术有限公司 起止年月：2007-08-01 ～ 2008-03-01公司性质：私营企业所属行业：机械/机电/设备/重工担任职位：助理工程师工作描述： PLC应用 氩弧焊机.三轴数控设备.铝塑复合管维修.调试离职原因：个人学习

教育背景

毕业院校：西安交通大学学历：本科获得学位: 毕业日期：2007-06-01专 业 一：工业自动化专 业 二：可编程序PLC 电子技术起始年月终止年月学校（机构）所学专业获得证书证书编号2001-09-012004-06-01新田二中理科高中毕业4329341412003-08-012007-06-01西安电气自动化656101040326792642007-03-012008-01-01广东工业大学维修电工高级电工/电工08190010032000052008-08-202008-11-01广东轻工学院CEAC电气智能技术应用电气智能技术应用高级工程师4500300

语言能力

工作能力及其他专长

自动化项目设计 安装 调试 熟悉西门子S7-200 S7-300PLC及MPI DP 以太网通讯 熟悉应用三菱 施耐德 欧姆龙 富士等PLC 交流伺服 变频器 步进电机 软启动器应用 调试

个人自我评价 工作积极 能吃苦耐劳 学习上进 团结同事 有错必改 (编辑：8021) ;