具体实施方式

为使本申请实施例的目的、技术方案和优点更加清楚，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整的描述，显然，所描述的实施例是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。

为了有效提高地质统计学反演横向的控制性，并回避虚拟井点密集而产生的大量运算，本申请提供一种基于相约束可变网格的地质统计学反演方法，也即：由地震相确定的变网格约束虚拟井确定方法；同时，在数据准备阶段，采用新的切比雪夫近似公式来描述正态分布的反函数，以提高正态变换的效率。

基于上述内容，本申请提供一种基于相约束可变网格的地质统计学反演方法的实施例，参见图1，所述基于相约束可变网格的地质统计学反演方法具体包含有如下内容：

步骤100：根据目标区域的地震相数据，生成所述目标区域对应的可变网格，并自该可变网格中确定多个目标点并由预设数据抽取多个虚拟井，其中，所述预设数据包括：确定性反演结果。

可以理解的是，所述目标区域是指地震勘探的作业区域，所述预设数据并不仅限于确定性反演结果，还可以包含有其它可用数据，具体可以根据实际应用情形设置。

在步骤100中，所述确定性反演结果是根据预设的确定性地震反演方式获取的，所述确定性反演方式是指根据目标区域中的已经确定的实钻井进行的地震反演(seismicinversion)。其中，所述地震反演是指把常规的界面型反射剖面转换成岩层型的测井剖面，将地震资料变成可与测井资料直接对比的形式，实现这种转换的处理过程。地震反演是利用地表观测地震资料，以已知地质规律和钻井、测井资料为约束，对地下岩层空间结构和物理性质进行成像(求解)的过程，广义的地震反演包含了地震处理解释的整个内容。

步骤200：对所述目标区域中待进行高斯模拟处理的实钻井和虚拟井进行正态变换，得到对应的正态变换结果数据。

步骤300：确定各个所述目标点预设距离范围内的已知属性值的均值和方差，并基于该均值和方差构建高斯分布，得到该高斯分布的累积概率分布。

步骤400：基于所述累积概率分布模拟预设距离范围内未知点储层参数的实现，将该未知点储层参数的实现与地震数据最佳匹配的实现作为地质统计学反演结果。

在步骤400中，可以输出用于获取所述目标区域的储层数据的地质统计学反演结果值用户的客户端设备或显示设备上。

从上述描述可知，本申请实施例提供的基于相约束可变网格的地质统计学反演方法，通过根据目标区域的地震相数据，生成所述目标区域对应的可变网格，并自该可变网格中确定多个目标点并由预设数据抽取多个虚拟井，其中，所述预设数据包括：确定性反演结果；对所述目标区域中待进行高斯模拟处理的实钻井和虚拟井进行正态变换，得到对应的正态变换结果数据；确定各个所述目标点预设距离范围内的已知属性值的均值和方差，并基于该均值和方差构建高斯分布，得到该高斯分布的累积概率分布；基于所述累积概率分布模拟预设距离范围内未知点储层参数的实现，将该未知点储层参数的实现与地震数据最佳匹配的实现作为地质统计学反演结果，能够将正交等距网格横向控制约束改进为由相带控制的可变网格约束，能够有效提高地质统计学反演横向的控制性，并能够有效降低由于虚拟井点密集而导致的数据运算量，同时改进了输入储层参数正态变换实施方式，进而能够有效提高地质统计学反演的效率及准确性，以有效提高基于地质统计学反演结果获取储层数据的准确性及效率，进而能够有效且快速的进行储层预测，并有效节省人力、物力成本。

在一种基于相约束可变网格的地质统计学反演方法的实施例中，参见图2，所述基于相约束可变网格的地质统计学反演方法中的步骤100具体包含有如下内容：

步骤110：获取所述目标区域的地震相数据。

步骤120：对目标区域实施预设的实施高分辨率确定性地震反演，得到所述目标区域的确定性地震反演结果。

步骤130：根据目标区域的地震相数据，生成所述目标区域对应的可变网格，其中，所述可变网格的网格大小根据道间距和网格阶数确定。

步骤140：自所述可变网格中确定多个目标点并由所述确定性反演结果抽取多个虚拟井。

