一种古代湖岸线的定量预测方法、设备及系统与流程

[0001]
本发明关于地质勘探技术领域，特别是关于古湖岸线的勘探技术，具体的讲是一种古代湖岸线的定量预测方法、系统、计算机设备以及计算机可读存储介质。


背景技术：

[0002]
目前，在对古代湖岸线进行位置判定主要包括以下几种方法：
[0003]
1)古地貌特征古冲沟及湖缘峡谷可以作为识别古湖岸线的指示标志；
[0004]
2)介(贝)壳滩是重要的滨岸标志物；
[0005]
3)水流交错层理、波痕、粉砂岩中的小型波状层理、泥岩中的氧化条带、植物根迹以及泥裂、雨痕等暴露构造(kom atsu et al.,2001；姜在兴,2003)、碳酸盐岩溶蚀带等；
[0006]
4)泥炭层、煤层、蒸发岩、沿着沿岸流方向成层分布的分选和磨圆较好的砾石层、砂岩层等，目前常采用岩心观察到的煤、炭质泥岩、含碳泥岩来判断湖岸发育的大致位置：
[0007]
5)地震反射特征，地震地层分析被广泛地应用于重建湖平面变化,这是因为这种变化控制了地震层序的横向展布、上超关系、剥蚀面和古湖岸线的产状(ansel mett,i 2006；rabineau et al.,2006；mchugh et al.,2008)；
[0008]
6)地球化学标志识别古湖岸线主要是通过化学指标判断沉积环境的氧化、还原条件和水介质性质,进而根据沉积环境恢复古湖岸线的位置(rosenet al.,2002)。
[0009]
以上方法都存在数据点获取困难，及数据较少等原因，很难在纵向及平面上相对准确的刻画，大多只能以少数点来勾勒大致形态，离实际应用要求相差很远。
[0010]
因此，如何提供一种新的方案，其能够克服上述缺陷是本领域亟待解决的技术难题。


技术实现要素：

