基于无阻流量的含气饱和度校正方法和装置

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN112630853A(43)申请公布日2021.04.09(21)申请号201910954997.7(22)申请日2019.10.09(71)申请人中国石油化工股份有限公司地址100728北京市朝阳区朝阳门北大街22号申请人中国石油化工股份有限公司石油物探技术研究院(72)发明人张克非　李呈呈　林正良　胡华锋　(74)专利代理机构北京思创毕升专利事务所11218代理人孙向民　廉莉莉(51)Int.Cl.G01V9/00(2006.01)权利要求书2页说明书5页附图3页(54)发明名称基于无阻流量的含气饱和度校正方法和装置(57)摘要本申请公开了基于无阻流量的含气饱和度校正方法和装置。该方法包括：基于地质信息和 工程 路基工程安全技术交底工程项目施工成本控制工程量增项单年度零星工程技术标正投影法基本原理 先验信息确定工区内无阻流量与含气饱和度之间的关系；确定无阻流量和含气饱和度估计值之间的关系曲线；根据无阻流量计算每口井的含气饱和度估计值；根据测井数据得到每口井的原始含气饱和度平均值；基于含气饱和度估计值和原始含气饱和度平均值，得到每口井的含气饱和度校正系数；采用含气饱和度校正系数对每口井的原始含气饱和度进行校正。根据本申请，可以获得与开发生产资料及地质认识相符的含气饱和度曲线，能有效指导含气性预测等后续一系列工作，在气藏发育地区具有极大的推广应用价值。CN112630853ACN112630853A权　利　要　求　书1/2页1.一种基于无阻流量的含气饱和度校正方法，其特征在于，所述方法包括：基于地质信息和工程先验信息确定工区内无阻流量与含气饱和度之间的关系；根据无阻流量与含气饱和度之间的关系，确定无阻流量和含气饱和度估计值之间的关系曲线；根据无阻流量计算每口井的含气饱和度估计值；根据测井数据得到每口井的原始含气饱和度平均值；基于含气饱和度估计值和原始含气饱和度平均值，得到每口井的含气饱和度校正系数；采用含气饱和度校正系数对每口井的原始含气饱和度进行校正。2.根据权利要求1所述的方法，其特征在于，所述方法还包括：根据工区内各个井的无阻流量和原始含气饱和度平均值之间的关系，确定是否对原始含气饱和度进行校正。3.根据权利要求1所述的方法，其特征在于，所述基于地质信息和工程先验信息确定工区内无阻流量与含气饱和度之间的关系，包括：确定当无阻流量处于不同取值范围时对应的含气饱和度的取值范围。4.根据权利要求1所述的方法，其特征在于，所述根据测井数据得到每口井的原始含气饱和度平均值，包括：针对每口井，对测井数据在目标层的各个采样点测得的原始含气饱和度求算术平均，得到原始含气饱和度平均值。5.根据权利要求1所述的方法，其特征在于，所述确定无阻流量和含气饱和度估计值之间的关系曲线，包括：以含气饱和度估计值为无阻流量的幂函数为基础来拟合无阻流量与含气饱和度估计值之间的关系。6.一种基于无阻流量的含气饱和度校正装置，其特征在于，所述装置包括：先验信息获取单元，用于基于地质信息和工程先验信息确定工区内无阻流量与含气饱和度之间的关系；关系曲线确定单元，用于根据无阻流量与含气饱和度之间的关系，确定无阻流量和含气饱和度估计值之间的关系曲线；估计值计算单元，用于根据无阻流量计算每口井的含气饱和度估计值；原始平均值计算单元，用于根据测井数据得到每口井的原始含气饱和度平均值；校正系数确定单元，用于基于含气饱和度估计值和原始含气饱和度平均值，得到每口井的含气饱和度校正系数；校正单元，用于采用含气饱和度校正系数对每口井的原始含气饱和度进行校正。7.根据权利要求6所述的装置，其特征在于，所述装置还包括：数据评估单元，用于根据工区内各个井的无阻流量和原始含气饱和度平均值之间的关系，确定是否对原始含气饱和度进行校正。8.根据权利要求6所述的装置，其特征在于，所述先验信息获取单元具体用于：确定当无阻流量处于不同取值范围时对应的含气饱和度的取值范围。9.根据权利要求6所述的装置，其特征在于，所述原始平均值计算单元具体用于：2CN112630853A权　利　要　求　书2/2页针对每口井，对测井数据在目标层的各个采样点测得的原始含气饱和度求算术平均，得到原始含气饱和度平均值。10.根据权利要求6所述的装置，其特征在于，所述关系曲线确定单元具体用于：以含气饱和度估计值为无阻流量的幂函数为基础来拟合无阻流量与含气饱和度估计值之间的关系。3CN112630853A说　明　书1/5页基于无阻流量的含气饱和度校正方法和装置技术领域[0001]本发明属于石油勘探领域，更具体地，涉及一种基于无阻流量的饱和度校正方法和一种基于无阻流量的饱和度校正装置。背景技术[0002]含气饱和度，指在原始状态下，储层内天然气体积占连通孔隙体积的百分数。在石油勘探中，含气饱和度是一个非常重要的参数，其对于含气性预测等一系列工作具有重大指导意义。目前含气饱和度大多是根据阿尔奇公式进行计算，由于公式中的地层水矿化度等参数较难获取，测井解释人员只能根据经验进行计算，因此往往导致含气饱和度的计算不够准确。[0003]因此，需要一种能够获取更准确的含气饱和度的技术 方案 气瓶 现场处置方案 .pdf气瓶 现场处置方案 .doc见习基地管理方案.doc关于群访事件的化解方案建筑工地扬尘治理专项方案下载 。发明内容[0004]有鉴于此，本申请提出了一种能够对现有的含气饱和度进行校正的方法。本申请还提出了相应的装置。[0005]根据本申请的一方面，提供了一种基于无阻流量的含气饱和度校正方法，所述方法包括：基于地质信息和工程先验信息确定工区内无阻流量与含气饱和度之间的关系；根据无阻流量与含气饱和度之间的关系，确定无阻流量和含气饱和度估计值之间的关系曲线；根据无阻流量计算每口井的含气饱和度估计值；根据测井数据得到每口井的原始含气饱和度平均值；基于含气饱和度估计值和原始含气饱和度平均值，得到每口井的含气饱和度校正系数；采用含气饱和度校正系数对每口井的原始含气饱和度进行校正。[0006]在一种可能的实施方式中，所述方法还包括：根据工区内各个井的无阻流量和原始含气饱和度平均值之间的关系，确定是否对原始含气饱和度进行校正。[0007]在一种可能的实施方式中，所述基于地质信息和工程先验信息确定工区内无阻流量与含气饱和度之间的关系，包括：确定当无阻流量处于不同取值范围时对应的含气饱和度的取值范围。[0008]在一种可能的实施方式中，所述根据测井数据得到每口井的原始含气饱和度平均值，包括：针对每口井，对测井数据在目标层的各个采样点测得的原始含气饱和度求算术平均，得到原始含气饱和度平均值。[0009]在一种可能的实施方式中，所述确定无阻流量和含气饱和度估计值之间的关系曲线，包括：以含气饱和度估计值为无阻流量的幂函数为基础来拟合无阻流量与含气饱和度估计值之间的关系。[0010]根据本申请的另一方面，提供了一种基于无阻流量的含气饱和度校正装置，所述装置包括：先验信息获取单元，用于基于地质信息和工程先验信息确定工区内无阻流量与含气饱和度之间的关系；关系曲线确定单元，用于根据无阻流量与含气饱和度之间的关系，确定无阻流量和含气饱和度估计值之间的关系曲线；估计值计算单元，用于根据无阻流量4CN112630853A说　明　书2/5页计算每口井的含气饱和度估计值；原始平均值计算单元，用于根据测井数据得到每口井的原始含气饱和度平均值；校正系数确定单元，用于基于含气饱和度估计值和原始含气饱和度平均值，得到每口井的含气饱和度校正系数；校正单元，用于采用含气饱和度校正系数对每口井的原始含气饱和度进行校正。[0011]在一种可能的实施方式中，所述装置还包括：数据评估单元，用于根据工区内各个井的无阻流量和原始含气饱和度平均值之间的关系，确定是否对原始含气饱和度进行校正。[0012]在一种可能的实施方式中，所述先验信息获取单元具体用于：确定当无阻流量处于不同取值范围时对应的含气饱和度的取值范围。[0013]在一种可能的实施方式中，所述原始平均值计算单元具体用于：针对每口井，对测井数据在目标层的各个采样点测得的原始含气饱和度求算术平均，得到原始含气饱和度平均值。[0014]在一种可能的实施方式中，所述关系曲线确定单元具体用于：以含气饱和度估计值为无阻流量的幂函数为基础来拟合无阻流量与含气饱和度估计值之间的关系。[0015]根据本申请，针对常规的含气饱和度曲线计算，以无阻流量与含气饱和度的正相关关系为基础，在保证含气饱和度曲线方差不变的情况下进行均值校正，获得与开发生产资料及地质认识相符的含气饱和度曲线，能有效指导含气性预测等后续一系列工作，在气藏发育地区具有极大的推广应用价值。附图说明[0016]通过结合附图对本申请示例性实施方式进行更详细的描述，本申请的上述以及其它目的、特征和优势将变得更加明显，其中，在本申请示例性实施方式中，相同的参考标号通常代表相同部件。[0017]图1示出根据本申请的一个实施例的基于无阻流量的含气饱和度校正方法的流程图。[0018]图2示出根据本申请的一个示例性实施例的校正路径示意图。[0019]图3示出根据本申请的一个实施例的基于无阻流量的含气饱和度校正装置的结构框图。[0020]图4示出某工区在校正前不同产气井的无阻流量和原始含气饱和度的直方图。[0021]图5示出在某工区根据地质信息和工程先验信息得到的无阻流量和含气饱和度估计值之间的关系曲线。[0022]图6示出根据本申请进行校正后某工区的不同产气井的无阻流量和含气饱和度的直方图。具体实施方式[0023]下面将参照附图更详细地描述本申请的优选实施方式。虽然附图中显示了本申请的优选实施方式，然而应该理解，可以以各种形式实现本申请而不应被这里阐述的实施方式所限制。相反，提供这些实施方式是为了使本申请更加透彻和完整，并且能够将本申请的范围完整地传达给本领域的技术人员。5CN112630853A说　明　书3/5页[0024]请参见图1。图1示出根据本申请的一个实施例的基于无阻流量的含气饱和度校正方法的流程图。如图所示，所述方法包括下列步骤。[0025]步骤102，基于地质信息和工程先验信息确定工区内无阻流量A与含气饱和度sg之间的关系。[0026]此处无阻流量A和含气饱和度sg之间的关系是该工区的工作人员根据之前的经验总结得到的。由于采纳工作人员经对该工区长年累月的勘探和认识(即对地质信息和工程先验信息的了解)而总结出的经验来确定无阻流量A和含气饱和度sg之间的关系，使得得到的无阻流量A和含气饱和度sg间的关系能准确表达该工区的实际情况，为后续进行校正打下良好基础。[0027]在一种可能的实施方式中，步骤102可以包括确定当无阻流量A处于不同取值范围时对应的含气饱和度sg的取值范围。换言之，根据经验得到的无阻流量A和含气饱和度sg之间的关系，可以是一种大致的范围对应关系。[0028]步骤104，根据无阻流量A与含气饱和度sg之间的关系，确定无阻流量A和含气饱和度估计值sg_ave之间的关系曲线。[0029]在一种可能的实现方式中，可以以含气饱和度估计值sg_ave为无阻流量A的幂函数为基础来拟合无阻流量A与含气饱和度估计值sg_ave之间的关系。[0030]由于无阻流量A是敞开井口的放喷方法得到的产量，可用来衡量气井的生产能力及进行产气井间生产能力的比较，因此无阻流量A与含气饱和度sg之间应该是正相关的。而在此基础上，发明人进行了进一步研究，并对从多个工区采集的数据进行深入分析，发现将含气饱和度作为无阻流量A的幂函数，可以较准确地表达二者之间的关系。因此，可以以此为基础来拟合二者的关系曲线。[0031]步骤106，根据无阻流量A计算每口井的含气饱和度估计值sg_ave。[0032]针对每口井，可以把其无阻流量A带入步骤104中得到的无阻流量A和含气饱和度估计值sg_ave之间的关系，以得到该井的含气饱和度估计值sg_ave。[0033]步骤108，根据测井数据得到每口井的原始含气饱和度平均值sg_orign_ave。[0034]在一种可能的实施方式中，针对每口井，对测井数据在目标层的各个采样点测得的原始含气饱和度sg_orign求算术平均，得到原始含气饱和度平均值sg_orign_ave。[0035]步骤110，基于含气饱和度估计值sg_ave和原始含气饱和度平均值sg_orign_ave，得到每口井的含气饱和度校正系数sg_check。[0036]可以令sg_check＝sg_ave-sg_orign_ave，得到sg_check。[0037]步骤112，采用含气饱和度校正系数对每口井的原始含气饱和度进行校正。[0038]在对原始含气饱和度进行校正后，可制作不同产气井的无阻流量和校正后含气饱和度的直方图，以确认无阻流量与含气饱和度对应关系符合正相关。[0039]在上述实施例中，针对常规的含气饱和度曲线计算，以无阻流量与含气饱和度的正相关关系为基础，在保证含气饱和度曲线方差不变的情况下进行均值校正，获得与开发生产资料及地质认识相符的含气饱和度曲线，能有效指导含气性预测等后续一系列工作，在气藏发育地区具有极大的推广应用价值[0040]在一种可能的实现方式中，在上述步骤102前，可先根据工区内各个井的无阻流量A和含气饱和度sg_orign之间的关系，确定是否对含气饱和度sg_orign进行校正。例如，可6CN112630853A说　明　书4/5页以制作不同产气井的无阻流量A和原始含气饱和度平均值sg_orign_ave的直方图，如果无阻流量A与原始含气饱和度平均值sg_orign_ave符合正相关，则可认为原始含气饱和度sg_orign是基本准确的，可不做校正；如果二者不符合正相关，则可认为原始含气饱和度sg_orign有误，可对其进行校正。[0041]图2示出根据本申请的一个示例性实施例的校正路径示意图。如图所示，可根据工区的地质信息、工程先验信息(202)得到无阻流量和含气饱和度的实际关系(204)。将无阻流量(206)带入无阻流量和含气饱和度的实际关系(204)，得到每口井的含气饱和度估计值sg_ave(208)。根据原始含气饱和度(210)得到原始含气饱和度平均值sg_orign_ave(212)。sg_ave(208)和sg_orign_ave(212)二者相减，得到含气饱和度校正系数sg_check(214)。采用sg_check(214)对原始含气饱和度(210)进行校正，得到校正后的含气饱和度(216)。[0042]图3示出根据本申请的一个实施例的基于无阻流量的含气饱和度校正装置的结构框图。如图所示，该装置包括先验信息获取单元302、关系曲线确定单元304、估计值计算单元306、原始平均值计算单元308、校正系数确定单元310和校正单元312。[0043]先验信息获取单元302用于基于地质信息和工程先验信息确定工区内无阻流量与含气饱和度之间的关系。[0044]关系曲线确定单元304用于根据无阻流量与含气饱和度之间的关系，确定无阻流量和含气饱和度估计值之间的关系曲线。[0045]估计值计算单元306用于根据无阻流量计算每口井的含气饱和度估计值。[0046]原始平均值计算单元308用于根据测井数据得到每口井的原始含气饱和度平均值。[0047]校正系数确定单元310用于基于含气饱和度估计值和原始含气饱和度平均值，得到每口井的含气饱和度校正系数。[0048]校正单元312用于采用含气饱和度校正系数对每口井的原始含气饱和度进行校正。[0049]在一种可能的实施方式中，所述装置还包括：数据评估单元，用于根据工区内各个井的无阻流量和原始含气饱和度平均值之间的关系，确定是否对原始含气饱和度进行校正。[0050]在一种可能的实施方式中，所述先验信息获取单元具体用于：确定当无阻流量处于不同取值范围时对应的含气饱和度的取值范围。[0051]在一种可能的实施方式中，所述原始平均值计算单元具体用于：针对每口井，对测井数据在目标层的各个采样点测得的原始含气饱和度求算术平均，得到原始含气饱和度平均值。[0052]在一种可能的实施方式中，所述关系曲线确定单元具体用于：以含气饱和度估计值为无阻流量的幂函数为基础来拟合无阻流量与含气饱和度估计值之间的关系[0053]应用示例[0054]在本应用示例中，先基于工区测井解释成果，得到工区内不同产气井的无阻流量A和原始含气饱和度sg_orign的直方图，如图4所示。从图4中可以看出，得到的无阻流量A和含气饱和度sg_orign之间不是正相关关系，因此，可认为需要对图2中的含气饱和度sg_orign进行校正。7CN112630853A说　明　书5/5页[0055]在本示例中，本工区的作业人员经过对该工区长年累月的勘探和认识，认为本工区无阻流量A和含气饱和度sg之间具有如下关系：[0056]当无阻流量A>10，含气饱和度sg>70％；[0057]当3<无阻流量A<10，55％<含气饱和度sg<70％；[0058]当无阻流量A<3，40％<含气饱和度sg<55％。[0059]基于上述关系，并以含气饱和度估计值sg_ave为无阻流量A的幂函数为基础，得到如图5所示的无阻流量A和含气饱和度估计值sg_ave(即图5中纵坐标所表示的含气饱和度)之间的关系曲线：sg_ave＝32*A1/6+20。然后将各口井的无阻流量A带入sg_ave＝32*A1/6+20，得到每口井的含气饱和度估计值sg_ave。[0060]此外，针对每口井，对测井数据在目标层的各个采样点测得的原始含气饱和度sg_orign求算术平均，得到原始含气饱和度平均值sg_orign_ave，并令sg_check＝sg_ave-sg_orign_ave，获得每口井的含气饱和度校正系数sg_check。[0061]采用含气饱和度校正系数sg_check对每口井的原始含气饱和度sg_orign进行校正，从而可以在保证原始含气饱和度曲线方差不变的情况下进行均值校正。可根据校正结果制作不同产气井的无阻流量A和校正后的含气饱和度sg_checked的直方图，如图6所示。从图6中可以看出，根据本申请校正后的含气饱和度sg_checked与无阻流量A具有良好的正相关对应关系较好，符合地质认识。[0062]本申请可以是系统、方法和/或计算机程序产品。计算机程序产品可以包括计算机可读存储介质，其上载有用于使处理器实现本申请的各个方面的计算机可读程序指令。[0063]计算机可读存储介质可以是可以保持和存储由指令执行设备使用的指令的有形设备。计算机可读存储介质例如可以是――但不限于――电存储设备、磁存储设备、光存储设备、电磁存储设备、半导体存储设备或者上述的任意合适的组合。计算机可读存储介质的更具体的例子(非穷举的列表)包括：便携式计算机盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、静态随机存取存储器(SRAM)、便携式压缩盘只读存储器(CD-ROM)、数字多功能盘(DVD)、记忆棒、软盘、机械编码设备、例如其上存储有指令的打孔卡或凹槽内凸起结构、以及上述的任意合适的组合。这里所使用的计算机可读存储介质不被解释为瞬时信号本身，诸如无线电波或者其他自由传播的电磁波、通过波导或其他传输媒介传播的电磁波(例如，通过光纤电缆的光脉冲)、或者通过电线传输的电信号。[0064]这里参照根据本申请实施例的方法、装置(系统)和计算机程序产品的流程图和/或框图描述了本申请的各个方面。应当理解，流程图和/或框图的每个方框以及流程图和/或框图中各方框的组合，都可以由计算机可读程序指令实现。[0065]以上已经描述了本申请的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择，旨在最好地解释各实施例的原理、实际应用或对市场中的技术的技术改进，或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。8CN112630853A说　明　书　附　图1/3页图1图29CN112630853A说　明　书　附　图2/3页图3图410CN112630853A说　明　书　附　图3/3页图5图611