具体来说，地震相可以由地质特征和地球物理如利用地震属性来描述，通常划分为三到五种，通过地震相来确定不同网格的大小。可变网格的具体实现可以描述为，在整个工作区域由下式确定初始网格和加密网格，即：

g＝dx×2^r

其中，g为网格大小，dx为道间距，r为阶数，取值为1，2，3，4和5。初始网格可以通过较大的阶数如r＝4或者r＝5来实现。在有利相带确定的有利储层区域，阶数r＝1或者r＝2。

在一种基于相约束可变网格的地质统计学反演方法的实施例中，所述基于相约束可变网格的地质统计学反演方法中的步骤200还可以具体包含有如下内容：

步骤210：应用正态分布反函数的预设近似函数，对所述目标区域中待进行高斯模拟处理的实钻井和虚拟井进行正态变换，得到对应的正态变换结果数据。

其中，所述预设近似函数ftf(t)如公式一所示：

在公式一中，t＝1.0–ru，且ru为(0,1)间的均匀随机数；a₁＝2.505922；a₃＝-15.73223；a₅＝23.54337；b₂＝-7.337743；b₄＝14.97266；b₆＝-6.016088。

具体来说，油气开发阶段常用地质统计学反演方法来精细描述油气储层，测井数据是地质统计学反演方法高分辨信息的来源。测井数据通常不满足正态分布特征，需要对测井数据进行正态变换。常规正正态变换采用正态分布反函数的Hastings近似：

这里，ru是(0,1)均匀分布的随机数，如果ru大于0.5，则ru等于1.0减去原ru值；其他参数a₀＝2.155517，a₁＝0.802853，a₂＝0.010328，b₁＝1.432788，b₂＝0.189269，b₃＝0.001308。

而在本申请的一个实施例中，则是采用正态分布反函数的Byars-Roscoe近似来实施正态变换：

这里，t＝1.0–ru，且ru为(0,1)间的均匀随机数。其他参数a₁＝2.505922，a₃＝-15.73223，a₅＝23.54337，b₂＝-7.337743，b₄＝14.97266，b₆＝-6.016088。

在一种基于相约束可变网格的地质统计学反演方法的实施例中，参见图3，所述基于相约束可变网格的地质统计学反演方法中的步骤300具体包含有如下内容：

步骤310：应用预设的克里金法确定各个所述目标点预设距离范围内的已知属性值的均值和方差。

步骤320：根据各个所述目标点预设距离范围内的已知属性值的均值和方差构建高斯分布，得到该高斯分布的累积概率分布。

在一种基于相约束可变网格的地质统计学反演方法的实施例中，所述基于相约束可变网格的地质统计学反演方法中的步骤400的第一种实现方式具体包含有如下内容：

步骤411：在均匀分布的随机数(0，1)之间进行一次随机抽样，得到对应的抽样结果。

在一种基于相约束可变网格的地质统计学反演方法的实施例中，参见图4，所述基于相约束可变网格的地质统计学反演方法中的步骤400的第二种实现方式具体包含有如下内容：

步骤421：在均匀分布的随机数(0，1)之间进行多次随机抽样，得到对应的多个抽样结果。

步骤422：基于预设的褶积理论数据在每个所述目标点各自对应的多个储层参数实现数据中选取与所述目标区域对应的地震数据的匹配度最高的一个实现，作为各个所述目标点唯一对应的储层参数。

步骤423：将各个所述目标点各自对应的所述储层参数作为用于获取所述目标区域的地质统计学反演结果。

综上所述，本申请实现了一种基于相约束可变网格的地质统计学反演方法，能够有效提高地质统计学反演横向的控制性，并能够有效降低由于虚拟井点密集而导致的数据运算量，同时改进了输入储层参数正态变换实施方式，进而能够有效提高地质统计学反演的效率及准确性，以有效提高基于地质统计学反演结果获取储层数据的准确性及效率，进而能够有效且快速的进行储层预测，并有效节省人力、物力成本。

