一种页岩薄储层的预测方法、装置、设备及存储介质与流程

1.本发明涉及地球物理勘探领域，特别地，涉及一种页岩薄储层的预测方法、装置、设备及存储介质。


背景技术：

2.随着油气勘探理论的不断深入和开发技术的不断提高，页岩气作为一种非常规的天然气资源，已经成为我国天然气储量增长的一个十分重要的勘探领域，展示了其巨大的资源潜力。伴随我国页岩勘探开发的不断推进，海相页岩气已实现规模开发，但海陆过渡相页岩气仍处在探索之中，我国海陆过渡相页岩气资源丰富，勘探潜力大，是我国页岩气勘探的一个重要领域。
3.然而海陆过渡相页岩层多与煤层交互发育，页岩层往往发育于厚煤层之下，页岩层厚度薄、非均质性强。用常规的地震预测方法进行页岩薄储层预测的过程中，由于煤层的强振幅地震反射特征，致使其下伏的页岩薄储层地震反射很弱甚至为空白地震反射特征，从而干扰页岩薄储层的识别和预测，给此类页岩的勘探和开发造成了很大的困难。因此，用常规的地震预测方法来预测页岩薄储层的准确度低且难度大，亟需一种针对海陆过渡相煤层强屏蔽下的页岩薄储层的预测方法。


技术实现要素：

