多层密度界面的拟BP神经网络反演方法教材.doc

多层密度界面的拟BP神经网络反演方法 刘展1 赵文举2 相鹏1 中国石油大学（华东）地球资源与信息学院，东营，257061 东方地球物理勘探有限责任公司综合物化探事业部，涿州，072751） 摘要 提出一种根据重力异常反演三维密度界面分布的反演模式。该模式将拟BP神经网络与重力反演理论结合，与传统神经网络不同的是拟BP神经网络不需要进行训练，而是直接求取隐层中的物性值。该模式应用于合成数据集可以发现三维密度界面能够被很好的复原。最后，利用该方法反演了冲绳海槽南部第三系底与莫霍面深度。 关键词 三维重力反演，密度界面，拟BP神经网络，冲绳海槽南，莫霍面，第三系基底1、引言 根据重力异常求取三维密度界面的几何形态是重力数据解释工作的一个主要目标。目前存在很多种不同的算法，例如，Oldenburg（1974）对Parker（1973）提出的非均匀层状介质正演公式重新推导得到根据已知重力异常求取密度界面深度的反演公式。Rao等（1999）利用邻接直立棱柱体模型根据重力异常或者基底构造求取三维密度界面深度。 人工神经网络已经被成功地应用于地球物理数据处理和反演问题当中。例如，测井数据解释（Wiener等，1991；Huang等1996），反射地震数据处理（Ashida，1996），近地表电磁成像模式识别（Poulton等，1992），密度界面反演（Taylor，Vasco，1991；朱自强，1995）等。尽管取得了进步，但是众所周知的是神经网络的反演结果很大程度上取决于训练数据，所以训练数据的选择是决定神经网络性能的关键。当训练数据与观测数据的模式存在很大差异时，神经网络会求出不合理的结果。管志宁（1998）将BP神经网络与重磁异常反演相结合提出了一种新的反演算法—拟BP神经网络。隐层中的物性值可以直接求出而不需要传统神经网络的训练过程。 本文提出一种迭代拟BP神经网络三维密度界面反演算法。首先，简要回顾一下三维密度界面的正演模型；然后详细介绍拟BP神经网络三维密度界面反演算法；接着利用合成数据分析算法的性能；最后用该方法求取南冲绳海槽盆地的第三系基底和莫霍面深度。 2、方法原理 2.1 三维密度界面正演方法 对于三维密度界面，我们建立如图1所示的坐标系以及模型剖分方案，采用一组平行于X轴及Y轴的等间距(和)平面，将各密度层划分为数个小长方体的组合。只要这种剖分足够细，就可以精细地拟合该三维模型。 图1 三维密度界面模型及剖分方案 采样点位置，剩余密度 对图1所示的模型，我们设X方向上的计算间距、计算点数分别为dx和N，Y方向上的计算间距、计算点数分别为dy和M；X方向上的长方体宽度、长方体个数分别为和J，Y方向上的长方体宽度、长方体个数分别为和I；X方向上的计算点编号为n=0,1,…,N-1；Y方向上的计算点编号为m=0,1,…,M-1；X方向上的长方体编号为j=0,1,…,J-1；Y方向上的长方体编号为i=0,1,…,I-1；第一计算点坐标为(dx/2,dy/2,0)；界面数为K，界面编号为k=1,2,…,K；各密度界面的剩余密度为,,…,；则编号为nm的计算点处的重力异常可用如下离散公式计算： (1) 式中，为第k层编号为ij的矩形下底的深度，为第k层编号为ij的矩形上底的深度，且。 这样，当n从0变化到N-1，m从0变化到M-1时，利用公式（1）便可计算出所有计算点的重力异常△。 2.2 密度界面的拟BP神经网络反演算法 根据BP神经网络及其网络映射定理，任意一连续函数（映射）都可以由三层网络来实现。因此，选择拟BP神经网络由3层组成（图2），即包括输入层、隐层和输出层，取输入层与输出层的神经元个数相等，即等于参加反演计算的测点总数。设三维空间划分单元的底面深度为隐层神经元，故隐层的神经元个数即等于三维空间划分单元的总数。若将实测重力异常场值加到输入层各神经元，输出层输出即为各测点的理论重力异常值（隐层各划分单元正演的场值），并将输出层的输出反馈到输入层，以便进行实测场与理论场的比较(管志宁,1998)。 图2 密度界面反演的拟BP神经网络结构 对于如图1所示的三维密度界面模型，设第kij个单元的第t次迭代深度值为，第t+1次迭代的深度值为，其调整值为，则有 ， （2） 设和分别为实测重力异常和理论重力异常，对和分别作S型函数变换，可得 （3） （4） 为神经元的阈值，为控制系数。 定义网络的误差函数为 （5） 由（5）式，对E求