[发明专利]一种天然气水合物饱和度计算方法及系统

说明书

本发明涉及一种天然气水合物饱和度计算方法及系统，该方法包括：基于已分类的样本数据及贝叶斯算法建立判别函数，并通过所述判别函数对样本数据和待判别测井数据进行分类；通过决策超平面在多维空间中对不同组样本的划分，建立天然气水合物饱和度与判别函数的关系，以计算天然气水合物饱和度。通过该方案可以综合有效利用各测井曲线，提高天然气水合物饱和度计算的准确性。

技术领域

本发明涉及油气勘探领域，尤其涉及一种天然气水合物饱和度计算方法及系统。

背景技术

在天然气勘探过程中，需要预先探明天然气位置及储量，对于天然气水合物开采，由于开采难度更大，更需要提前探测其位置及存储量。由于天然气水合物的赋存机理与常规油气不同，主要是通过低渗透性岩层来抑制孔隙中的烃类气体逸散，因此天然气水合物储层多为粉砂岩、泥岩和油页岩等低孔低渗储层，这导致天然气水合物的测井识别与储层参数计算难度加大。

对于天然气水合物的储量，可以通过储层的孔隙度和饱和度来衡量，一般孔隙度计算相对容易，而天然气水合物的饱和度计算通常是重难点。在天然气水合物的饱和度计算模型中，一种是将水合物看作岩石中的矿物颗粒成分，对于温压条件已经发生变化的冻土深度以下的地层而言，天然气水合物会由固态直接气化，而将水合物看作矿物颗粒成分难以表征实际测井环境；另一种是将天然气水合物当作孔隙中的流体，从而方便引入油气解释中的饱和度计算模型，这类饱和度计算方法通常有电阻率法和声波时差法。电阻率法有阿尔奇公式、修正阿尔奇公式、双水模型、印度尼西亚公式等，声波时差法有时间平均方程、修正伍德方程、等效介质理论、热弹性理论、BGTL理论、K-T方程等，对于每一种饱和度算法都有其适用范围，不同算法的饱和度计算结果相差非常大。这些常规饱和度计算模型通常只用到某一项关键测井参数，而基于石油或天然气的单一测井曲线计算饱和度的方法难以全面表征水合物的实际赋存情况，影响饱和度计算结果的准确性。

发明内容

有鉴于此，本发明实施例提供了一种天然气水合物饱和度计算方法及系统，以解决现有天然气水合物饱和度计算准确性不高的问题。

在本发明实施例的第一方面，提供了一种天然气水合物饱和度计算方法，包括：

基于已分类的样本数据及贝叶斯算法建立判别函数，并通过所述判别函数对样本数据和待判别测井数据进行分类；

通过决策超平面在多维空间中对不同组样本的划分，建立天然气水合物饱和度与判别函数的关系，以计算天然气水合物饱和度。

在本发明实施例的第二方面，提供了一种天然气水合物饱和度计算系统，包括：

分类模块，用于基于已分类的样本数据及贝叶斯算法建立判别函数，并通过所述判别函数对样本数据和待判别测井数据进行分类；

计算模块，用于通过决策超平面在多维空间中对不同组样本的划分，建立天然气水合物饱和度与判别函数的关系，以计算天然气水合物饱和度。

在本发明实施例中，基于已分类的样本数据及贝叶斯算法建立判别函数，并通过所述判别函数对样本数据和待判别测井数据进行分类；通过决策超平面在多维空间中对不同组样本的划分，建立天然气水合物饱和度与判别函数的关系，以计算天然气水合物饱和度。可以有效综合利用各测井曲线信息，提高计算精度，解决了现有天然气水合物饱和度计算准确度不高的问题。同时，该方案实现方式简单，计算速度快，可以及时提供关键气层位和取芯测试深度的建议。

附图说明

为了更清楚地说明本发明实施例中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单介绍，显而易见地，下面描述的附图仅仅是本发明的一些实施例，对本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获取其他附图。

图1为本发明的一个实施例提供的天然气水合物饱和度计算方法的流程示意图；