从软件层面来说，为了有效提高地质统计学反演横向的控制性，并回避虚拟井点密集而产生的大量运算，本申请提供一种用于实现所述基于相约束可变网格的地质统计学反演方法中全部或部分内容的基于相约束可变网格的地质统计学反演装置的实施例，参见图5，所述基于相约束可变网格的地质统计学反演装置具体包含有如下内容：

可变网格生成模块10，用于根据目标区域的地震相数据，生成所述目标区域对应的可变网格，并自该可变网格中确定多个目标点并由预设数据抽取多个虚拟井，其中，所述预设数据包括：确定性反演结果。

正态变换模块20，用于对所述目标区域中待进行高斯模拟处理的实钻井和虚拟井进行正态变换，得到对应的正态变换结果数据。

累积概率分布模块30，用于确定各个所述目标点预设距离范围内的已知属性值的均值和方差，并基于该均值和方差构建高斯分布，得到该高斯分布的累积概率分布。

反演结果获取模块40，用于基于所述累积概率分布模拟预设距离范围内未知点储层参数的实现，将该未知点储层参数的实现与地震数据最佳匹配的实现作为地质统计学反演结果。

本申请提供的基于相约束可变网格的地质统计学反演装置的实施例具体可以用于执行上述实施例中的基于相约束可变网格的地质统计学反演方法的实施例的处理流程，其功能在此不再赘述，可以参照上述方法实施例的详细描述。

从上述描述可知，本申请实施例提供的基于相约束可变网格的地质统计学反演装置，通过根据目标区域的地震相数据，生成所述目标区域对应的可变网格，并自该可变网格中确定多个目标点并由预设数据抽取多个虚拟井，其中，所述预设数据包括：确定性反演结果；对所述目标区域中待进行高斯模拟处理的实钻井和虚拟井进行正态变换，得到对应的正态变换结果数据；确定各个所述目标点预设距离范围内的已知属性值的均值和方差，并基于该均值和方差构建高斯分布，得到该高斯分布的累积概率分布；基于所述累积概率分布模拟预设距离范围内未知点储层参数的实现，将该未知点储层参数的实现与地震数据最佳匹配的实现作为地质统计学反演结果，能够将正交等距网格横向控制约束改进为由相带控制的可变网格约束，能够有效提高地质统计学反演横向的控制性，并能够有效降低由于虚拟井点密集而导致的数据运算量，同时改进了输入储层参数正态变换实施方式，进而能够有效提高地质统计学反演的效率及准确性，以有效提高基于地质统计学反演结果获取储层数据的准确性及效率，进而能够有效且快速的进行储层预测，并有效节省人力、物力成本。

在一种基于相约束可变网格的地质统计学反演装置的实施例中，所述基于相约束可变网格的地质统计学反演装置中的可变网格生成模块10用于执行下述内容：

步骤110：获取所述目标区域的地震相数据。

步骤120：对目标区域实施预设的实施高分辨率确定性地震反演，得到所述目标区域的确定性地震反演结果。

步骤130：根据目标区域的地震相数据，生成所述目标区域对应的可变网格，其中，所述可变网格的网格大小根据道间距和网格阶数确定。

步骤140：自所述可变网格中确定多个目标点并由所述确定性反演结果抽取多个虚拟井。

具体来说，地震相可以由地质特征和地球物理如利用地震属性来描述，通常划分为三到五种，通过地震相来确定不同网格的大小。发明的可变网格的具体实现可以描述为，在整个工作区域由下式确定初始网格和加密网格，即：

g＝dx×2^r

其中，g为网格大小，dx为道间距，r为阶数，取值为1，2，3，4和5。初始网格可以通过较大的阶数如r＝4或者r＝5来实现。在有利相带确定的有利储层区域，阶数r＝1或者r＝2。

在一种基于相约束可变网格的地质统计学反演装置的实施例中，所述基于相约束可变网格的地质统计学反演装置中正态变换模块20用于执行下述内容：

步骤210：应用正态分布反函数的预设近似函数，对所述目标区域中待进行高斯模拟处理的实钻井和虚拟井进行正态变换，得到对应的正态变换结果数据。

其中，所述预设近似函数ftf(t)如公式一所示：

在公式一中，t＝1.0–ru，且ru为(0,1)间的均匀随机数；a₁＝2.505922；a₃＝-15.73223；a₅＝23.54337；b₂＝-7.337743；b₄＝14.97266；b₆＝-6.016088。