[0011]
有鉴于此，本发明实施例提供了一种古代湖岸线的定量预测方法、系统、计算机设备以及计算机可读存储介质，根据对泥炭碳质泥岩层、煤层这样的特殊岩性识别来达到定量刻画湖岸线位置，不但解决了数据点资料获取难，也解决了小层级别湖岸线变迁演化过程的重现问题。
[0012]
本发明的目的之一是，提供一种古代湖岸线的定量预测方法，包括：
[0013]
获取沉积学及地震属性平面图；
[0014]
根据所述沉积学及地震属性平面图建立网格数据的预测样本集；
[0015]
构建沉积体类型判别模型；
[0016]
根据所述预测样本集通过所述沉积体类型判别模型识别沉积体类型。
[0017]
在本发明的优选实施方式中，对等时小层精细划分后的钻井进行岩心观察，得到岩心观察结果包括：
[0018]
根据区域研究结果，对等时小层精细划分后的钻井在三角洲平原与三角洲前缘亚相过渡带选取砂岩不发育的取心井；
[0019]
将所述取心井中观察到的煤层、炭质泥岩以及含碳泥岩挑选出来；
[0020]
记录所述煤层、炭质泥岩以及含碳泥岩各自的深度及对应的测井曲线特征，作为岩心观察结果。
[0021]
在本发明的优选实施方式中，根据所述岩心观察结果建立岩性区分模板包括：
[0022]
根据所述深度及对应的测井曲线特征得到所述煤层、炭质泥岩以及含碳泥岩各自的测井电性参数数据；
[0023]
根据所述测井电性参数数据进行平面投点分类，建立测井电性参数与煤层、炭质泥岩、含碳泥岩的区分模板。
[0024]
在本发明的优选实施方式中，根据所述岩性区分模板对钻井数据进行分析处理包括：
[0025]
根据煤层、炭质泥岩以及含碳泥岩的岩性成因及分布特征，剔除河道砂岩及滑塌岩发育的层段，选择出河道间及沼泽相的泥岩为对象数据；
[0026]
根据所述岩性区分模板确定出门槛值；
[0027]
根据所述门槛值在单井上计算识别出煤层、炭质泥岩以及含碳泥岩，计算并统计炭地比、炭厚度、炭频率及炭评价综合指数作为分析处理结果。
[0028]
在本发明的优选实施方式中，根据分析处理结果以及所述岩性区分模板预测古代湖岸线的位置信息包括：
[0029]
根据所述分析处理结果得到连井剖面上单井各小层中煤层、炭质泥岩以及含碳泥岩的累积厚度集中发育部位及频率变化；
[0030]
根据所述岩性区分模板建立数学模型；
[0031]
根据所述数学模型结合所述累积厚度集中发育部位及频率变化判别出古代湖岸线的位置信息。
[0032]
在本发明的优选实施方式中，所述方法还包括：
[0033]
获取所述位置信息附近的多个钻井进行岩心观察，得到钻井观察结果；
[0034]
对所述钻井观察结果进行检验，确定出为古湖岸线发育位置的吻合度；
[0035]
当所述吻合度达到预设的阈值时，将所述位置信息确定为古代湖岸线的位置信息。
[0036]
本发明的目的之一是，提供一种古代湖岸线的定量预测系统，包括：
[0037]
小层划分模块，用于对研究区内的钻井进行等时小层精细划分；
[0038]
岩心观察模块，用于对等时小层精细划分后的钻井进行岩心观察，得到岩心观察结果；
[0039]
模板建立模块，用于根据所述岩心观察结果建立岩性区分模板；
[0040]
分析处理模块，用于根据所述岩性区分模板对钻井数据进行分析处理；
[0041]
湖岸线预测模块，用于根据分析处理结果以及所述岩性区分模板预测古代湖岸线的位置信息。
[0042]
在本发明的优选实施方式中，所述岩心观察模块包括：
[0043]
取心井选取模块，用于根据区域研究结果，对等时小层精细划分后的钻井在三角洲平原与三角洲前缘亚相过渡带选取砂岩不发育的取心井；
[0044]
岩性挑选模块，用于将所述取心井中观察到的煤层、炭质泥岩以及含碳泥岩挑选
出来；
[0045]
数据记录模块，用于记录所述煤层、炭质泥岩以及含碳泥岩各自的深度及对应的测井曲线特征，作为岩心观察结果。
[0046]
在本发明的优选实施方式中，所述模板建立模块包括：
[0047]
参数数据确定模块，用于根据所述深度及对应的测井曲线特征得到所述煤层、炭质泥岩以及含碳泥岩各自的测井电性参数数据；
[0048]
区分模板建立模块，用于根据所述测井电性参数数据进行平面投点分类，建立测井电性参数与煤层、炭质泥岩、含碳泥岩的区分模板。
[0049]
在本发明的优选实施方式中，所述分析处理模块包括：
[0050]
层段剔除模块，用于根据煤层、炭质泥岩以及含碳泥岩的岩性成因及分布特征，剔除河道砂岩及滑塌岩发育的层段，选择出河道间及沼泽相的泥岩为对象数据；
[0051]
门槛确定模块，用于根据所述岩性区分模板确定出门槛值；
[0052]
识别计算模块，用于根据所述门槛值在单井上计算识别出煤层、炭质泥岩以及含碳泥岩，计算并统计炭地比、炭厚度、炭频率及炭评价综合指数作为分析处理结果。
[0053]
在本发明的优选实施方式中，所述湖岸线预测模块包括：
[0054]
变化确定模块，用于根据所述分析处理结果得到连井剖面上单井各小层中煤层、炭质泥岩以及含碳泥岩的累积厚度集中发育部位及频率变化；
[0055]
模型建立模块，用于根据所述岩性区分模板建立数学模型；
[0056]
湖岸线判别模块，用于根据所述数学模型结合所述累积厚度集中发育部位及频率变化判别出古代湖岸线的位置信息。
[0057]
在本发明的优选实施方式中，所述系统还包括：
[0058]
钻井获取模块，用于获取所述位置信息附近的多个钻井进行岩心观察，得到钻井观察结果；
[0059]
结果检验模块，用于对所述钻井观察结果进行检验，确定出为古湖岸线发育位置的吻合度；
[0060]
位置确定模块，用于当所述吻合度达到预设的阈值时，将所述位置信息确定为古代湖岸线的位置信息。
[0061]
本发明的目的之一是，提供一种计算机设备，包括：适于实现各指令的处理器以及存储设备，所述存储设备存储有多条指令，所述指令适于由处理器加载并执行一种古代湖岸线的定量预测方法。
[0062]
本发明的目的之一是，提供一种计算机可读存储介质，存储有计算机程序，所述计算机程序用于执行一种古代湖岸线的定量预测方法。
[0063]
本发明的有益效果在于，提供了一种古代湖岸线的定量预测方法、系统、计算机设备以及计算机可读存储介质，根据对泥碳质泥岩炭层、煤层这样的特殊岩性识别来达到定量刻画湖岸线位置，不但解决了数据点资料获取难，也解决了小层级别湖岸线变迁演化过程的重现问题。
[0064]
为让本发明的上述和其他目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附图式，作详细说明如下。
附图说明
[0065]
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0066]
图1为本发明实施例提供的一种古代湖岸线的定量预测系统的结构示意图；
[0067]
图2为本发明实施例提供的一种古代湖岸线的定量预测系统中岩心观察模块的结构示意图；
[0068]
图3为本发明实施例提供的一种古代湖岸线的定量预测系统中模板建立模块的结构示意图；
[0069]
图4为本发明实施例提供的一种古代湖岸线的定量预测系统中分析处理模块的结构示意图；
[0070]
图5为本发明实施例提供的一种古代湖岸线的定量预测系统中湖岸线预测模块的结构示意图；
[0071]
图6为本发明实施例提供的一种古代湖岸线的定量预测系统的实施方式二的结构示意图；
[0072]
图7为本发明实施例提供的一种古代湖岸线的定量预测方法的流程图；
[0073]
图8为图7中的步骤s102的具体流程图；
[0074]
图9为图7中的步骤s103的具体流程图；
[0075]
图10为图7中的步骤s104的具体流程图；
[0076]
图11为图7中的步骤s105的具体流程图；
[0077]
图12为本发明实施例提供的一种古代湖岸线的定量预测方法的实施方式二的流程图；
[0078]
图13为本发明提供的具体实施例中古代湖岸线的定量预测方法的流程图；
[0079]
图14(a)、图14(b)为本发明提供的具体实施例中鄂尔多斯盆地某地区延长组长81段岩性划分模板图；
[0080]
图15(a)、图15(b)为本发明提供的具体实施例中鄂尔多斯盆地某地区延长组长81段古湖岸线变迁模式剖面图。
具体实施方式
[0081]
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0082]
本领域技术人员知道，本发明的实施方式可以实现为一种系统、装置、方法或计算机程序产品。因此，本发明公开可以具体实现为以下形式，即：完全的硬件、完全的软件(包括固件、驻留软件、微代码等)，或者硬件和软件结合的形式。
[0083]
下面参考本发明的若干代表性实施方式，详细阐释本发明的原理和精神。
[0084]
本发明涉及古代湖岸线的定量预测方案，是针对目前根据少部分岩心观察结果及
相模式来人为定性判断湖岸线方法的改进及完善。
[0085]
具体的，本发明提供一种古代湖岸线的定量预测系统，图1为该系统的结构示意图，请参见图1，所述古代湖岸线的定量预测系统包括：
[0086]
小层划分模块100，用于对研究区内的钻井进行等时小层精细划分。
[0087]
在本发明的一种实施方式中，基于高分辨率层序地层学对研究区内井进行2～10米级等时小层精细划分，确保沉积期为等时。
[0088]
岩心观察模块200，用于对等时小层精细划分后的钻井进行岩心观察，得到岩心观察结果。
[0089]
图2为本发明实施例的岩心观察模块200的结构示意图，请参阅图2，岩心观察模块200包括：
[0090]
取心井选取模块210，用于根据区域研究结果，对等时小层精细划分后的钻井在对应小层内三角洲平原与三角洲前缘亚相过渡带选取砂岩不发育的取心井；
[0091]
岩性挑选模块220，用于将所述取心井中观察到的煤层、炭质泥岩以及含碳泥岩挑选出来；
[0092]
数据记录模块230，用于记录所述煤层、炭质泥岩以及含碳泥岩各自的深度及对应的测井曲线特征，作为岩心观察结果。
[0093]
也即，在本发明的一种实施方式中，根据区域研究结果，在平原与前缘亚相过渡带选取砂岩相对不发育的取心井，将钻井取心中观察到的煤层、炭质泥岩、含碳泥岩挑选出来，记录深度及对应测井曲线特征，主要进行煤层、炭质泥岩、含碳泥岩这三种岩性的识别与精细刻画，同时与测井曲线(gr+den)进行校正与归位。
[0094]
请参见图1，所述古代湖岸线的定量预测系统还包括：
[0095]
模板建立模块300，用于根据所述岩心观察结果建立岩性区分模板。
[0096]
图3为本发明实施例中模板建立模块300的结构示意图，请参见图3，所述模板建立模块300包括：
[0097]
参数数据确定模块310，用于根据所述深度及对应的测井曲线特征得到所述煤层、炭质泥岩以及含碳泥岩各自的测井电性参数数据；
[0098]
区分模板建立模块320，用于根据所述测井电性参数数据进行平面投点分类，建立测井电性参数与煤层、炭质泥岩、含碳泥岩的区分模板。
[0099]
也即，在本发明的一种实施方式中，根据岩心观察所得到的煤层、炭质泥岩、含碳泥岩三种岩性的gr、den、ac、cnl数据进行平面投点分类，选择区分效果最好的模板作为定量区分这三种岩性的参数，建立测井电性参数与煤层、炭质泥岩、含碳泥岩这三种岩性的区分模板。图14(a)、图14(b)为本发明提供的具体实施例中鄂尔多斯盆地某地区延长组长81段岩性划分模板图。
[0100]
请参见图1，所述古代湖岸线的定量预测系统还包括：
[0101]
分析处理模块400，用于根据所述岩性区分模板对钻井数据进行分析处理。
[0102]
图4为本发明实施例中分析处理模块400的结构示意图，请参见图4，所述分析处理模块400包括：
[0103]
层段剔除模块410，用于根据煤层、炭质泥岩以及含碳泥岩的岩性成因及分布特征，剔除河道砂岩及滑塌岩发育的层段，选择出河道间及沼泽相的泥岩为对象数据；
[0104]
门槛确定模块420，用于根据所述岩性区分模板确定出门槛值；
[0105]
识别计算模块430，用于根据所述门槛值在单井上计算识别出煤层、炭质泥岩以及含碳泥岩，计算并统计炭地比、炭厚度、炭频率及炭评价综合指数作为分析处理结果。
[0106]
也即，在本发明的一种实施方式中，首先剔除测井解释出的砂岩段，然后根据模板确定的门槛值，在单井上计算识别出煤层、炭质泥岩、含碳泥岩这三种岩性，计算并统计出“炭地比”、“炭厚度”、“炭频率”及“炭评价综合指数”。
[0107]
请参见图1，所述古代湖岸线的定量预测系统还包括：
[0108]
湖岸线预测模块500，用于根据分析处理结果以及所述岩性区分模板预测古代湖岸线的位置信息。
[0109]
图5为本发明实施例中湖岸线预测模块500的结构示意图，请参见图5，所述湖岸线预测模块500包括：
[0110]
变化确定模块510，用于根据所述分析处理结果得到连井剖面上单井各小层中煤层、炭质泥岩以及含碳泥岩的累积厚度集中发育部位及频率变化；
[0111]
模型建立模块520，用于根据所述岩性区分模板建立数学模型；
[0112]
湖岸线判别模块530，用于根据所述数学模型结合所述累积厚度集中发育部位及频率变化判别出古代湖岸线的位置信息。
[0113]
也即，在本发明的一种实施方式中，根据测井计算出的煤层、炭质泥岩、含碳泥岩这三种岩性的累积厚度及频率，根据数学模型、三种岩性的累积厚度及频率得出的古湖岸线位置。
[0114]
图6为本发明实施例提供的一种古代湖岸线的定量预测系统的实施方式二的结构示意图，请参阅图6，该系统在本发明的实施方式二中还包括：
[0115]
钻井获取模块600，用于获取所述位置信息附近的多个钻井进行岩心观察，得到钻井观察结果；
[0116]
结果检验模块700，用于对所述钻井观察结果进行检验，确定出为古湖岸线发育位置的吻合度；
[0117]
位置确定模块800，用于当所述吻合度达到预设的阈值时，将所述位置信息确定为古代湖岸线的位置信息。
[0118]
也即，在本发明的一种实施方式中，选取古代湖岸线的位置信息附近的钻井进行岩心观察，根据实际观察结果判定其是否为古湖岸线发育位置，如果检验结果吻合度超过预设的阈值(诸如70％)，则进行确定结果，若吻合度低于预设的阈值(诸如70％)，则重新回到模板建立模块，换另一批井重新进行。
[0119]
如上即为本发明提供的一种古代湖岸线的定量预测系统，主要是根据对泥炭碳质泥岩层、煤层这样的特殊岩性识别来达到定量刻画湖岸线位置，不但解决了数据点资料获取难，也解决了小层级别湖岸线变迁演化过程的重现问题。
[0120]
此外，尽管在上文详细描述中提及了系统的若干单元模块，但是这种划分仅仅并非强制性的。实际上，根据本发明的实施方式，上文描述的两个或更多单元的特征和功能可以在一个单元中具体化。同样，上文描述的一个单元的特征和功能也可以进一步划分为由多个单元来具体化。以上所使用的术语“模块”和“单元”，可以是实现预定功能的软件和/或硬件。尽管以下实施例所描述的模块较佳地以软件来实现，但是硬件，或者软件和硬件的组
合的实现也是可能并被构想的。
[0121]
在介绍了本发明示例性实施方式的古代湖岸线的定量预测系统之后，接下来，参考附图对本发明示例性实施方式的方法进行介绍。该方法的实施可以参见上述整体的实施，重复之处不再赘述。
[0122]
本发明涉及古代湖岸线的定量预测方案，通过煤层、炭质泥岩、含碳泥岩三种岩性的横向及纵向发育厚度及频率，判断古湖岸线发育带，从而判别湖岸线的演化规律。
[0123]
具体的，本发明提供一种古代湖岸线的定量预测方法，图7为该方法的流程示意图，请参见图7，所述古代湖岸线的定量预测方法包括：
[0124]
s101：对研究区内的钻井进行等时小层精细划分。
[0125]
在本发明的一种实施方式中，基于高分辨率层序地层学对研究区内井进行2～10米级等时小层精细划分，确保沉积期为等时。
[0126]
s102：对等时小层精细划分后的钻井进行岩心观察，得到岩心观察结果。
[0127]
图8为本发明实施例的步骤s102的流程示意图，请参阅图8，步骤s102包括：
[0128]
s201：根据区域研究结果，对等时小层精细划分后的钻井在对应小层内三角洲平原与三角洲前缘亚相过渡带选取砂岩不发育的取心井；
[0129]
s202：将所述取心井中观察到的煤层、炭质泥岩以及含碳泥岩挑选出来；
[0130]
s203：记录所述煤层、炭质泥岩以及含碳泥岩各自的深度及对应的测井曲线特征，作为岩心观察结果。
[0131]
也即，在本发明的一种实施方式中，根据区域研究结果，在平原与前缘亚相过渡带选取砂岩相对不发育的取心井，将钻井取心中观察到的煤层、炭质泥岩、含碳泥岩挑选出来，记录深度及对应测井曲线特征，主要进行煤层、炭质泥岩、含碳泥岩这三种岩性的识别与精细刻画，同时与测井曲线(gr+den)进行校正与归位。
[0132]
请参见图7，所述古代湖岸线的定量预测方法还包括：
[0133]
s103：根据所述岩心观察结果建立岩性区分模板。
[0134]
图9为本发明实施例中步骤s103的流程示意图，请参见图9，该步骤包括：
[0135]
s301：根据所述深度及对应的测井曲线特征得到所述煤层、炭质泥岩以及含碳泥岩各自的测井电性参数数据；
[0136]
s302：根据所述测井电性参数数据进行平面投点分类，建立测井电性参数与煤层、炭质泥岩、含碳泥岩的区分模板。
[0137]
也即，在本发明的一种实施方式中，根据岩心观察所得到的煤层、炭质泥岩、含碳泥岩三种岩性的gr、den、ac、cnl数据进行平面投点分类，选择区分效果最好的模板作为定量区分这三种岩性的参数，建立测井电性参数与煤层、炭质泥岩、含碳泥岩这三种岩性的区分模板。图14(a)、图14(b)为本发明提供的具体实施例中鄂尔多斯盆地某地区延长组长81段岩性划分模板图。
[0138]
请参见图7，所述古代湖岸线的定量预测方法还包括：
[0139]
s104：根据所述岩性区分模板对钻井数据进行分析处理。
[0140]
图10为本发明实施例中步骤s104的流程示意图，请参见图10，该步骤包括：
[0141]
s401：根据煤层、炭质泥岩以及含碳泥岩的岩性成因及分布特征，剔除河道砂岩及滑塌岩发育的层段，选择出河道间及沼泽相的泥岩为对象数据；
[0142]
s402：根据所述岩性区分模板确定出门槛值；
[0143]
s403：根据所述门槛值在单井上计算识别出煤层、炭质泥岩以及含碳泥岩，计算并统计炭地比、炭厚度、炭频率及炭评价综合指数作为分析处理结果。
[0144]
也即，在本发明的一种实施方式中，首先剔除测井解释出的砂岩段，然后根据模板确定的门槛值，在单井上计算识别出煤层、炭质泥岩、含碳泥岩这三种岩性，计算并统计出“炭地比”、“炭厚度”、“炭频率”及“炭评价综合指数”。
[0145]
请参见图7，所述古代湖岸线的定量预测方法还包括：
[0146]
s105：根据分析处理结果以及所述岩性区分模板预测古代湖岸线的位置信息。
[0147]
图11为本发明实施例中步骤s105的流程示意图，请参见图11，该步骤s105包括：
[0148]
s501：根据所述分析处理结果得到连井剖面上单井各小层中煤层、炭质泥岩以及含碳泥岩的累积厚度集中发育部位及频率变化；
[0149]
s502：根据所述岩性区分模板建立数学模型；
[0150]
s503：根据所述数学模型结合所述累积厚度集中发育部位及频率变化判别出古代湖岸线的位置信息。
[0151]
也即，在本发明的一种实施方式中，根据测井计算出的煤层、炭质泥岩、含碳泥岩这三种岩性的累积厚度及频率，根据数学模型、三种岩性的累积厚度及频率得出的古湖岸线位置。
[0152]
图12为本发明实施例提供的一种古代湖岸线的定量预测方法的实施方式二的流程示意图，请参阅图12，该方法在本发明的实施方式二中还包括：
[0153]
s106：获取所述位置信息附近的多个钻井进行岩心观察，得到钻井观察结果；
[0154]
s107：对所述钻井观察结果进行检验，确定出为古湖岸线发育位置的吻合度；
[0155]
s108：当所述吻合度达到预设的阈值时，将所述位置信息确定为古代湖岸线的位置信息。
[0156]
也即，在本发明的一种实施方式中，选取古代湖岸线的位置信息附近的钻井进行岩心观察，根据实际观察结果判定其是否为古湖岸线发育位置，如果检验结果吻合度超过预设的阈值(诸如70％)，则进行确定结果，若吻合度低于预设的阈值(诸如70％)，则重新回到模板建立模块，换另一批井重新进行。
[0157]
如上即为本发明提供的一种古代湖岸线的定量预测方法，主要是根据对泥炭碳质泥岩层、煤层这样的特殊岩性识别来达到定量刻画湖岸线位置，不但解决了数据点资料获取难，也解决了小层级别湖岸线变迁演化过程的重现问题。
[0158]
应当注意，尽管在附图中以特定顺序描述了本发明方法的操作，但是，这并非要求或者暗示必须按照该特定顺序来执行这些操作，或是必须执行全部所示的操作才能实现期望的结果。附加地或备选地，可以省略某些步骤，将多个步骤合并为一个步骤执行，和/或将一个步骤分解为多个步骤执行。
[0159]
本发明还提供了一种计算机设备，包括：适于实现各指令的处理器以及存储设备，所述存储设备存储有多条指令，所述指令适于由处理器加载并执行一种古代湖岸线的定量预测方法。
[0160]
本发明还提供了一种计算机可读存储介质，存储有计算机程序，所述计算机程序用于执行一种古代湖岸线的定量预测方法。
[0161]
下面结合具体的实施例，详细介绍本发明的技术方案。图13为本发明提供的具体实施例中古代湖岸线的定量预测方法的流程图，若研究区域较大，则岩心观察井需在数量及平面分布上有保证，请参阅图13，在该实施例中该方法包括：
[0162]
1、对研究区钻井进行小层对比。
[0163]
基于高分辨率层序地层学对研究区内井进行2～10米级等时小层精细划分。
[0164]
2、根据钻井分层结果进行岩心观察。
[0165]
将钻井取心中观察到的煤层、炭质泥岩、含碳泥岩挑选出来，记录深度及对应测井曲线特征。
[0166]
3、根据岩心观察结果制作特殊岩性模板。
[0167]
根据岩心观察结果，建立测井电性参数与煤层、炭质泥岩、含碳泥岩这三种岩性的区分模板，图14(a)、图14(b)为本发明提供的具体实施例中鄂尔多斯盆地某地区延长组长81段岩性划分模板图，研究区内根据岩心观察确定的岩性，其对应的参数进行投点划分后所得结果，目的是为了在未进行岩心观察的井上定量区分出煤层、炭质泥岩、含碳泥岩这三种岩性。
[0168]
4、根据模板参数建立小层级别湖岸线定量识别数学模型；
[0169][0170]
其中，hax为胡岸线定量识别综合指数，ct为炭质岩厚度，ct/s为炭质岩与地层厚度比值；cf为炭质岩频率。
[0171]
5、根据数学模型判断古湖岸线位置。
[0172]
根据煤层、炭质泥岩、含碳泥岩这三种岩性成因及分布特征，首先剔除河道砂岩及滑塌岩发育的层段，优先选择河道间及沼泽相的泥岩为对象数据，然后根据模板确定的门槛值，在单井上计算识别出煤层、炭质泥岩、含碳泥岩这三种岩性，计算并统计出“炭地比”、“炭厚度”、“炭频率”及“炭评价综合指数”。
[0173]
6、与判断的湖岸线位置附近钻井岩心观察结果相对比进行结果验证。
[0174]
根据测井计算出的煤层、炭质泥岩、含碳泥岩这三种岩性的累积厚度及频率，湖岸带(线)判定结果检验：最后根据连井剖面上单井各小层三种岩性的累积厚度集中发育部位及频率变化结合数学模型判定古湖岸线位置。
[0175]
选取其附近的钻井进行岩心观察，根据实际观察结果判定其是否为古湖岸线发育位置，例如，通过古生物或古代浅水植物等特征，确定古湖岸线发育位置后与数学模型计算得出的结果进行校正，校正后的综合指数结果为该研究区小层内的湖岸线定量识别标准，如果新井综合指数大于标准指数(如检验结果吻合度超过70％)，则进行到第7步，确定结果，若低于标准指数(如吻合度低于70％)，则重新回到第三步，换另一批井重新进行。
[0176]
7、确定结果。
[0177]
图15(a)、图15(b)为本发明提供的具体实施例中鄂尔多斯盆地某地区延长组长81段古湖岸线变迁模式剖面图。
[0178]
由该实施例可见，在鄂尔多斯盆地陇东地区延长组长81段，通过该方法精细刻画了该区三个小层对应三期湖岸线的空间分布特征，为该区沉积及储层演化过程的研究奠定
了基础。
[0179]
综上所述，本发明提供的一直古代湖岸线的定量预测方法、系统、计算机设备以及计算机可读存储介质，创新性发现通过煤层、泥炭碳质泥岩、含碳泥岩这三种特殊岩性的频率、累积厚度、炭地比形成的综合指数可以判断古湖岸线(带)位置，也可以反映古湖盆边部底形的陡缓，不但解决了数据点资料获取难，也解决了小层级别湖岸线变迁演化过程的重现问题。
[0180]
本发明的难点在于岩心观察确定实际岩性后，对应的测井曲线归位及模板识别的准确度，创新点在于建立了由煤层、泥炭碳质泥岩、含碳泥岩这三种特殊岩性的频率、累积厚度、炭地比形成的综合指数模型，需要一定数量及相对均匀分布的钻井岩心观察来代表区域特征。
[0181]
对于一个技术的改进可以很明显地区分是硬件上的改进(例如，对二极管、晶体管、开关等电路结构的改进)还是软件上的改进(对于方法流程的改进)。然而，随着技术的发展，当今的很多方法流程的改进已经可以视为硬件电路结构的直接改进。设计人员几乎都通过将改进的方法流程编程到硬件电路中来得到相应的硬件电路结构。因此，不能说一个方法流程的改进就不能用硬件实体模块来实现。例如，可编程逻辑器件(programmable logic device,pld)(例如现场可编程门阵列(field programmable gate array，fpga))就是这样一种集成电路，其逻辑功能由用户对器件编程来确定。由设计人员自行编程来把一个数字系统“集成”在一片pld上，而不需要请芯片制造厂商来设计和制作专用的集成电路芯片。而且，如今，取代手工地制作集成电路芯片，这种编程也多半改用“逻辑编译器(logic compiler)”软件来实现，它与程序开发撰写时所用的软件编译器相类似，而要编译之前的原始代码也得用特定的编程语言来撰写，此称之为硬件描述语言(hardware description language，hdl)，而hdl也并非仅有一种，而是有许多种，如abel(advanced boolean expression language)、ahdl(altera hardware description language)、confluence、cupl(cornell university programming language)、hdcal、jhdl(java hardware description language)、lava、lola、myhdl、palasm、rhdl(ruby hardware description language)等，目前最普遍使用的是vhdl(very-high-speed integrated circuit hardware description language)与verilog2。本领域技术人员也应该清楚，只需要将方法流程用上述几种硬件描述语言稍作逻辑编程并编程到集成电路中，就可以很容易得到实现该逻辑方法流程的硬件电路。
[0182]
控制器可以按任何适当的方式实现，例如，控制器可以采取例如微处理器或处理器以及存储可由该(微)处理器执行的计算机可读程序代码(例如软件或固件)的计算机可读介质、逻辑门、开关、专用集成电路(application specific integrated circuit，asic)、可编程逻辑控制器和嵌入微控制器的形式，控制器的例子包括但不限于以下微控制器：arc 625d、atmel at91sam、microchip pic18f26k20以及silicone labs c8051f320，存储器控制器还可以被实现为存储器的控制逻辑的一部分。
[0183]
本领域技术人员也知道，除了以纯计算机可读程序代码方式实现控制器以外，完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件，而对其内包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者
甚至，可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
[0184]
上述实施例阐明的系统、装置、模块或单元，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。
[0185]
为了描述的方便，描述以上装置时以功能分为各种单元分别描述。当然，在实施本申请时可以把各单元的功能在同一个或多个软件和/或硬件中实现。
[0186]
通过以上的实施方式的描述可知，本领域的技术人员可以清楚地了解到本申请可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在存储介质中，如rom/ram、磁碟、光盘等，包括若干指令用以使得一台计算机系统(可以是个人计算机，服务器，或者网络系统等)执行本申请各个实施例或者实施例的某些部分所述的方法。
[0187]
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
[0188]
本申请可用于众多通用或专用的计算机系统环境或配置中。例如：个人计算机、服务器计算机、手持系统或便携式系统、平板型系统、多处理器系统、基于微处理器的系统、置顶盒、可编程的消费电子系统、网络pc、小型计算机、大型计算机、包括以上任何系统或系统的分布式计算环境等等。
[0189]
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请，在这些分布式计算环境中，由通过通信网络而被连接的远程处理系统来执行任务。在分布式计算环境中，程序模块可以位于包括存储系统在内的本地和远程计算机存储介质中。
[0190]
虽然通过实施例描绘了本申请，本领域普通技术人员知道，本申请有许多变形和变化而不脱离本申请的精神，希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。