[0058] 对于复杂地层，如褶皱、互层、地层倒转等，为得到 值，需进行多次地层相对关系的判断，而后求和：
[0059]
[0060] 式中， 表示第k号钻孔中第a次遇到某一组地层时的判断值； 值为1时，表示地层i位于地层j之上；值为‑1时，表示地层i位于地层j之下；值为0时，取 例如地层为CDC，则 求和后为0，此时应取
[0061] (3)综合所有钻孔数据和地层种类数对各地层逐一进行层序判断，获得各个钻孔中地层之间的位置关系值 以及所有钻孔数据得出的各地层之间的位置关系值Cij，并记录计算结果。例如地层A和地层B的位移关系值判断过程如下：
[0062] 钻孔1中A位于B之上，则
[0063] 钻孔2中A缺失，则
[0064] 钻孔3中A缺失，则
[0065] 钻孔4中A位于B之上，则
[0066]
[0067] 同理，依次可以计算其它各地层之间的位置关系值，至此，各元素Cij值已确定，更新上述二维数组，获得地层关系表，该地层关系表如下表1所示。
[0068] 表1地层关系表
[0069] CAB CAC CAD CAE CBC CBD CBE CCD CCE CDE2 2 1 ‑1 3 3 1 1 0 2
[0070] 更新后的二维数组为：
[0071]
[0072] (4)经过上述比较可初步确定各地层的相对位置关系，但之间可能存在矛盾。例如CAB＝4及CCE＝4说明A→B→E，而CAE＝4则E→A，与之前的A→B→E相互矛盾。此时以相同地层层序出现频率高的原则，依据Cij绝对值进行如下处理：
[0073] (4.1)依次取|Cij|较大值有条件地存入一个地层层序容器 中，当欲新加入的层序Cij与 中已有层序存在矛盾，则该Cij舍弃，否则加入，并将Cij值置为0，重复本步骤，可以得到更新后的地层关系表，如表2所示，该表中排除了地层层序相互矛盾的情况。
[0074] 表2更新后的地层关系表
[0075] CBC CBD CAB CAC CDE CAD CBE CCD3 3 2 2 2 1 1 1
[0076] (4.2)任意初始化地层层序，如ECBAD；
[0077] (4.3)依据上述更新后的地层关系表，对ECBAD的层序进行调整，调整过程如下所示：
[0078]
[0079]
[0080] (4.4)最终更新后的获得区域内地层层序表，上述中地层层序确定为ABCDE，并结果是唯一的。
[0081] S3、基于所述地层层序表，识别区域内复杂地质体形式和尖灭部位。
[0082] 该步骤在生成的地层层序表基础上，进行复杂地质体的判断，并优先生成复杂地质体，保证三维模型的准确性。
[0083] 如图3为例，钻孔Z1数据出现明显的CDC排列，即以D为中心上下两层地层相同；钻孔Z2数据出现CDCBC排列，即分别以D、B为中心，上下地层相同；钻孔Z3数据出现ABCBAB排列，即分别以C、A为中心，上下地层相同。图4为图3所对应的真实地层，该图中可以看出有典型的对称结构，如CDC、ABCBAB等。
[0084] (1)寻找钻孔中的对称序列，反映在数据形式上，以Z1为例，即为：数据互为相反数。
[0085] (2)寻找地层尖灭部位，以图3为例，Z1、Z2钻孔中均存在以D为中心的对称结构，上下地层均为C，Z3钻孔中不再出现，反映在数据形式上即为：即D地层尖灭与Z2与Z3钻孔之间。
[0086] S4、基于所述区域内复杂地质体的形式及尖灭部位，并结合所述地层层序表对复杂地质体进行优先建模，逐渐向周边扩展模型范围，直至生成完整的地质模型。
[0087] (1)以尖灭处为中心，向周边扩展，根据所述复杂地质体形式和地层层序表，采用空间插值法建模，优先生成复杂地质体，本实施例中的空间插值法为克里金插值方法。
[0088] (2)从复杂地质体向四周依次建立各地层的三维模型，并集成地层名称、力学参数等属性参数，将各地层模型组合为完整的地质模型。
[0089] 此外，地质建模时空间插值法也可选择三角剖分法、距离倒数加权法、最小曲率法、最近邻点法等。
[0090] 上述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read‑Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
[0091] 本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD‑ROM、光学存储器等)上实施的计算机程序产品的形式。本发明实施例中的方案可以采用各种计算机语言实现，例如，面向对象的程序设计语言Java和直译式脚本语言JavaScript等。
[0092] 本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
[0093] 这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
[0094] 这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
[0095] 尽管已描述了本发明的优选实施例，但本领域内的技术人员一旦得知了基本创造性概念，则可对这些实施例作出另外的变更和修改。所以，所附权利要求意欲解释为包括优选实施例以及落入本发明范围的所有变更和修改。
[0096] 显然，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也意图包含这些改动和变型在内。