4.本文实施例的目的在于提供一种页岩薄储层的预测方法、装置、设备及存储介质，以提高对煤层强屏蔽下页岩薄储层预测的准确度。
5.为达到上述目的，一方面，本文实施例提供了一种页岩薄储层的预测方法，所述方法包括：
6.获取工区内每口井对应的钻井数据、多种测井曲线、所述工区内目的层的地震解释层位数据，以及所述工区内的地震纯波偏移数据；
7.将所述钻井数据和所述多种测井曲线进行综合分析，确定所述工区内目的层段页岩薄储层的敏感参数及其取值范围；
8.利用所述目的层的地震解释层位数据、各层位的接触关系和所述敏感参数对应的测井曲线建立地震波形指示模拟的初始模型；
9.通过所述初始模型、所述地震纯波偏移数据和所述敏感参数对应的测井曲线进行模拟参数测试，确定最优模拟参数；
10.通过所述最优模拟参数进行地震波形指示模拟，确定所述工区内目的层段页岩薄储层的敏感参数模拟数据体；
11.提取所述敏感参数模拟数据体中位于敏感参数的取值范围内的数据，并将其确定为所述页岩薄储层的剖面预测结果。
12.优选的，所述利用所述目的层的地震解释层位数据、各层位的接触关系和所述敏感参数对应的测井曲线建立地震波形指示模拟的初始模型之前，还包括：
13.通过合成地震记录井—震标定方法，将所述多种测井曲线从深度域转换为时间域，为地震波形指示模拟的初始模型的建立提供准确的时深关系。
14.优选的，所述目的层的地震解释层位数据，包括：
15.通过所述合成地震记录井—震标定方法，确定目的层顶、底的地震解释层位数据；
16.或，通过所述合成地震记录井—震标定方法，确定目的层顶、底的地震解释层位数据以及目的层内部的地震解释层位数据。
17.优选的，所述将所述钻井数据和所述多种测井曲线进行综合分析，确定所述工区内目的层段页岩薄储层的敏感参数及其取值范围，包括：
18.将所述钻井数据和所述多种测井曲线进行综合分析，确定页岩薄储层对应的多种测井曲线的特征；
19.根据页岩薄储层对应的多种测井曲线的特征，确定每口井内页岩薄储层的分布；
20.将所述工区内目的层段所有井对应的多种测井曲线进行交会分析，将交会分析的结果与页岩薄储层的分布对应，确定所述工区内页岩薄储层的敏感参数及其取值范围。
21.优选的，所述工区内的多口井均匀分布在工区内。
22.优选的，所述方法还包括：
23.根据所述页岩薄储层的所述剖面预测结果，计算所述工区平面上各点的所述页岩薄储层的累计时间厚度；
24.将所述工区平面上各点的所述页岩薄储层的累计时间厚度进行平面成图，获得页岩薄储层的平面预测结果。
25.优选的，所述计算所述工区平面上各点的所述页岩薄储层的累计时间厚度，包括：
26.将平面上各点处每套所述页岩薄储层的底层和顶层对应的时间值相减，得到该点处每套所述页岩薄储层的时间厚度；
27.将各点对应的每套所述页岩薄储层的时间厚度累加，得到该点处页岩薄储层的累计时间厚度。
28.另一方面，本文实施例提供了一种页岩薄储层的预测装置，所述装置包括：
29.获取模块：获取工区内每口井对应的钻井数据、多种测井曲线、所述工区内目的层的地震解释层位数据，以及所述工区内的地震纯波偏移数据；
30.敏感参数确定模块：将所述钻井数据和所述多种测井曲线进行综合分析，确定所述工区内目的层段页岩薄储层的敏感参数及其取值范围；
31.模型建立模块：利用所述目的层的地震解释层位数据、各层位的接触关系和所述敏感参数对应的测井曲线建立地震波形指示模拟的初始模型；
32.最优模拟参数确定模块：通过所述初始模型、所述地震纯波偏移数据和所述敏感参数对应的测井曲线进行模拟参数测试，确定最优模拟参数；
33.敏感参数模拟数据体确定模块：通过所述最优模拟参数进行地震波形指示模拟，确定所述工区内目的层段页岩薄储层的敏感参数模拟数据体；
34.结果获取模块：提取所述敏感参数模拟数据体中位于敏感参数的取值范围内的数据，并将其确定为所述页岩薄储层的剖面预测结果。
35.又一方面，本文实施例还提供了一种计算机设备，包括存储器、处理器、以及存储在所述存储器上的计算机程序，所述计算机程序被所述处理器运行时，执行根据上述任意
一项所述方法的指令。
36.又一方面，本文实施例还提供了了一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被计算机设备的处理器运行时，执行根据上述任意一项所述方法的指令。
37.由以上本文实施例提供的技术方案可见，本文实施例通过将钻井数据和多种测井曲线数据进行对比分析，即可确定工区内页岩薄储层的敏感参数及其取值范围，其中敏感参数对应的测井曲线可以反映某一储层为页岩薄储层，敏感参数的取值范围可以作为区分页岩薄储层和其他储层的依据，进而能够通过敏感参数及其取值范围将页岩薄储层和其他储层区分开。
38.在地震波形指示模拟的初始模型建立后，通过模拟参数测试，确定最优模拟参数，通过最优模拟参数进行地震波形指示模拟，确定工区内的敏感参数模拟数据体。提取敏感参数模拟数据体中位于敏感参数的取值范围内的数据，通过敏感参数及其取值范围能够将页岩薄储层取出，该数据即为页岩薄储层的剖面预测结果，进而识别出煤层强屏蔽下的页岩薄储层。解决了对页岩薄储层进行预测时准确度较低的问题。因此，本文实施例提高了对页岩薄储层进行预测时的准确度。
39.为让本文的上述和其他目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附图式，作详细说明如下。
附图说明
40.为了更清楚地说明本文实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本文的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
41.图1示出了本文实施例提供的一种页岩薄储层的预测方法的流程示意图；
42.图2示出了本文实施例提供的井—震标定的步骤示意图；
43.图3示出了本文实施例提供的确定页岩薄储层的敏感参数及其取值范围的步骤示意图；
44.图4示出了本文实施例提供的对页岩薄储层进行平面成图的步骤示意图；
45.图5示出了本文实施例提供的计算时间厚度的步骤示意图；
46.图6示出了本文实施例提供的用于体现目的层其中一口井内的多种测井曲线图；
47.图7示出了本文实施例提供的自然伽马和补偿中子的交会图；
48.图8示出了本文实施例提供的用于体现目的层的地震反射示意图；
49.图9示出了本文实施例提供的初始模型示意图；
50.图10示出了本文实施例提供的用于体现页岩薄储层的剖面预测结果的示意图；
51.图11示出了本文实施例提供的页岩薄储层的平面图；
52.图12示出了本文实施例提供的一种页岩薄储层的预测装置的模块结构示意图；
53.图13示出了本文实施例提供的计算机设备的结构示意图。
54.附图符号说明：
55.100、获取模块；
56.200、敏感参数确定模块；
57.300、模型建立模块；
58.400、最优模拟参数确定模块；
59.500、敏感参数模拟数据体确定模块；
60.600、结果获取模块；
61.1302、计算机设备；
62.1304、处理器；
63.1306、存储器；
64.1308、驱动机构；
65.1310、输入/输出模块；
66.1312、输入设备；
67.1314、输出设备；
68.1316、呈现设备；
69.1318、图形用户接口；
70.1320、网络接口；
71.1322、通信链路；
72.1324、通信总线。
具体实施方式
73.下面将结合本文实施例中的附图，对本文实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本文一部分实施例，而不是全部的实施例。基于本文中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本文保护的范围。
74.随着我国页岩勘探开发的不断推进，海相页岩气已实现规模开发，但海陆过渡相页岩气仍处在探索之中，我国海陆过渡相页岩气资源丰富，勘探潜力大，是我国页岩气勘探的一个重要领域。但海陆过渡相页岩层多与煤层交互发育，页岩层往往发育于厚煤层之下，页岩层厚度薄、非均质性强，用常规的地震预测方法进行页岩薄储层勘探的过程中，由于煤层的强振幅地震反射特征，致使其下伏的页岩薄储层地震反射很弱甚至为空白地震反射特征，从而干扰页岩薄储层的识别和预测，给此类页岩的勘探和开发造成了很大的困难。因此，用常规的地震预测方法来预测海陆过渡相煤层强屏蔽下页岩薄储层的准确度较低且难度大。
75.为了解决上述问题，本文实施例提供了一种页岩薄储层的预测方法，能够提高页岩薄储层预测的准确度。图1是本文实施例提供的一种页岩薄储层的预测方法的步骤示意图，本说明书提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的劳动可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的系统或装置产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行。
76.参照图1，所述方法可以包括以下步骤：
77.s101：获取工区内每口井对应的钻井数据、多种测井曲线、所述工区内目的层的地
震解释层位数据，以及所述工区内的地震纯波偏移数据。
78.s102：将所述钻井数据和所述多种测井曲线进行综合分析，确定所述工区内目的层段页岩薄储层的敏感参数及其取值范围。
79.s103：利用所述目的层的地震解释层位数据、各层位的接触关系和所述敏感参数对应的测井曲线建立地震波形指示模拟的初始模型。
80.s104：通过所述初始模型、所述地震纯波偏移数据和所述敏感参数对应的测井曲线进行模拟参数测试，确定最优模拟参数。
81.s105：通过所述最优模拟参数进行地震波形指示模拟，确定所述工区内目的层段页岩薄储层的敏感参数模拟数据体。
82.s106：提取所述敏感参数模拟数据体中位于敏感参数的取值范围内的数据，并将其确定为所述页岩薄储层的剖面预测结果。
83.通过将钻井数据和多种测井曲线数据进行综合分析，可以确定工区内页岩薄储层的敏感参数及其取值范围，其中敏感参数可以有效地反映某一储层为页岩薄储层，敏感参数的取值范围可以作为区分页岩薄储层和其他储层的依据，进而能够通过敏感参数及其取值范围将页岩薄储层和其他储层区分开。在地震波形指示模拟的初始模型建立后，通过模拟参数测试，确定最优模拟参数，以最优模拟参数进行地震波形指示模拟，确定工区内的敏感参数模拟数据体。提取敏感参数模拟数据体中位于敏感参数的取值范围内的数据，该数据即为页岩薄储层的剖面预测结果，进而可以得到较为准确的页岩薄储层的预测信息。
84.在本文实施例中，工区内的多口井可以均匀分布在工区内，多口井均匀分布后能够更加广泛的获取工区整个范围内的钻井数据和测井曲线，进而能够提高对页岩薄储层进行预测的准确度。
85.在本文实施例中，每口井对应的钻井数据可包括地质分层数据和岩性数据等，多种测井曲线可包括声波(ac)曲线、自然伽马(gr)曲线、电阻率(rt)曲线、补偿中子(cnl)曲线和密度(den)曲线等。
86.参照图2，在一些实施例中，步骤s103之前，还包括以下步骤：
87.s201：通过合成地震记录井—震标定方法，将所述多种测井曲线从深度域转换为时间域，为地震波形指示模拟的初始模型的建立提供准确的时深关系。
88.在本文实施例中，对每口井来说，利用对应的波阻抗曲线和钻井数据计算出地震反射系数，然后将地震反射系数与子波褶积得到该井的合成地震记录。地震纯波偏移数据可以提供该井处对应的地震反射数据。利用合成地震记录井—震标定方法，将每口井的合成地震记录和地震反射数据进行对比标定，建立深度域地质分层与时间域地震反射层位之间的关系，将深度域的测井曲线转换到时间域，获得准确的时深关系，为地震波形指示模拟的初始模型的建立提供基础。
89.在一些实施例中，所述目的层的地震解释层位数据，包括：
90.通过所述合成地震记录井—震标定方法，确定目的层顶、底的地震解释层位数据；
91.或，通过所述合成地震记录井—震标定方法，确定目的层顶、底的地震解释层位数据以及目的层内部的地震解释层位数据。
92.在本文实施例中，目的层即为判断为页岩薄储层的地质层。通过合成地震记录井—震标定方法，建立深度域的地质分层与时间域的地震反射层位之间的对应关系，可以
确定与深度域上目的层段对应的时间域的地震反射特征。当页岩薄储层较薄，即目的层内部的地震反射特征难以反映出来时，因此只需确定目的层顶、底的地震解释层位数据即可。当页岩薄储层较厚，即目的层内部的地震反射特征可以反映出来时，可以确定目的层顶、底的地震解释层位数据以及目的层内部的地震解释层位数据。
93.参照图3，在一些实施例中，步骤s102包括以下子步骤：
94.s301：将所述钻井数据和所述多种测井曲线进行综合分析，确定页岩薄储层对应的多种测井曲线的特征。
95.s302：根据页岩薄储层对应的多种测井曲线的特征，确定每口井内页岩薄储层的分布。
96.s303：将所述工区内目的层段所有井对应的多种测井曲线进行交会分析，将交会分析的结果与页岩薄储层的分布对应，确定所述工区内页岩薄储层的敏感参数及其取值范围。
97.在本文实施例中，首先，通过钻井数据对已钻页岩薄储层的多种测井曲线进行分析，确定页岩薄储层对应的多种测井曲线的具体特征，例如：自然伽马的高低，声波的大小等。然后，根据与页岩薄储层对应的多种测井曲线的特征，即可分辨每口井内页岩薄储层的具体分布情况。
98.最后，对工区内所有井的多种测井曲线分岩性进行交会分析，交会分析可以为两两交会也可以为三者进行交会，具体的交会方式本文不做限定。根据交会分析的结果，确定最能反映页岩薄储层的测井曲线，该测井曲线为敏感参数对应的测井曲线，其中，敏感参数至少为一个。确定敏感参数之后，进一步确定敏感参数的取值范围，通过该取值范围能够将页岩薄储层与其他储层区分开，在敏感参数的取值范围确定的过程中，对工区范围内的所有井进行页岩薄储层对应敏感参数的取值分析，求取最具代表性的取值范围，可以通过求平均数或求众数的方式来求取最具代表性的取值范围，具体的求范围方式本文不做限定。
99.在本文实施例中，利用目的层的地震解释层位数据、获取到的敏感参数对应的测井曲线，按照各层位的接触关系，其中各层位的接触关系可以为整合接触、平行顶底的接触关系，具体的接触关系根据工区的地质条件确定，本文不做限定，以建立地震波形指示模拟的初始模型。
100.利用初始模型、工区内的地震纯波偏移数据和敏感参数对应的测井曲线进行模拟参数测试，其中模拟参数为平滑系数、截止频率和有效样本数。在进行地震波形指示模拟之前，需要通过测试确定三个参数的最优值，确定最优值的原则是：通过不断调整各参数的取值，优选出能够使模拟的结果与每口井敏感参数对应的测井曲线实测结果误差最小的三个参数的值，为最优模拟参数。
101.其中平滑系数取值范围为2～3，该值越大，其模拟结果的平滑性越好，但误差会增加；有效样本数取值范围为3～6，随着其取值的增大，模拟误差逐渐减小，但达到一定数值后，随着取值的增大误差基本保持不变；截止频率分为低频部分和高频部分，其中低频部分的高通频率范围为5～9hz，高截频率范围为9～12hz；高频部分的低截频率范围为45～60hz，低通频率范围为60～80hz，高通频率范围为150～260hz，高截频率范围为230～350hz。截止频率的取值要根据地震纯波偏移数据的振幅谱进行选择，其取值原则是有效地补偿地震纯波偏移数据所缺失的低频部分和高频部分，在此原则下，再对截止频率值进行
细微调节，找到使模拟结果误差达到最小的参数。
102.通过最优模拟参数进行敏感参数的地震波形指示模拟，确定工区内目的层段页岩薄储层的敏感参数模拟数据体，即在目的层上各个位置处对应的敏感参数的取值。提取敏感参数模拟数据体中位于敏感参数的取值范围内的数据，即可在剖面上将页岩薄层识别出来，可以形成以时间刻度为纵坐标，水平位置为横坐标的剖面图，进而得到页岩薄层的剖面预测结果，该预测结果可以体现页岩薄储层在剖面上的分布情况。
103.参照图4，在本文实施例中，在通过上述步骤得到页岩薄储层的剖面预测结果之后，方法还包括以下步骤：
104.s401：根据所述页岩薄储层的所述剖面预测结果，计算所述工区平面上各点的所述页岩薄储层的累计时间厚度。
105.s402：将所述工区平面上各点的所述页岩薄储层的累计时间厚度进行平面成图，获得页岩薄储层的平面预测结果。
106.参照图5，在一些实施例中，所述计算所述工区平面上各点的所述页岩薄储层的累计时间厚度，包括以下子步骤：
107.s501：将平面上各点处每套所述页岩薄储层的底层和顶层对应的时间值相减，得到该点处每套所述页岩薄储层的时间厚度。
108.s502：将各点对应的每套所述页岩薄储层的时间厚度累加，得到该点处页岩薄储层的累计时间厚度。
109.在本文实施例中，通过页岩薄储层的剖面预测结果进行页岩薄储层的平面预测的过程中，对于工区平面上各点处对应的多套页岩薄储层，首先，将每套页岩薄储层的底层和顶层对应的时间值相减，得到该套页岩薄储层的时间厚度；然后，将各点处对应的多套页岩薄储层的时间厚度累加，得到对应该点出页岩薄储层的累计时间厚度。最后，将工区平面上各点对应的累计时间厚度进行平面成图，得到页岩薄储层的平面预测结果。
110.在一些实施例中，以山西组的山2-3段为目的层进行预测为例：
111.参照图6，在工区内设置有50口井，获取工区内每口井对应的钻井数据、多种测井曲线、工区内目的层的地震解释层位数据，以及工区内的地震纯波偏移数据。通过钻井资料和测井曲线综合分析，识别页岩薄储层的测井曲线特征，认为页岩薄储层具有高自然伽马(gr)、高电阻率(rt)、高声波时差(ac)、高补偿中子(cnl)和低密度(den)的特征，即在页岩薄储层发育层段，一般“自然伽马值明显增高，电阻率增加，声波时差明显变大，补偿中子增大，密度明显降低，且一般伴有井径扩径”的特征。
112.参照图7，通过对密度、声波时差、声波阻抗、自然伽马、电阻率和补偿中子等测井曲线分岩性进行的两两交会分析，最终认为自然伽马和补偿中子的交会图能够更好地反映页岩薄储层，两者结合，能够将页岩薄储层和非储层区分开，因此自然伽马和补偿中子为本实施例中页岩薄储层的敏感参数。并且确定了敏感参数的取值范围为自然伽马gr》120api，补偿中子cnl》19。
113.参照图8，通过50口井的合成地震记录井—震标定，建立深度域的地质分层与时间域的地震反射层位之间的对应关系，识别了目的层顶、底的地震反射特征。目的层顶界(top反射层)为一个单相位、高连续、强振幅的波峰地震反射，顶界这种强振幅高连续的地震反射特征主要由于在目的层顶界发育一套煤层而导致。目的层底界(base反射层)也是一个单
相位、中高连续、中强振幅的波谷地震反射，在目的层底界也发育一套煤层。因此，受目的层顶界煤层强反射的影响，目的层内部的页岩薄层则为弱—空白地震反射，较难识别。
114.参照图9，利用目的层的顶、底地震解释层位和敏感参数对应的测井曲线(自然伽马和补偿中子测井曲线)，按照整合接触、平行顶底的接触关系，建立了地震波形指示模拟的初始模型。并以地震纯波偏移数据、初始模型和敏感参数对应的测井曲线(自然伽马和补偿中子测井曲线)为基础，进行模拟参数测试，参数试验的目的是，确定最优模拟参数，使得模拟结果达到最佳。模拟参数包括平滑系数、截止频率和有效样本数，其中平滑系数取值范围为2～3，该值越大，其模拟结果的平滑性越好，但误差会增加；有效样本数取值范围为3～6个，随着其取值的增大，模拟误差逐渐减小，但达到一定数值后，随着取值的增大误差基本保持不变；截止频率分为低频部分和高频部分，其中低频部分的高通频率范围为5～9hz，高截频率范围为9～12hz；高频部分的低截频率范围为45～60hz，低通频率范围为60～80hz，高通频率范围为150～260hz，高截频率范围为230～350hz。截止频率的取值要根据地震数据的振幅谱进行选择，其取值原则是有效地补偿地震数据所缺失的低频部分和高频部分，在此原则下，再对截止频率值进行细微调节，找到使模拟结果误差达到最小的参数。通过测试，确定最优模拟参数如下：平滑系数为2，有效样本数为5，低频部分的高通频率为7hz，低频部分的高截频率为10hz；高频部分的低截频率为50hz，高频部分的低通频率为70hz，高频部分的高通频率为200hz，高频部分的高截频率为270hz。利用如上模拟参数，进行地震波形指示模拟，分别获得了自然伽马模拟数据体和补偿中子模拟数据体。
115.参照图10，利用地震波形指示模拟得到的自然伽马模拟数据体和补偿中子模拟数据体，根据页岩薄储层对应的敏感参数的取值范围，即自然伽马gr》120api，补偿中子cnl》19，提取两个模拟数据体中对应此取值范围内的数据，即为页岩薄储层的剖面预测结果。
116.参照图11，在目的层山西组山2-3段，将平面上同一位置处所预测的剖面上的各套页岩薄储层的时间厚度进行累加，得到工区内平面上各点的页岩薄储层的累计时间厚度，并对其进行平面成图，得到页岩薄储层的平面预测结果。页岩薄储层在工区的东北部和中部较为发育，a区、b区和c区为有利页岩薄储层发育区，储层平均时间厚度大于5毫秒。
117.基于上述所述的一种页岩薄储层的预测方法，本文实施例还提供一种页岩薄储层的预测装置。所述的装置可以包括使用了本文实施例所述方法的系统(包括分布式系统)、软件(应用)、模块、组件、服务器、客户端等并结合必要的实施硬件的装置。基于同一创新构思，本文实施例提供的一个或多个实施例中的装置如下面的实施例所述。由于装置解决问题的实现方案与方法相似，因此本文实施例具体的装置的实施可以参见前述方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。
118.具体地，图12是本文实施例提供的一种页岩薄储层的预测装置一个实施例的模块结构示意图，参照图12所示，本文实施例提供的一种页岩薄储层的预测装置包括：获取模块100、敏感参数确定模块200、模型建立模块300、最优模拟参数确定模块400、敏感参数模拟数据体确定模块500、结果获取模块600。
119.获取模块100，获取工区内每口井对应的钻井数据、多种测井曲线、所述工区内目的层的地震解释层位数据，以及所述工区内的地震纯波偏移数据；
120.敏感参数确定模块200，将所述钻井数据和所述多种测井曲线进行综合分析，确定所述工区内目的层段页岩薄储层的敏感参数及其取值范围；
121.模型建立模块300，利用所述目的层的地震解释层位数据、各层位的接触关系和所述敏感参数对应的测井曲线建立地震波形指示模拟的初始模型；
122.最优模拟参数确定模块400，通过所述初始模型、所述地震纯波偏移数据和所述敏感参数对应的测井曲线进行模拟参数测试，确定最优模拟参数；
123.敏感参数模拟数据体确定模块500，通过所述最优模拟参数进行地震波形指示模拟，确定所述工区内目的层段页岩薄储层的敏感参数模拟数据体；
124.结果获取模块600，提取所述敏感参数模拟数据体中位于敏感参数的取值范围内的数据，并将其确定为所述页岩薄储层的剖面预测结果。
125.本文一实施例中，参照图13所示，还提供一种计算机设备1302。计算机设备1302可以包括一个或多个处理器1304，诸如一个或多个中央处理单元(cpu)或图形处理器(gpu)，每个处理单元可以实现一个或多个硬件线程。计算机设备1302还可以包括任何存储器1306，其用于存储诸如代码、设置、数据等之类的任何种类的信息，一具体实施方式中，存储器1306上并可在处理器1304上运行的计算机程序，所述计算机程序被所述处理器1304运行时，可以执行根据上述方法的指令。非限制性的，比如，存储器1306可以包括以下任一项或多种组合：任何类型的ram，任何类型的rom，闪存设备，硬盘，光盘等。更一般地，任何存储器都可以使用任何技术来存储信息。进一步地，任何存储器可以提供信息的易失性或非易失性保留。进一步地，任何存储器可以表示计算机设备1302的固定或可移除部件。在一种情况下，当处理器1304执行被存储在任何存储器或存储器的组合中的相关联的指令时，计算机设备1302可以执行相关联指令的任一操作。计算机设备1302还包括用于与任何存储器交互的一个或多个驱动机构1308，诸如硬盘驱动机构、光盘驱动机构等。
126.计算机设备1302还可以包括输入/输出模块1310(i/o)，其用于接收各种输入(经由输入设备1312)和用于提供各种输出(经由输出设备1314)。一个具体输出机构可以包括呈现设备1316和相关联的图形用户接口1318(gui)。在其他实施例中，还可以不包括输入/输出模块1310(i/o)、输入设备1312以及输出设备1314，仅作为网络中的一台计算机设备。计算机设备402还可以包括一个或多个网络接口1320，其用于经由一个或多个通信链路1322与其他设备交换数据。一个或多个通信总线1324将上文所描述的部件耦合在一起。
127.通信链路1322可以以任何方式实现，例如，通过局域网、广域网(例如，因特网)、点对点连接等、或其任何组合。通信链路1322可以包括由任何协议或协议组合支配的硬连线链路、无线链路、路由器、网关功能、名称服务器等的任何组合。
128.对应于图1-图5中的方法，本文实施例还提供了一种计算机可读存储介质，该计算机可读存储介质上存储有计算机程序，该计算机程序被处理器运行时执行上述方法的步骤。
129.本文实施例还提供一种计算机可读指令，其中当处理器执行所述指令时，其中的程序使得处理器执行如图1至图5所示的方法。
130.应理解，在本文的各种实施例中，上述各过程的序号的大小并不意味着执行顺序的先后，各过程的执行顺序应以其功能和内在逻辑确定，而不应对本文实施例的实施过程构成任何限定。
131.还应理解，在本文实施例中，术语“和/或”仅仅是一种描述关联对象的关联关系，表示可以存在三种关系。例如，a和/或b，可以表示：单独存在a，同时存在a和b，单独存在b这三种情况。另外，本文中字符“/”，一般表示前后关联对象是一种“或”的关系。
132.本领域普通技术人员可以意识到，结合本文中所公开的实施例描述的各示例的单元及算法步骤，能够以电子硬件、计算机软件或者二者的结合来实现，为了清楚地说明硬件和软件的可互换性，在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本文的范围。
133.所属领域的技术人员可以清楚地了解到，为了描述的方便和简洁，上述描述的系统、装置和单元的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
134.在本文所提供的几个实施例中，应该理解到，所揭露的系统、装置和方法，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另外，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口、装置或单元的间接耦合或通信连接，也可以是电的，机械的或其它的形式连接。
135.所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本文实施例方案的目的。
136.另外，在本文各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以是两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。
137.所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本文的技术方案本质上或者说对现有技术做出贡献的部分，或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本文各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：u盘、移动硬盘、只读存储器(rom，read-only memory)、随机存取存储器(ram，random access memory)、磁碟或者光盘等各种可以存储程序代码的介质。
138.本文中应用了具体实施例对本文的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本文的方法及其核心思想；同时，对于本领域的一般技术人员，依据本文的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本文的限制。