具体来说，油气开发阶段常用地质统计学反演方法来精细描述油气储层，测井数据是地质统计学反演方法高分辨信息的来源。测井数据通常不满足正态分布特征，需要对测井数据进行正态变换。常规正正态变换采用正态分布反函数的Hastings近似：

这里，ru是(0,1)均匀分布的随机数，如果ru大于0.5，则ru等于1.0减去原ru值；其他参数a₀＝2.155517，a₁＝0.802853，a₂＝0.010328，b₁＝1.432788，b₂＝0.189269，b₃＝0.001308。

而在本申请的一个实施例中，则是采用正态分布反函数的Byars-Roscoe近似来实施正态变换：

这里，t＝1.0–ru，且ru为(0,1)间的均匀随机数。其他参数a₁＝2.505922，a₃＝-15.73223，a₅＝23.54337，b₂＝-7.337743，b₄＝14.97266，b₆＝-6.016088。

在一种基于相约束可变网格的地质统计学反演装置的实施例中，所述基于相约束可变网格的地质统计学反演装置中的累积概率分布模块30具体用于包含有如下内容：

步骤310：应用预设的克里金法确定各个所述目标点预设距离范围内的已知属性值的均值和方差。

步骤320：根据各个所述目标点预设距离范围内的已知属性值的均值和方差构建高斯分布，得到该高斯分布的累积概率分布。

在一种基于相约束可变网格的地质统计学反演装置的实施例中，所述基于相约束可变网格的地质统计学反演装置中的反演结果获取模块40的具体用于包含有如下内容：

步骤411：在均匀分布的随机数(0，1)之间进行一次随机抽样，得到对应的抽样结果。

在一种基于相约束可变网格的地质统计学反演装置的实施例中，所述基于相约束可变网格的地质统计学反演装置中的反演结果获取模块40的还具体用于包含有如下内容：

步骤421：在均匀分布的随机数(0，1)之间进行多次随机抽样，得到对应的多个抽样结果。

步骤422：基于预设的褶积理论数据在每个所述目标点各自对应的多个储层参数实现数据中选取与所述目标区域对应的地震数据的匹配度最高的一个实现，作为各个所述目标点唯一对应的储层参数。

步骤423：将各个所述目标点各自对应的所述储层参数作为用于获取所述目标区域的地质统计学反演结果。

综上所述，本申请提供了一种基于相约束可变网格的地质统计学反演装置，能够有效提高地质统计学反演横向的控制性，并能够有效降低由于虚拟井点密集而导致的数据运算量，同时改进了输入储层参数正态变换实施方式，进而能够有效提高地质统计学反演的效率及准确性，以有效提高基于地质统计学反演结果获取储层数据的准确性及效率，进而能够有效且快速的进行储层预测，并有效节省人力、物力成本。

为进一步说明本方案，本申请还提供一种应用基于相约束可变网格的地质统计学反演装置实现基于相约束可变网格的地质统计学反演方法的具体应用实例，其中，地震相可以由地质特征和地球物理如利用地震属性来描述，通常划分为三到五种，通过地震相来确定不同网格的大小。发明的可变网格的具体实现可以描述为，在整个工作区域由下式确定初始网格和加密网格，即：

g＝dx×2^r

其中，g为网格大小，dx为道间距，r为阶数，取值为1，2，3，4和5。初始网格可以通过较大的阶数如r＝4或者r＝5来实现。在有利相带确定的有利储层区域，阶数r＝1或者r＝2。

油气开发阶段常用地质统计学反演方法来精细描述油气储层，测井数据是地质统计学反演方法高分辨信息的来源。测井数据通常不满足正态分布特征，需要对测井数据进行正态变换。常规正正态变换采用正态分布反函数的Hastings近似：

这里，ru是(0,1)均匀分布的随机数，如果ru大于0.5，则ru等于1.0减去原ru值；其他参数a₀＝2.155517，a₁＝0.802853，a₂＝0.010328，b₁＝1.432788，b₂＝0.189269，b₃＝0.001308。

本应用实例采用正态分布反函数的Byars-Roscoe近似来实施正态变换：

这里，t＝1.0–ru，且ru为(0,1)间的均匀随机数。其他参数a₁＝2.505922，a₃＝-15.73223，a₅＝23.54337，b₂＝-7.337743，b₄＝14.97266，b₆＝-6.016088。

在一种具体应用举例中，参见图6，所述基于相约束可变网格的地质统计学反演方法具体可以包含有如下内容：

1.叠研究区域实施高分辨率确定性地震反演；

2.根据沉积相及确定性反演结果，采用可变网格抽出虚拟井；

3.参与序贯高斯模拟的实钻井及虚拟井属性值按照Byars-Roscoe近似进行正态变换；

4.由克里金算法计算待估计点一定范围内的已知属性值的均值和方差，并由此均值和方差构建高斯分布，得到高斯分布的累积概率分布；

5.在均匀分布的随机数(0，1)间随机抽样，由高斯分布的累积概率分布得到未知点储层参数的一个实现；

6.在均匀分布的随机数(0，1)间多次随机抽样，由高斯分布的累积概率分布未知点储层参数的多个实现，由褶积理论判断这些实现，输出与地震数据最佳匹配的实现，作为地质统计学反演的结果。

参见图7所示的基于地震相约束的可变网格虚拟井抽取结果示意图，由相带确定的曲线圈定区域，是储层发育有利区，可采用本应用实例所确定的调节阶数方法使得虚拟井适当加密。

参见图8所示的正态变换前的储层属性和图9所示的正态变换前的储层属性分布直方图(双峰，非正态分布)，在执行正态变换操作后，正态变换后的储层属性参见图10，正态变换之后的储层属性分布直方图(单峰，服从正态分布特征)参见图11。

基于上述内容，常规地质统计学反演方法反演效果参见图12，本申请提供的基于相约束可变网格的地质统计学反演方法反演效果参见图13，由此可知，基于相约束可变网格的地质统计学反演方法结果更符合地质规律，更有利于储层的精细描述。

本申请应用实例提供的基于相约束可变网格的地质统计学反演方法，较好地解决了地质统计学反演纵横向的井间控制问题，使得反演结果更合理，降低地下地质情况的多解性，利于油气储层精细描述；同时，本应用实例将地质统计学反演必要步骤中的正态变换由Hastings近似改变为Byars-Roscoe近似后，提高了正态变换的效率。

从硬件层面来说，本申请提供一种用于实现所述基于相约束可变网格的地质统计学反演方法中的全部或部分内容的电子设备的实施例，所述电子设备具体包含有如下内容：

图14为本申请实施例的电子设备9600的装置构成的示意框图。如图14所示，该电子设备9600可以包括中央处理器9100和存储器9140；存储器9140耦合到中央处理器9100。值得注意的是，该图14是示例性的；还可以使用其他类型的结构，来补充或代替该结构，以实现电信功能或其他功能。

在一实施例中，基于相约束可变网格的地质统计学反演功能可以被集成到中央处理器中。其中，中央处理器可以被配置为进行如下控制：

步骤100：根据目标区域的地震相数据，生成所述目标区域对应的可变网格，并自该可变网格中确定多个目标点并由预设数据抽取多个虚拟井，其中，所述预设数据包括：确定性反演结果。

可以理解的是，所述目标区域是指地震勘探的作业区域，所述预设数据并不仅限于确定性反演结果，还可以包含有其它可用数据，具体可以根据实际应用情形设置。

在步骤100中，所述确定性反演结果是根据预设的确定性地震反演方式获取的，所述确定性反演方式是指根据目标区域中的已经确定的实钻井进行的地震反演(seismicinversion)。其中，所述地震反演是指把常规的界面型反射剖面转换成岩层型的测井剖面，将地震资料变成可与测井资料直接对比的形式，实现这种转换的处理过程。地震反演是利用地表观测地震资料，以已知地质规律和钻井、测井资料为约束，对地下岩层空间结构和物理性质进行成像(求解)的过程，广义的地震反演包含了地震处理解释的整个内容。

步骤200：对所述目标区域中待进行高斯模拟处理的实钻井和虚拟井进行正态变换，得到对应的正态变换结果数据。

步骤300：确定各个所述目标点预设距离范围内的已知属性值的均值和方差，并基于该均值和方差构建高斯分布，得到该高斯分布的累积概率分布。

步骤400：基于所述累积概率分布模拟预设距离范围内未知点储层参数的实现，将该未知点储层参数的实现与地震数据最佳匹配的实现作为地质统计学反演结果。

在步骤400中，可以输出用于获取所述目标区域的储层数据的地质统计学反演结果值用户的客户端设备或显示设备上。

从上述描述可知，本申请实施例提供的电子设备，通过根据目标区域的地震相数据，生成所述目标区域对应的可变网格，并自该可变网格中确定多个目标点并由预设数据抽取多个虚拟井，其中，所述预设数据包括：确定性反演结果；对所述目标区域中待进行高斯模拟处理的实钻井和虚拟井进行正态变换，得到对应的正态变换结果数据；确定各个所述目标点预设距离范围内的已知属性值的均值和方差，并基于该均值和方差构建高斯分布，得到该高斯分布的累积概率分布；基于所述累积概率分布模拟预设距离范围内未知点储层参数的实现，将该未知点储层参数的实现与地震数据最佳匹配的实现作为地质统计学反演结果，能够将正交等距网格横向控制约束改进为由相带控制的可变网格约束，能够有效提高地质统计学反演横向的控制性，并能够有效降低由于虚拟井点密集而导致的数据运算量，同时改进了输入储层参数正态变换实施方式，进而能够有效提高地质统计学反演的效率及准确性，以有效提高基于地质统计学反演结果获取储层数据的准确性及效率，进而能够有效且快速的进行储层预测，并有效节省人力、物力成本。

在另一个实施方式中，基于相约束可变网格的地质统计学反演装置可以与中央处理器9100分开配置，例如可以将基于相约束可变网格的地质统计学反演装置配置为与中央处理器9100连接的芯片，通过中央处理器的控制来实现基于相约束可变网格的地质统计学反演功能。

如图14所示，该电子设备9600还可以包括：通信模块9110、输入单元9120、音频处理器9130、显示器9160、电源9170。值得注意的是，电子设备9600也并不是必须要包括图14中所示的所有部件；此外，电子设备9600还可以包括图14中没有示出的部件，可以参考现有技术。

如图14所示，中央处理器9100有时也称为控制器或操作控件，可以包括微处理器或其他处理器装置和/或逻辑装置，该中央处理器9100接收输入并控制电子设备9600的各个部件的操作。

其中，存储器9140，例如可以是缓存器、闪存、硬驱、可移动介质、易失性存储器、非易失性存储器或其它合适装置中的一种或更多种。可储存上述与失败有关的信息，此外还可存储执行有关信息的程序。并且中央处理器9100可执行该存储器9140存储的该程序，以实现信息存储或处理等。

输入单元9120向中央处理器9100提供输入。该输入单元9120例如为按键或触摸输入装置。电源9170用于向电子设备9600提供电力。显示器9160用于进行图像和文字等显示对象的显示。该显示器例如可为LCD显示器，但并不限于此。

该存储器9140可以是固态存储器，例如，只读存储器(ROM)、随机存取存储器(RAM)、SIM卡等。还可以是这样的存储器，其即使在断电时也保存信息，可被选择性地擦除且设有更多数据，该存储器的示例有时被称为EPROM等。存储器9140还可以是某种其它类型的装置。存储器9140包括缓冲存储器9141(有时被称为缓冲器)。存储器9140可以包括应用/功能存储部9142，该应用/功能存储部9142用于存储应用程序和功能程序或用于通过中央处理器9100执行电子设备9600的操作的流程。

存储器9140还可以包括数据存储部9143，该数据存储部9143用于存储数据，例如联系人、数字数据、图片、声音和/或任何其他由电子设备使用的数据。存储器9140的驱动程序存储部9144可以包括电子设备的用于通信功能和/或用于执行电子设备的其他功能(如消息传送应用、通讯录应用等)的各种驱动程序。

通信模块9110即为经由天线9111发送和接收信号的发送机/接收机9110。通信模块(发送机/接收机)9110耦合到中央处理器9100，以提供输入信号和接收输出信号，这可以和常规移动通信终端的情况相同。

基于不同的通信技术，在同一电子设备中，可以设置有多个通信模块9110，如蜂窝网络模块、蓝牙模块和/或无线局域网模块等。通信模块(发送机/接收机)9110还经由音频处理器9130耦合到扬声器9131和麦克风9132，以经由扬声器9131提供音频输出，并接收来自麦克风9132的音频输入，从而实现通常的电信功能。音频处理器9130可以包括任何合适的缓冲器、解码器、放大器等。另外，音频处理器9130还耦合到中央处理器9100，从而使得可以通过麦克风9132能够在本机上录音，且使得可以通过扬声器9131来播放本机上存储的声音。

本申请的实施例还提供能够实现上述实施例中的基于相约束可变网格的地质统计学反演方法中全部步骤的一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序，该计算机程序被处理器执行时实现上述实施例中的执行主体为服务器或客户端的基于相约束可变网格的地质统计学反演方法的全部步骤，例如，所述处理器执行所述计算机程序时实现下述步骤：

步骤100：根据目标区域的地震相数据，生成所述目标区域对应的可变网格，并自该可变网格中确定多个目标点并由预设数据抽取多个虚拟井，其中，所述预设数据包括：确定性反演结果。

可以理解的是，所述目标区域是指地震勘探的作业区域，所述预设数据并不仅限于确定性反演结果，还可以包含有其它可用数据，具体可以根据实际应用情形设置。

在步骤100中，所述确定性反演结果是根据预设的确定性地震反演方式获取的，所述确定性反演方式是指根据目标区域中的已经确定的实钻井进行的地震反演(seismicinversion)。其中，所述地震反演是指把常规的界面型反射剖面转换成岩层型的测井剖面，将地震资料变成可与测井资料直接对比的形式，实现这种转换的处理过程。地震反演是利用地表观测地震资料，以已知地质规律和钻井、测井资料为约束，对地下岩层空间结构和物理性质进行成像(求解)的过程，广义的地震反演包含了地震处理解释的整个内容。

步骤200：对所述目标区域中待进行高斯模拟处理的实钻井和虚拟井进行正态变换，得到对应的正态变换结果数据。

步骤300：确定各个所述目标点预设距离范围内的已知属性值的均值和方差，并基于该均值和方差构建高斯分布，得到该高斯分布的累积概率分布。

步骤400：基于所述累积概率分布模拟预设距离范围内未知点储层参数的实现，将该未知点储层参数的实现与地震数据最佳匹配的实现作为地质统计学反演结果。

在步骤400中，可以输出用于获取所述目标区域的储层数据的地质统计学反演结果值用户的客户端设备或显示设备上。

从上述描述可知，本申请实施例提供的计算机可读存储介质，通过根据目标区域的地震相数据，生成所述目标区域对应的可变网格，并自该可变网格中确定多个目标点并由预设数据抽取多个虚拟井，其中，所述预设数据包括：确定性反演结果；对所述目标区域中待进行高斯模拟处理的实钻井和虚拟井进行正态变换，得到对应的正态变换结果数据；确定各个所述目标点预设距离范围内的已知属性值的均值和方差，并基于该均值和方差构建高斯分布，得到该高斯分布的累积概率分布；基于所述累积概率分布模拟预设距离范围内未知点储层参数的实现，将该未知点储层参数的实现与地震数据最佳匹配的实现作为地质统计学反演结果，能够将正交等距网格横向控制约束改进为由相带控制的可变网格约束，能够有效提高地质统计学反演横向的控制性，并能够有效降低由于虚拟井点密集而导致的数据运算量，同时改进了输入储层参数正态变换实施方式，进而能够有效提高地质统计学反演的效率及准确性，以有效提高基于地质统计学反演结果获取储层数据的准确性及效率，进而能够有效且快速的进行储层预测，并有效节省人力、物力成本。

本领域内的技术人员应明白，本发明的实施例可提供为方法、装置、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。

本发明是参照根据本发明实施例的方法、设备(装置)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。

本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。