CN111005715A - 确定气井产量的方法、装置及存储介质 - Google Patents
确定气井产量的方法、装置及存储介质 Download PDFInfo
- Publication number
- CN111005715A CN111005715A CN201811168646.5A CN201811168646A CN111005715A CN 111005715 A CN111005715 A CN 111005715A CN 201811168646 A CN201811168646 A CN 201811168646A CN 111005715 A CN111005715 A CN 111005715A
- Authority
- CN
- China
- Prior art keywords
- coefficient
- gas well
- determining
- equation
- gas
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 43
- 238000003860 storage Methods 0.000 title claims abstract description 10
- 238000004519 manufacturing process Methods 0.000 claims abstract description 151
- 230000003068 static effect Effects 0.000 claims abstract description 125
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 123
- 238000012360 testing method Methods 0.000 claims description 96
- 238000004364 calculation method Methods 0.000 claims description 13
- 239000000126 substance Substances 0.000 claims description 5
- 239000007789 gas Substances 0.000 description 401
- 230000002093 peripheral effect Effects 0.000 description 10
- 230000001133 acceleration Effects 0.000 description 9
- 230000008569 process Effects 0.000 description 9
- 238000012545 processing Methods 0.000 description 9
- 238000004891 communication Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 4
- 238000013473 artificial intelligence Methods 0.000 description 2
- 239000000919 ceramic Substances 0.000 description 2
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000003345 natural gas Substances 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mining & Mineral Resources (AREA)
- Geology (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Databases & Information Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- Computational Mathematics (AREA)
- Algebra (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Measuring Volume Flow (AREA)
Abstract
本发明公开了一种确定气井产量的方法、装置及存储介质,属于气井开采领域。该方法包括:根据多个气井中每个气井对应的第一无阻流量,确定目标第三系数;获取目标气井的地层静压、井底流压和产气量;根据目标第三系数、目标气井的地层静压、井底流压和产气量,通过第三产能方程得到无阻流量,将所得到的无阻流量确定为第二无阻流量;根据第一产能方程、第二产能方程、第二无阻流量、目标气井的地层静压、井底流压和产气量,确定目标气井对应的第一系数和第二系数;根据第一产能方程、目标气井对应的第一系数和第二系数,确定目标气井的合理产量。本发明可以提高确定目标气井的无阻流量的效率和确定出的目标气井的气井产量的准确度。
Description
技术领域
本发明涉及气井开采领域,特别涉及一种确定气井产量的方法、装置及存储介质。
背景技术
为了了解气井的产气能力,常通过产能试井求得气井无阻流量,再根据气井无阻流量确定气井产量。
目前,在通过产能试井求取气井无阻流量,并据此确定气井产量时,先获取气井的多组气井产能试井数据,再根据该气井的多组气井产能试井数据确定产能方程中的系数。然后,令该产能方程中的井底流压为0.101MPa(兆帕),并将确定出的系数和该气井的地层静压代入该产能方程中,得到该气井的无阻流量。之后,将该气井的无阻流量乘以预设的比例系数,得到该气井的合理产量。
然而,上述方式中预设的比例系数通常是根据试井经验所得,因而会导致据此确定出的气井产量的准确度较低。
发明内容
本发明实施例提供了一种确定气井产量的方法、装置及存储介质,可以解决相关技术中确定出的气井产量的准确度较低的问题。所述技术方案如下:
根据本发明实施例的第一方面,提供一种确定气井产量的方法,所述方法包括:
获取多个气井中每个气井的地层静压和多组试井数据,所述多组试井数据中每组试井数据包括井底流压和产气量;
对于所述多个气井中的每个气井,根据第一产能方程、所述气井的地层静压和多组试井数据,确定所述气井对应的第一系数和第二系数;
根据所述气井的地层静压、所述气井对应的第一系数和第二系数,通过第二产能方程得到无阻流量,将所得到的无阻流量确定为所述气井对应的第一无阻流量;
根据所述多个气井中每个气井对应的第一无阻流量,确定目标第三系数;
获取目标气井的地层静压、井底流压和产气量;
根据所述目标第三系数、所述目标气井的地层静压、井底流压和产气量,通过第三产能方程得到无阻流量,将所得到的无阻流量确定为第二无阻流量;
根据所述第一产能方程、所述第二产能方程、所述第二无阻流量、所述目标气井的地层静压、井底流压和产气量,确定所述目标气井对应的第一系数和第二系数;
根据所述第一产能方程、所述目标气井对应的第一系数和第二系数,确定所述目标气井的合理产量;
其中,所述第一产能方程为pR 2-pwf 2=Aqg+Bqg 2,所述第二产能方程为pR 2-0.1012=AqAOF+BqAOF 2,所述第三产能方程为pD=αqD+(1-α)qD 2,所述A为第一系数,所述B为第二系数,所述α为第三系数;或者,所述第一产能方程为qg=C(pR 2-pwf 2)n,所述第二产能方程为qAOF=C(pR 2-0.1012)n,所述第三产能方程为lgqD=klgpD+b,所述C为第一系数,所述n为第二系数,所述K为第三系数中的第一子系数,所述b为第三系数中的第二子系数;
可选地,所述根据第一产能方程、所述气井的地层静压和多组试井数据,确定所述气井对应的第一系数和第二系数,包括:
根据所述气井的地层静压和多组试井数据,确定所述第一公式对应的第一曲线;
将所述第一曲线的斜率确定为所述气井对应的第二系数,将所述第一曲线的截距确定为所述气井对应的第一系数。
可选地,所述根据第一产能方程、所述气井的地层静压和多组试井数据,确定所述气井对应的第一系数和第二系数,包括:
根据所述气井的地层静压和多组试井数据,确定所述第二公式对应的第二曲线;
将所述第二曲线的斜率的倒数确定为所述气井对应的第二系数;
获取所述第二曲线上的任一坐标点的坐标值中的第一数值和第二数值,所述第一数值为pR 2-pwf 2,所述第二数值为qg;
根据所述第一数值、所述第二数值和所述气井对应的第二系数,通过所述第一产能方程得到第一系数,将所得到的第一系数确定为所述气井对应的第一系数。
可选地,所述根据所述多个气井中每个气井对应的第一无阻流量,确定目标第三系数,包括:
当所述第三产能方程为PD=αqD+(1-α)qD 2时,对于所述多个气井中的每个气井,根据所述气井对应的第一无阻流量、第一系数和第二系数,通过第三公式得到第三系数,将所得到的第三系数确定为所述气井对应的第三系数,所述第三公式为
将所述多个气井对应的第三系数的平均值确定为目标第三系数。
可选地,所述根据所述多个气井中每个气井对应的第一无阻流量,确定目标第三系数,包括:
当所述第三产能方程为lgqD=klgpD+b时,对于所述多个气井中的每个气井,根据所述气井的地层静压、多组试井数据和所述气井对应的第一无阻流量,确定所述气井对应的多组参数,所述多组参数中的每组参数包括pD和qD;
根据所述多个气井中每个气井对应的多组参数,确定多个坐标点,所述多个坐标点中每个坐标点的坐标值为(lgpD,lgqD);
对所述多个坐标点进行线性拟合,得到第三曲线;
将所述第三曲线的斜率确定为目标第三系数中的第一子系数,将所述第三曲线的截距确定为目标第三系数中的第二子系数。
可选地,所述根据所述第一产能方程、所述第二产能方程、所述第二无阻流量、所述目标气井的地层静压、井底流压和产气量,确定所述目标气井对应的第一系数和第二系数,包括:
当所述第一产能方程为pR 2-pwf 2=Aqg+Bqg 2,所述第二产能方程为pR 2-0.1012=AqAOF+BqAOF 2时,对所述第一产能方程和所述第二产能方程的联立方程式求解,得到第四公式和第五公式,所述第四公式为所述第五公式为
根据所述第二无阻流量、所述目标气井的地层静压、井底流压和产气量,通过所述第四公式得到第一系数,并通过所述第五公式得到第二系数,将所得到的第一系数确定为所述目标气井对应的第一系数,将所得到的第二系数确定为所述目标气井对应的第二系数。
可选地,所述根据所述第一产能方程、所述第二产能方程、所述第二无阻流量、所述目标气井的地层静压、井底流压和产气量,确定所述目标气井对应的第一系数和第二系数,包括:
当所述第一产能方程为qg=C(pR 2-pwf 2)n,所述第二产能方程为qAOF=C(pR 2-0.1012)n时,对所述第一产能方程和所述第二产能方程的联立方程式求解,得到第六公式和第七公式,所述第六公式为所述第七公式为
根据所述第二无阻流量、所述目标气井的地层静压、井底流压和产气量,通过所述第六公式得到第一系数,并通过所述第七公式得到第二系数,将所得到的第一系数确定为所述目标气井对应的第一系数,将所得到的第二系数确定为所述目标气井对应的第二系数。
根据本发明实施例的第二方面,提供一种确定气井产量的装置,所述装置包括:
第一获取模块,用于获取多个气井中每个气井的地层静压和多组试井数据,所述多组试井数据中每组试井数据包括井底流压和产气量;
第一确定模块,用于对于所述多个气井中的每个气井,根据第一产能方程、所述气井的地层静压和多组试井数据,确定所述气井对应的第一系数和第二系数;
第二确定模块,用于根据所述气井的地层静压、所述气井对应的第一系数和第二系数,通过第二产能方程得到无阻流量,将所得到的无阻流量确定为所述气井对应的第一无阻流量;
第三确定模块,用于根据所述多个气井中每个气井对应的第一无阻流量,确定目标第三系数;
第二获取模块,用于获取目标气井的地层静压、井底流压和产气量;
第四确定模块,用于根据所述目标第三系数、所述目标气井的地层静压、井底流压和产气量,通过第三产能方程得到无阻流量,将所得到的无阻流量确定为第二无阻流量;
第五确定模块,用于根据所述第一产能方程、所述第二产能方程、所述第二无阻流量、所述目标气井的地层静压、井底流压和产气量,确定所述目标气井对应的第一系数和第二系数;
第六确定模块,用于根据所述第一产能方程、所述目标气井对应的第一系数和第二系数,确定所述目标气井的合理产量;
其中,所述第一产能方程为pR 2-pwf 2=Aqg+Bqg 2,所述第二产能方程为pR 2-0.1012=AqAOF+BqAOF 2,所述第三产能方程为pD=αqD+(1-α)qD 2,所述A为第一系数,所述B为第二系数,所述α为第三系数;或者,所述第一产能方程为qg=C(pR 2-pwf 2)n,所述第二产能方程为qAOF=C(pR 2-0.1012)n,所述第三产能方程为lgqD=klgpD+b,所述C为第一系数,所述n为第二系数,所述K为第三系数中的第一子系数,所述b为第三系数中的第二子系数;
可选地,所述第一确定模块,包括:
第一确定子模块,用于根据所述气井的地层静压和多组试井数据,确定所述第一公式对应的第一曲线;
第二确定子模块,用于将所述第一曲线的斜率确定为所述气井对应的第二系数,将所述第一曲线的截距确定为所述气井对应的第一系数。
可选地,所述第一确定模块,包括:
第二计算子模块,用于当所述第一产能方程为qg=C(pR 2-pwf 2)n时,将所述第一产能方程的两端同取对数,得到第二公式,所述第二公式为所述y2=lg(pR 2-pwf 2),所述x2=lgqg;
第三确定子模块,用于根据所述气井的地层静压和多组试井数据,确定所述第二公式对应的第二曲线;
第四确定子模块,用于将所述第二曲线的斜率的倒数确定为所述气井对应的第二系数;
第一获取子模块,用于获取所述第二曲线上的任一坐标点的坐标值中的第一数值和第二数值,所述第一数值为pR 2-pwf 2,所述第二数值为qg;
第五确定子模块,用于根据所述第一数值、所述第二数值和所述气井对应的第二系数,通过所述第一产能方程得到第一系数,将所得到的第一系数确定为所述气井对应的第一系数。
可选地,所述第三确定模块,包括:
第六确定子模块,用于当所述第三产能方程为PD=αqD+(1-α)qD 2时,对于所述多个气井中的每个气井,根据所述气井对应的第一无阻流量、第一系数和第二系数,通过第三公式得到第三系数,将所得到的第三系数确定为所述气井对应的第三系数,所述第三公式为
第七确定子模块,用于将所述多个气井对应的第三系数的平均值确定为目标第三系数。
可选地,所述第三确定模块,包括:
第八确定子模块,用于当所述第三产能方程为lgqD=klgpD+b时,对于所述多个气井中的每个气井,根据所述气井的地层静压、多组试井数据和所述气井对应的第一无阻流量,确定所述气井对应的多组参数,所述多组参数中的每组参数包括pD和qD;
第九确定子模块,用于根据所述多个气井中每个气井对应的多组参数,确定多个坐标点,所述多个坐标点中每个坐标点的坐标值为(lgpD,lgqD);
拟合子模块,用于对所述多个坐标点进行线性拟合,得到第三曲线;
第十确定子模块,用于将所述第三曲线的斜率确定为目标第三系数中的第一子系数,将所述第三曲线的截距确定为目标第三系数中的第二子系数。
可选地,所述第五确定模块,包括:
第三计算子模块,用于当所述第一产能方程为pR 2-pwf 2=Aqg+Bqg 2,所述第二产能方程为pR 2-0.1012=AqAOF+BqAOF 2时,对所述第一产能方程和所述第二产能方程的联立方程式求解,得到第四公式和第五公式,所述第四公式为所述第五公式为
第十一确定子模块,用于根据所述第二无阻流量、所述目标气井的地层静压、井底流压和产气量,通过所述第四公式得到第一系数,并通过所述第五公式得到第二系数,将所得到的第一系数确定为所述目标气井对应的第一系数,将所得到的第二系数确定为所述目标气井对应的第二系数。
可选地,所述第五确定模块,包括:
第四计算子模块,用于当所述第一产能方程为qg=C(pR 2-pwf 2)n,所述第二产能方程为qAOF=C(pR 2-0.1012)n时,对所述第一产能方程和所述第二产能方程的联立方程式求解,得到第六公式和第七公式,所述第六公式为所述第七公式为
第十二确定子模块,用于根据所述第二无阻流量、所述目标气井的地层静压、井底流压和产气量,通过所述第六公式得到第一系数,并通过所述第七公式得到第二系数,将所得到的第一系数确定为所述目标气井对应的第一系数,将所得到的第二系数确定为所述目标气井对应的第二系数。
根据本发明实施例的第三方面,提供一种确定气井产量的装置,所述装置包括处理器、存储器及存储在所述存储器上并可在所述处理器上运行的程序代码,所述处理器执行所述程序代码时实现上述第一方面所述的确定气井产量的方法。
根据本发明实施例的第四方面,提供了一种计算机可读存储介质,所述存储介质上存储有指令,所述指令被处理器执行时实现上述第一方面所述的确定气井产量的方法的步骤。
本发明实施例提供的技术方案至少可以带来以下有益效果:本发明实施例中,获取多个气井中每个气井的地层静压和多组试井数据,该多组试井数据中每组试井数据包括井底流压和产气量。然后,对于该多个气井中的每个气井,根据第一产能方程、这个气井的地层静压和多组试井数据,确定这个气井对应的第一系数和第二系数,再根据这个气井的地层静压、这个气井对应的第一系数和第二系数,通过第二产能方程得到无阻流量,将所得到的无阻流量确定为这个气井对应的第一无阻流量。然后,根据该多个气井中每个气井对应的第一无阻流量,确定目标第三系数。之后,先获取目标气井的地层静压、井底流压和产气量,再根据目标第三系数、目标气井的地层静压、井底流压和产气量,通过第三产能方程得到无阻流量,将所得到的无阻流量确定为第二无阻流量。由于在确定目标第三系数后,仅需要获取目标气井的地层静压以及一组井底流压和产气量就可以确定第二无阻流量,即目标气井的无阻流量,所以提高了确定目标气井的无阻流量的效率,降低了确定目标气井的无阻流量的成本。最后,根据第一产能方程、第二产能方程、第二无阻流量、目标气井的地层静压、井底流压和产气量,确定目标气井对应的第一系数和第二系数,再根据第一产能方程、目标气井对应的第一系数和第二系数,确定目标气井的合理产量,从而提高了确定出的气井产量的准确度。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的一种确定气井产量的方法的流程图;
图2是本发明实施例提供的一种确定气井产量的装置的结构示意图;
图3是本发明实施例提供的另一种确定气井产量的装置的结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明实施方式作进一步地详细描述。
图1是本发明实施例提供的一种确定气井产量的方法的流程图。参见图1,该方法包括以下步骤:
步骤101:获取多个气井中每个气井的地层静压和多组试井数据,该多组试井数据中每组试井数据包括井底流压和产气量。
需要说明的是,地层静压也称油层压力,是指在关井后,待压力恢复到稳定状态时所测得的油层中部的压力。井底流压指的是油、气从地层流到井底后剩余的压力。试井指的是气井产能试井,是在地层静压一定的条件下,测试气井在不同的井底流压时的产气量。
另外,该多个气井可以为位于同一区域的多个气井,如可以为位于同一省、市等的多个气井。
步骤102:对于该多个气井中的每个气井,根据第一产能方程、这个气井的地层静压和多组试井数据,确定这个气井对应的第一系数和第二系数。
需要说明的是,第一产能方程可以为二项式产能方程,此时,第一产能方程为pR 2-pwf 2=Aqg+Bqg 2,A为第一系数,B为第二系数;或者,第一产能方程也可以为指数式产能方程,此时,第一产能方程为qg=C(pR 2-pwf 2)n,C为第一系数,n为第二系数。其中,pR为地层静压,pwf为井底流压,qg为产气量。
其中,根据第一产能方程、这个气井的地层静压和多组试井数据,确定这个气井对应的第一系数和第二系数的操作可以通过以下两种方式实现:
第一种方式:当第一产能方程为pR 2-pwf 2=Aqg+Bqg 2时,将第一产能方程的两端同除以qg,得到第一公式,第一公式为y1=A+Bx1,其中,x1=qg;根据这个气井的地层静压和多组试井数据,确定第一公式对应的第一曲线;将第一曲线的斜率确定为这个气井对应的第二系数,将第一曲线的截距确定为这个气井对应的第一系数。
其中,根据这个气井的地层静压和多组试井数据,确定第一公式对应的第一曲线的实现过程可以是:将这个气井的地层静压和多组试井数据中的每组试井数据代入第一公式中,得到多个坐标值;对该多个坐标值指示的多个坐标点进行线性拟合,得到第一公式对应的第一曲线。
需要说明的是,该多个坐标值与这个气井的多组试井数据一一对应,该多个坐标值中的某个坐标值是根据这个气井的地层静压和这个坐标值对应的一组试井数据得到。
第二种方式:当第一产能方程为qg=C(pR 2-pwf 2)n时,将第一产能方程的两端同取对数,得到第二公式,第二公式为y2=lg(pR 2-pwf 2),x2=lgqg;根据这个气井的地层静压和多组试井数据,确定第二公式对应的第二曲线;将第二曲线的斜率的倒数确定为这个气井对应的第二系数;获取第二曲线上的任一坐标点的坐标值中的第一数值和第二数值,第一数值为pR 2-pwf 2,第二数值为qg;根据第一数值、第二数值和这个气井对应的第二系数,通过第一产能方程得到第一系数,将所得到的第一系数确定为这个气井对应的第一系数。
其中,根据这个气井的地层静压和多组试井数据,确定第二公式对应的第二曲线的实现过程可以是:将这个气井的地层静压和多组试井数据中的每组试井数据代入第二公式中,得到多个坐标值;对该多个坐标值指示的多个坐标点进行线性拟合,得到第二公式对应的第二曲线。
需要说明的是,该多个坐标值与这个气井的多组试井数据一一对应,该多个坐标值中的某个坐标值是根据这个气井的地层静压和这个坐标值对应的一组试井数据得到。
步骤103:根据这个气井的地层静压、这个气井对应的第一系数和第二系数,通过第二产能方程得到无阻流量,将所得到的无阻流量确定为这个气井对应的第一无阻流量。
需要说明的是,由于第一产能方程可以为二项式产能方程,也可以为指数式产能方程,而第二产能方程是令第一产能方程中的井底流压为0.101MPa时得到的,所以,当第一产能方程为二项式产能方程时,第二产能方程可以为pR 2-0.1012=AqAOF+BqAOF 2;当第一产能方程为指数式产能方程时,第二产能方程可以为qAOF=C(pR 2-0.1012)n。其中,qAOF为无阻流量。
另外,当第二产能方程为pR 2-0.1012=AqAOF+BqAOF 2时,可以将这个气井的地层静压、这个气井对应的第一系数(A)和第二系数(B)代入第二产能方程中,得到关于qAOF的一元二次方程,求解该一元二次方程即可得到无阻流量。
另外,当第二产能方程为qAOF=C(pR 2-0.1012)n时,可以将这个气井的地层静压、这个气井对应的第一系数(C)和第二系数(n)代入第二产能方程中,求得无阻流量。
通过上述步骤102和步骤103可以得到该多个气井中每个气井对应的第一无阻流量,在得到该多个气井中每个气井对应的第一无阻流量后,可以继续执行如下步骤104。
步骤104:根据该多个气井中每个气井对应的第一无阻流量,确定目标第三系数。
需要说明的是,目标第三系数为第三产能方程中的系数。当第一产能方程为二项式产能方程时,第三产能方程可以为pD=αqD+(1-α)qD 2,α为第三系数;当第一产能方程为指数式产能方程时,第三产能方程为lgqD=klgpD+b,K为第三系数中的第一子系数,b为第三系数中的第二子系数。其中,pD为无因次压力,qD为无因次产量,
另外,当第一产能方程为二项式产能方程时,第三产能方程可以是令第二产能方程(pR 2-0.1012=AqAOF+BqAOF 2)中的pR 2-0.1012≈pR 2,并用第一产能方程(pR 2-pwf 2=Aqg+Bqg 2)除以第二产能方程(pR 2-0.1012=AqAOF+BqAOF 2)得到的。
具体地,步骤104的实现过程可以通过以下两种方式实现:
第一种方式:当第三产能方程为PD=αqD+(1-α)qD 2时,对于该多个气井中的每个气井,根据这个气井对应的第一无阻流量、第一系数和第二系数,通过第三公式得到第三系数,将所得到的第三系数确定为这个气井对应的第三系数,第三公式为将该多个气井对应的第三系数的平均值确定为目标第三系数。
第二种方式:当第三产能方程为lgqD=klgpD+b时,对于该多个气井中的每个气井,根据这个气井的地层静压、多组试井数据和这个气井对应的第一无阻流量,确定这个气井对应的多组参数,该多组参数中的每组参数包括pD和qD;根据该多个气井中每个气井对应的多组参数,确定多个坐标点,该多个坐标点中每个坐标点的坐标值为(lgpD,lgqD);对该多个坐标点进行线性拟合,得到第三曲线;将第三曲线的斜率确定为目标第三系数中的第一子系数,将第三曲线的截距确定为目标第三系数中的第二子系数。
需要说明的是,这个气井对应的多组参数与这个气井对应的多组试井数据一一对应,这个气井对应的多组参数中的某组参数是根据这个气井的地层静压、这个气井对应的第一无阻流量和这组参数对应的一组试井数据得到。
另外,该多个坐标点与这个气井对应的多组参数一一对应,该多个坐标点中某个坐标点的坐标值是根据这个坐标点对应的一组参数得到。
再者,实际应用中,第三曲线的曲线方程即为第三产能方程(lgqD=klgpD+b)。
步骤105:获取目标气井的地层静压、井底流压和产气量。
需要说明的是,目标气井是与该多个气井处于同一区域的气井,目标气井可以是该多个气井中的任一气井,也可以是除该多个气井之外的任一气井。目标气井的地层静压、井底流压和产气量可以通过对目标气井进行试井得到,且目标气井的井底流压和产气量可以为在地层静压一定的条件下得到的一组试井数据。
值得说明的是,本发明实施例中通过步骤104根据该多个气井中每个气井对应的第一无阻流量,确定目标第三系数后,通过步骤105仅需要获取目标气井的地层静压以及一组井底流压和产气量,后续即可通过步骤106据此求得目标气井的无阻流量(第二无阻流量)。这样,不仅可以提高确定目标气井的无阻流量的效率,而且可以降低确定目标气井的无阻流量的成本。
步骤106:根据目标第三系数、目标气井的地层静压、井底流压和产气量,通过第三产能方程得到无阻流量,将所得到的无阻流量确定为第二无阻流量。
需要说明的是,当第三产能方程为PD=αqD+(1-α)qD 2时,可以将第三产能方程看为关于qD的一元二次方程,继而可以得到该一元二次方程的根的表达式然后将目标气井的地层静压和井底流压代入中,得到pD,将pD和目标第三系数(α)代入该一元二次方程的根的表达式中,得到qD。之后,根据得到qAOF的表达式将qD和目标气井的产气量代入qAOF的表达式中,得到无阻流量作为第二无阻流量。
另外,当第三产能方程为lgqD=klgpD+b时,可以将代入第三产能方程中,得到qAOF的对数表达式然后对该qAOF的对数表达式进行变形,得到qAOF的指数表达式之后,将目标气井的地层静压和井底流压代入中,得到pD,将pD、目标第三系数中的第一子系数(K)、目标第三系数中的第二子系数(b)和目标气井的产气量代入该qAOF的指数表达式中,得到无阻流量作为第二无阻流量。
步骤107:根据第一产能方程、第二产能方程、第二无阻流量、目标气井的地层静压、井底流压和产气量,确定目标气井对应的第一系数和第二系数。
具体地,步骤107的实现过程可以通过以下两种方式实现:
第一种方式:当第一产能方程为pR 2-pwf 2=Aqg+Bqg 2,第二产能方程为pR 2-0.1012=AqAOF+BqAOF 2时,对第一产能方程和第二产能方程的联立方程式求解,得到第四公式和第五公式,第四公式为第五公式为根据第二无阻流量、目标气井的地层静压、井底流压和产气量,通过第四公式得到第一系数,并通过第五公式得到第二系数,将所得到的第一系数确定为目标气井对应的第一系数,将所得到的第二系数确定为目标气井对应的第二系数。
第二种方式:当第一产能方程为qg=C(pR 2-pwf 2)n,第二产能方程为qAOF=C(pR 2-0.1012)n时,对第一产能方程和第二产能方程的联立方程式求解,得到第六公式和第七公式,第六公式为第七公式为根据第二无阻流量、目标气井的地层静压、井底流压和产气量,通过第六公式得到第一系数,并通过第七公式得到第二系数,将所得到的第一系数确定为目标气井对应的第一系数,将所得到的第二系数确定为目标气井对应的第二系数。
步骤108:根据第一产能方程、目标气井对应的第一系数和第二系数,确定目标气井的合理产量。
需要说明的是,目标气井的合理产量也是目标气井的最大极限产量,如果目标气井以该合理产量进行生产,会达到稳产的效果,并且产气量较高。
具体地,步骤108的实现过程可以是:根据目标气井的地层静压、第一产能方程、目标气井对应的第一系数和第二系数,获取目标气井的多组第一试井数据,并获取目标气井的多组第二试井数据,第一试井数据包括第一井底流压和第一产气量,第二试井数据包括第二井底流压和第二产气量;根据该多组第一试井数据,确定多个第一坐标点,该多个第一坐标点中每个第一坐标点的坐标值为(第一产气量,第一井底流压);根据该多组第二试井数据,确定多个第二坐标点,该多个第二坐标点中每个第二坐标点的坐标值为(第二产气量,第二井底流压);将该多个第一坐标点连接起来,得到第四曲线;将该多个第二坐标点连接起来,得到第五曲线;将第四曲线与第五曲线的交点的横坐标确定为目标气井的合理产量。
需要说明的是,该多个第一坐标点与该多组第一试井数据一一对应,该多个第一坐标点中的某个第一坐标点的横坐标为这个第一坐标点对应的一组第一试井数据中包括的第一产气量,这个第一坐标点的纵坐标为这个第一坐标点对应的一组第一试井数据中包括的第一井底流压。
另外,该多个第二坐标点与该多组第二试井数据一一对应,该多个第二坐标点中的某个第二坐标点的横坐标为这个第二坐标点对应的一组第二试井数据中包括的第二产气量,这个第二坐标点的纵坐标为这个第二坐标点对应的一组第二试井数据中包括的第二井底流压。
再者,第四曲线为目标气井的流入动态曲线,该流入动态曲线对应天然气从地层流入目标气井的井筒中的过程。第五曲线为目标气井的流出动态曲线,该流出动态曲线对应天然气从目标气井的井筒流出至地面的过程。
其中,根据目标气井的地层静压、第一产能方程、目标气井对应的第一系数和第二系数,获取目标气井的多组第一试井数据的实现过程可以为:将目标气井的地层静压、目标气井对应的第一系数和第二系数代入第一产能方程中,根据第一产能方程确定目标气井在多个产气量中的每个产气量时的井底流压,得到多组第一试井数据。
其中,目标气井的多组第二试井数据是在目标气井中的油管尺寸和目标气井的井口压力一定的条件下,根据垂直管流法计算得到的。其中,根据垂直管流法计算目标气井的多组第二试井数据的操作可以参考相关技术,此处不做详细描述。
本发明实施例中,获取多个气井中每个气井的地层静压和多组试井数据,该多组试井数据中每组试井数据包括井底流压和产气量。然后,对于该多个气井中的每个气井,根据第一产能方程、这个气井的地层静压和多组试井数据,确定这个气井对应的第一系数和第二系数,再根据这个气井的地层静压、这个气井对应的第一系数和第二系数,通过第二产能方程得到无阻流量,将所得到的无阻流量确定为这个气井对应的第一无阻流量。然后,根据该多个气井中每个气井对应的第一无阻流量,确定目标第三系数。之后,先获取目标气井的地层静压、井底流压和产气量,再根据目标第三系数、目标气井的地层静压、井底流压和产气量,通过第三产能方程得到无阻流量,将所得到的无阻流量确定为第二无阻流量。由于在确定目标第三系数后,仅需要获取目标气井的地层静压以及一组井底流压和产气量就可以确定第二无阻流量,即目标气井的无阻流量,所以提高了确定目标气井的无阻流量的效率,降低了确定目标气井的无阻流量的成本。最后,根据第一产能方程、第二产能方程、第二无阻流量、目标气井的地层静压、井底流压和产气量,确定目标气井对应的第一系数和第二系数,再根据第一产能方程、目标气井对应的第一系数和第二系数,确定目标气井的合理产量,从而提高了确定出的气井产量的准确度。
图2是本发明实施例提供的一种确定气井产量的装置的结构示意图。参见图2,该装置包括:第一获取模块201、第一确定模块202、第二确定模块203、第三确定模块204、第二获取模块205、第四确定模块206、第五确定模块207和第六确定模块208。
第一获取模块201,用于获取多个气井中每个气井的地层静压和多组试井数据,多组试井数据中每组试井数据包括井底流压和产气量;
第一确定模块202,用于对于多个气井中的每个气井,根据第一产能方程、气井的地层静压和多组试井数据,确定气井对应的第一系数和第二系数;
第二确定模块203,用于根据气井的地层静压、气井对应的第一系数和第二系数,通过第二产能方程得到无阻流量,将所得到的无阻流量确定为气井对应的第一无阻流量;
第三确定模块204,用于根据多个气井中每个气井对应的第一无阻流量,确定目标第三系数;
第二获取模块205,用于获取目标气井的地层静压、井底流压和产气量;
第四确定模块206,用于根据目标第三系数、目标气井的地层静压、井底流压和产气量,通过第三产能方程得到无阻流量,将所得到的无阻流量确定为第二无阻流量;
第五确定模块207,用于根据第一产能方程、第二产能方程、第二无阻流量、目标气井的地层静压、井底流压和产气量,确定目标气井对应的第一系数和第二系数;
第六确定模块208,用于根据第一产能方程、目标气井对应的第一系数和第二系数,确定目标气井的合理产量;
其中,第一产能方程为pR 2-pwf 2=Aqg+Bqg 2,第二产能方程为pR 2-0.1012=AqAOF+BqAOF 2,第三产能方程为pD=αqD+(1-α)qD 2,A为第一系数,B为第二系数,α为第三系数;或者,第一产能方程为qg=C(pR 2-pwf 2)n,第二产能方程为qAOF=C(pR 2-0.1012)n,第三产能方程为lgqD=klgpD+b,C为第一系数,n为第二系数,K为第三系数中的第一子系数,b为第三系数中的第二子系数;
可选地,该第一确定模块202,包括:
第一确定子模块,用于根据气井的地层静压和多组试井数据,确定第一公式对应的第一曲线;
第二确定子模块,用于将第一曲线的斜率确定为气井对应的第二系数,将第一曲线的截距确定为气井对应的第一系数。
可选地,该第一确定模块202,包括:
第三确定子模块,用于根据气井的地层静压和多组试井数据,确定第二公式对应的第二曲线;
第四确定子模块,用于将第二曲线的斜率的倒数确定为气井对应的第二系数;
第一获取子模块,用于获取第二曲线上的任一坐标点的坐标值中的第一数值和第二数值,第一数值为pR 2-pwf 2,第二数值为qg;
第五确定子模块,用于根据第一数值、第二数值和气井对应的第二系数,通过第一产能方程得到第一系数,将所得到的第一系数确定为气井对应的第一系数。
可选地,该第三确定模块204,包括:
第六确定子模块,用于当第三产能方程为PD=αqD+(1-α)qD 2时,对于多个气井中的每个气井,根据气井对应的第一无阻流量、第一系数和第二系数,通过第三公式得到第三系数,将所得到的第三系数确定为气井对应的第三系数,第三公式为
第七确定子模块,用于将多个气井对应的第三系数的平均值确定为目标第三系数。
可选地,该第三确定模块204,包括:
第八确定子模块,用于当第三产能方程为lgqD=klgpD+b时,对于多个气井中的每个气井,根据气井的地层静压、多组试井数据和气井对应的第一无阻流量,确定气井对应的多组参数,多组参数中的每组参数包括pD和qD;
第九确定子模块,用于根据多个气井中每个气井对应的多组参数,确定多个坐标点,多个坐标点中每个坐标点的坐标值为(lgpD,lgqD);
拟合子模块,用于对多个坐标点进行线性拟合,得到第三曲线;
第十确定子模块,用于将第三曲线的斜率确定为目标第三系数中的第一子系数,将第三曲线的截距确定为目标第三系数中的第二子系数。
可选地,该第五确定模块207,包括:
第三计算子模块,用于当第一产能方程为pR 2-pwf 2=Aqg+Bqg 2,第二产能方程为pR 2-0.1012=AqAOF+BqAOF 2时,对第一产能方程和第二产能方程的联立方程式求解,得到第四公式和第五公式,第四公式为第五公式为
第十一确定子模块,用于根据第二无阻流量、目标气井的地层静压、井底流压和产气量,通过第四公式得到第一系数,并通过第五公式得到第二系数,将所得到的第一系数确定为目标气井对应的第一系数,将所得到的第二系数确定为目标气井对应的第二系数。
可选地,该第五确定模块207,包括:
第四计算子模块,用于当第一产能方程为qg=C(pR 2-pwf 2)n,第二产能方程为qAOF=C(pR 2-0.1012)n时,对第一产能方程和第二产能方程的联立方程式求解,得到第六公式和第七公式,第六公式为第七公式为
第十二确定子模块,用于根据第二无阻流量、目标气井的地层静压、井底流压和产气量,通过第六公式得到第一系数,并通过第七公式得到第二系数,将所得到的第一系数确定为目标气井对应的第一系数,将所得到的第二系数确定为目标气井对应的第二系数。
本发明实施例中,获取多个气井中每个气井的地层静压和多组试井数据,该多组试井数据中每组试井数据包括井底流压和产气量。然后,对于该多个气井中的每个气井,根据第一产能方程、这个气井的地层静压和多组试井数据,确定这个气井对应的第一系数和第二系数,再根据这个气井的地层静压、这个气井对应的第一系数和第二系数,通过第二产能方程得到无阻流量,将所得到的无阻流量确定为这个气井对应的第一无阻流量。然后,根据该多个气井中每个气井对应的第一无阻流量,确定目标第三系数。之后,先获取目标气井的地层静压、井底流压和产气量,再根据目标第三系数、目标气井的地层静压、井底流压和产气量,通过第三产能方程得到无阻流量,将所得到的无阻流量确定为第二无阻流量。由于在确定目标第三系数后,仅需要获取目标气井的地层静压以及一组井底流压和产气量就可以确定第二无阻流量,即目标气井的无阻流量,所以提高了确定目标气井的无阻流量的效率,降低了确定目标气井的无阻流量的成本。最后,根据第一产能方程、第二产能方程、第二无阻流量、目标气井的地层静压、井底流压和产气量,确定目标气井对应的第一系数和第二系数,再根据第一产能方程、目标气井对应的第一系数和第二系数,确定目标气井的合理产量,从而提高了确定出的气井产量的准确度。
关于上述实施例中的装置,其中各个模块执行操作的具体方式已经在有关方法实施例中进行了详细描述,此处将不做详细阐述说明。
需要说明的是:上述实施例提供的确定气井产量的装置在确定气井产量时,仅以上述各功能模块的划分进行举例说明,实际应用中,可以根据需要而将上述功能分配由不同的功能模块完成,即将装置的内部结构划分成不同的功能模块,以完成以上描述的全部或者部分功能。另外,上述实施例提供的确定气井产量的装置与确定气井产量的方法实施例属于同一构思,其具体实现过程详见方法实施例,这里不再赘述。
图3是本发明实施例提供的一种确定气井产量的装置的结构示意图。参见图3,该装置可以为终端300,终端300可以是:智能手机、平板电脑、MP3播放器(Moving PictureExperts Group Audio Layer III,动态影像专家压缩标准音频层面3)、MP4(MovingPicture Experts Group Audio Layer IV,动态影像专家压缩标准音频层面4)播放器、笔记本电脑或台式电脑。终端300还可能被称为用户设备、便携式终端、膝上型终端、台式终端等其他名称。
通常,终端300包括有:处理器301和存储器302。
处理器301可以包括一个或多个处理核心,比如4核心处理器、8核心处理器等。处理器301可以采用DSP(Digital Signal Processing,数字信号处理)、FPGA(Field-Programmable Gate Array,现场可编程门阵列)、PLA(Programmable Logic Array,可编程逻辑阵列)中的至少一种硬件形式来实现。处理器301也可以包括主处理器和协处理器,主处理器是用于对在唤醒状态下的数据进行处理的处理器,也称CPU(Central ProcessingUnit,中央处理器);协处理器是用于对在待机状态下的数据进行处理的低功耗处理器。在一些实施例中,处理器301可以集成有GPU(Graphics Processing Unit,图像处理器),GPU用于负责显示屏所需要显示的内容的渲染和绘制。一些实施例中,处理器301还可以包括AI(Artificial Intelligence,人工智能)处理器,该AI处理器用于处理有关机器学习的计算操作。
存储器302可以包括一个或多个计算机可读存储介质,该计算机可读存储介质可以是非暂态的。存储器302还可包括高速随机存取存储器以及非易失性存储器,比如一个或多个磁盘存储设备、闪存存储设备。在一些实施例中,存储器302中的非暂态的计算机可读存储介质用于存储至少一个指令,该至少一个指令用于被处理器301所执行以实现本申请中方法实施例提供的确定气井产量的方法。
在一些实施例中,终端300还可选包括有:外围设备接口303和至少一个外围设备。处理器301、存储器302和外围设备接口303之间可以通过总线或信号线相连。各个外围设备可以通过总线、信号线或电路板与外围设备接口303相连。具体地,外围设备包括:射频电路304、触摸显示屏305、摄像头306、音频电路307、定位组件308和电源309中的至少一种。
外围设备接口303可被用于将I/O(Input/Output,输入/输出)相关的至少一个外围设备连接到处理器301和存储器302。在一些实施例中,处理器301、存储器302和外围设备接口303被集成在同一芯片或电路板上;在一些其他实施例中,处理器301、存储器302和外围设备接口303中的任意一个或两个可以在单独的芯片或电路板上实现,本申请对此不加以限定。
射频电路304用于接收和发射RF(Radio Frequency,射频)信号,也称电磁信号。射频电路304通过电磁信号与通信网络以及其他通信设备进行通信。射频电路304将电信号转换为电磁信号进行发送,或者,将接收到的电磁信号转换为电信号。可选地,射频电路304包括:天线系统、RF收发器、一个或多个放大器、调谐器、振荡器、数字信号处理器、编解码芯片组、用户身份模块卡等。射频电路304可以通过至少一种无线通信协议来与其它终端进行通信。该无线通信协议包括但不限于:城域网、各代移动通信网络(2G、3G、4G及5G)、无线局域网和/或WiFi(Wireless Fidelity,无线保真)网络。在一些实施例中,射频电路304还可以包括NFC(Near Field Communication,近距离无线通信)有关的电路,本申请对此不加以限定。
显示屏305用于显示UI(User Interface,用户界面)。该UI可以包括图形、文本、图标、视频及其它们的任意组合。当显示屏305是触摸显示屏时,显示屏305还具有采集在显示屏305的表面或表面上方的触摸信号的能力。该触摸信号可以作为控制信号输入至处理器301进行处理。此时,显示屏305还可以用于提供虚拟按钮和/或虚拟键盘,也称软按钮和/或软键盘。在一些实施例中,显示屏305可以为一个,设置在终端300的前面板;在另一些实施例中,显示屏305可以为至少两个,分别设置在终端300的不同表面或呈折叠设计;在再一些实施例中,显示屏305可以是柔性显示屏,设置在终端300的弯曲表面上或折叠面上。甚至,显示屏305还可以设置成非矩形的不规则图形,也即异形屏。显示屏305可以采用LCD(Liquid Crystal Display,液晶显示屏)、OLED(Organic Light-Emitting Diode,有机发光二极管)等材质制备。
摄像头组件306用于采集图像或视频。可选地,摄像头组件306包括前置摄像头和后置摄像头。通常,前置摄像头设置在终端的前面板,后置摄像头设置在终端的背面。在一些实施例中,后置摄像头为至少两个,分别为主摄像头、景深摄像头、广角摄像头、长焦摄像头中的任意一种,以实现主摄像头和景深摄像头融合实现背景虚化功能、主摄像头和广角摄像头融合实现全景拍摄以及VR(Virtual Reality,虚拟现实)拍摄功能或者其它融合拍摄功能。在一些实施例中,摄像头组件306还可以包括闪光灯。闪光灯可以是单色温闪光灯,也可以是双色温闪光灯。双色温闪光灯是指暖光闪光灯和冷光闪光灯的组合,可以用于不同色温下的光线补偿。
音频电路307可以包括麦克风和扬声器。麦克风用于采集用户及环境的声波,并将声波转换为电信号输入至处理器301进行处理,或者输入至射频电路304以实现语音通信。出于立体声采集或降噪的目的,麦克风可以为多个,分别设置在终端300的不同部位。麦克风还可以是阵列麦克风或全向采集型麦克风。扬声器则用于将来自处理器301或射频电路304的电信号转换为声波。扬声器可以是传统的薄膜扬声器,也可以是压电陶瓷扬声器。当扬声器是压电陶瓷扬声器时,不仅可以将电信号转换为人类可听见的声波,也可以将电信号转换为人类听不见的声波以进行测距等用途。在一些实施例中,音频电路307还可以包括耳机插孔。
定位组件308用于定位终端300的当前地理位置,以实现导航或LBS(LocationBased Service,基于位置的服务)。定位组件308可以是基于美国的GPS(GlobalPositioning System,全球定位系统)、中国的北斗系统、俄罗斯的格雷纳斯系统或欧盟的伽利略系统的定位组件。
电源309用于为终端300中的各个组件进行供电。电源309可以是交流电、直流电、一次性电池或可充电电池。当电源309包括可充电电池时,该可充电电池可以支持有线充电或无线充电。该可充电电池还可以用于支持快充技术。
在一些实施例中,终端300还包括有一个或多个传感器310。该一个或多个传感器310包括但不限于:加速度传感器311、陀螺仪传感器312、压力传感器313、指纹传感器314、光学传感器315以及接近传感器316。
加速度传感器311可以检测以终端300建立的坐标系的三个坐标轴上的加速度大小。比如,加速度传感器311可以用于检测重力加速度在三个坐标轴上的分量。处理器301可以根据加速度传感器311采集的重力加速度信号,控制触摸显示屏305以横向视图或纵向视图进行用户界面的显示。加速度传感器311还可以用于游戏或用户的运动数据的采集。
陀螺仪传感器312可以检测终端300的机体方向及转动角度,陀螺仪传感器312可以与加速度传感器311协同采集用户对终端300的3D动作。处理器301根据陀螺仪传感器312采集的数据,可以实现如下功能:动作感应(比如根据用户的倾斜操作来改变UI)、拍摄时的图像稳定、游戏控制以及惯性导航。
压力传感器313可以设置在终端300的侧边框和/或触摸显示屏305的下层。当压力传感器313设置在终端300的侧边框时,可以检测用户对终端300的握持信号,由处理器301根据压力传感器313采集的握持信号进行左右手识别或快捷操作。当压力传感器313设置在触摸显示屏305的下层时,由处理器301根据用户对触摸显示屏305的压力操作,实现对UI界面上的可操作性控件进行控制。可操作性控件包括按钮控件、滚动条控件、图标控件、菜单控件中的至少一种。
指纹传感器314用于采集用户的指纹,由处理器301根据指纹传感器314采集到的指纹识别用户的身份,或者,由指纹传感器314根据采集到的指纹识别用户的身份。在识别出用户的身份为可信身份时,由处理器301授权该用户执行相关的敏感操作,该敏感操作包括解锁屏幕、查看加密信息、下载软件、支付及更改设置等。指纹传感器314可以被设置在终端300的正面、背面或侧面。当终端300上设置有物理按键或厂商Logo时,指纹传感器314可以与物理按键或厂商Logo集成在一起。
光学传感器315用于采集环境光强度。在一个实施例中,处理器301可以根据光学传感器315采集的环境光强度,控制触摸显示屏305的显示亮度。具体地,当环境光强度较高时,调高触摸显示屏305的显示亮度;当环境光强度较低时,调低触摸显示屏305的显示亮度。在另一个实施例中,处理器301还可以根据光学传感器315采集的环境光强度,动态调整摄像头组件306的拍摄参数。
接近传感器316也称距离传感器,通常设置在终端300的前面板。接近传感器316用于采集用户与终端300的正面之间的距离。在一个实施例中,当接近传感器316检测到用户与终端300的正面之间的距离逐渐变小时,由处理器301控制触摸显示屏305从亮屏状态切换为息屏状态;当接近传感器316检测到用户与终端300的正面之间的距离逐渐变大时,由处理器301控制触摸显示屏305从息屏状态切换为亮屏状态。
本领域技术人员可以理解,图3中示出的结构并不构成对终端300的限定,可以包括比图示更多或更少的组件,或者组合某些组件,或者采用不同的组件布置。
本领域技术人员在考虑说明书及实践这里公开的发明后,将容易想到本发明的其它实施方案。本申请旨在涵盖本发明的任何变型、用途或者适应性变化,这些变型、用途或者适应性变化遵循本发明的一般性原理并包括本发明未公开的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的,本发明的真正范围和精神由下面的权利要求指出。
应当理解的是,本发明并不局限于上面已经描述并在附图中示出的精确结构,并且可以在不脱离其范围进行各种修改和改变。本发明的范围仅由所附的权利要求来限制。
Claims (16)
1.一种确定气井产量的方法,其特征在于,所述方法包括:
获取多个气井中每个气井的地层静压和多组试井数据,所述多组试井数据中每组试井数据包括井底流压和产气量;
对于所述多个气井中的每个气井,根据第一产能方程、所述气井的地层静压和多组试井数据,确定所述气井对应的第一系数和第二系数;
根据所述气井的地层静压、所述气井对应的第一系数和第二系数,通过第二产能方程得到无阻流量,将所得到的无阻流量确定为所述气井对应的第一无阻流量;
根据所述多个气井中每个气井对应的第一无阻流量,确定目标第三系数;
获取目标气井的地层静压、井底流压和产气量;
根据所述目标第三系数、所述目标气井的地层静压、井底流压和产气量,通过第三产能方程得到无阻流量,将所得到的无阻流量确定为第二无阻流量;
根据所述第一产能方程、所述第二产能方程、所述第二无阻流量、所述目标气井的地层静压、井底流压和产气量,确定所述目标气井对应的第一系数和第二系数;
根据所述第一产能方程、所述目标气井对应的第一系数和第二系数,确定所述目标气井的合理产量;
其中,所述第一产能方程为pR 2-pwf 2=Aqg+Bqg 2,所述第二产能方程为pR 2-0.1012=AqAOF+BqAOF 2,所述第三产能方程为pD=αqD+(1-α)qD 2,所述A为第一系数,所述B为第二系数,所述α为第三系数;或者,所述第一产能方程为qg=C(pR 2-pwf 2)n,所述第二产能方程为qAOF=C(pR 2-0.1012)n,所述第三产能方程为lgqD=klgpD+b,所述C为第一系数,所述n为第二系数,所述K为第三系数中的第一子系数,所述b为第三系数中的第二子系数;
3.如权利要求1所述的方法,其特征在于,所述根据第一产能方程、所述气井的地层静压和多组试井数据,确定所述气井对应的第一系数和第二系数,包括:
根据所述气井的地层静压和多组试井数据,确定所述第二公式对应的第二曲线;
将所述第二曲线的斜率的倒数确定为所述气井对应的第二系数;
获取所述第二曲线上的任一坐标点的坐标值中的第一数值和第二数值,所述第一数值为pR 2-pwf 2,所述第二数值为qg;
根据所述第一数值、所述第二数值和所述气井对应的第二系数,通过所述第一产能方程得到第一系数,将所得到的第一系数确定为所述气井对应的第一系数。
4.如权利要求1所述的方法,其特征在于,所述根据所述多个气井中每个气井对应的第一无阻流量,确定目标第三系数,包括:
当所述第三产能方程为PD=αqD+(1-α)qD 2时,对于所述多个气井中的每个气井,根据所述气井对应的第一无阻流量、第一系数和第二系数,通过第三公式得到第三系数,将所得到的第三系数确定为所述气井对应的第三系数,所述第三公式为
将所述多个气井对应的第三系数的平均值确定为目标第三系数。
5.如权利要求1所述的方法,其特征在于,所述根据所述多个气井中每个气井对应的第一无阻流量,确定目标第三系数,包括:
当所述第三产能方程为lgqD=klgpD+b时,对于所述多个气井中的每个气井,根据所述气井的地层静压、多组试井数据和所述气井对应的第一无阻流量,确定所述气井对应的多组参数,所述多组参数中的每组参数包括pD和qD;
根据所述多个气井中每个气井对应的多组参数,确定多个坐标点,所述多个坐标点中每个坐标点的坐标值为(lgpD,lgqD);
对所述多个坐标点进行线性拟合,得到第三曲线;
将所述第三曲线的斜率确定为目标第三系数中的第一子系数,将所述第三曲线的截距确定为目标第三系数中的第二子系数。
6.如权利要求1所述的方法,其特征在于,所述根据所述第一产能方程、所述第二产能方程、所述第二无阻流量、所述目标气井的地层静压、井底流压和产气量,确定所述目标气井对应的第一系数和第二系数,包括:
当所述第一产能方程为pR 2-pwf 2=Aqg+Bqg 2,所述第二产能方程为pR 2-0.1012=AqAOF+BqAOF 2时,对所述第一产能方程和所述第二产能方程的联立方程式求解,得到第四公式和第五公式,所述第四公式为所述第五公式为
根据所述第二无阻流量、所述目标气井的地层静压、井底流压和产气量,通过所述第四公式得到第一系数,并通过所述第五公式得到第二系数,将所得到的第一系数确定为所述目标气井对应的第一系数,将所得到的第二系数确定为所述目标气井对应的第二系数。
7.如权利要求1所述的方法,其特征在于,所述根据所述第一产能方程、所述第二产能方程、所述第二无阻流量、所述目标气井的地层静压、井底流压和产气量,确定所述目标气井对应的第一系数和第二系数,包括:
当所述第一产能方程为qg=C(pR 2-pwf 2)n,所述第二产能方程为qAOF=C(pR 2-0.1012)n时,对所述第一产能方程和所述第二产能方程的联立方程式求解,得到第六公式和第七公式,所述第六公式为所述第七公式为
根据所述第二无阻流量、所述目标气井的地层静压、井底流压和产气量,通过所述第六公式得到第一系数,并通过所述第七公式得到第二系数,将所得到的第一系数确定为所述目标气井对应的第一系数,将所得到的第二系数确定为所述目标气井对应的第二系数。
8.一种确定气井产量的装置,其特征在于,所述装置包括:
第一获取模块,用于获取多个气井中每个气井的地层静压和多组试井数据,所述多组试井数据中每组试井数据包括井底流压和产气量;
第一确定模块,用于对于所述多个气井中的每个气井,根据第一产能方程、所述气井的地层静压和多组试井数据,确定所述气井对应的第一系数和第二系数;
第二确定模块,用于根据所述气井的地层静压、所述气井对应的第一系数和第二系数,通过第二产能方程得到无阻流量,将所得到的无阻流量确定为所述气井对应的第一无阻流量;
第三确定模块,用于根据所述多个气井中每个气井对应的第一无阻流量,确定目标第三系数;
第二获取模块,用于获取目标气井的地层静压、井底流压和产气量;
第四确定模块,用于根据所述目标第三系数、所述目标气井的地层静压、井底流压和产气量,通过第三产能方程得到无阻流量,将所得到的无阻流量确定为第二无阻流量;
第五确定模块,用于根据所述第一产能方程、所述第二产能方程、所述第二无阻流量、所述目标气井的地层静压、井底流压和产气量,确定所述目标气井对应的第一系数和第二系数;
第六确定模块,用于根据所述第一产能方程、所述目标气井对应的第一系数和第二系数,确定所述目标气井的合理产量;
其中,所述第一产能方程为pR 2-pwf 2=Aqg+Bqg 2,所述第二产能方程为pR 2-0.1012=AqAOF+BqAOF 2,所述第三产能方程为pD=αqD+(1-α)qD 2,所述A为第一系数,所述B为第二系数,所述α为第三系数;或者,所述第一产能方程为qg=C(pR 2-pwf 2)n,所述第二产能方程为qAOF=C(pR 2-0.1012)n,所述第三产能方程为lgqD=klgpD+b,所述C为第一系数,所述n为第二系数,所述K为第三系数中的第一子系数,所述b为第三系数中的第二子系数;
10.如权利要求8所述的装置,其特征在于,所述第一确定模块,包括:
第二计算子模块,用于当所述第一产能方程为qg=C(pR 2-pwf 2)n时,将所述第一产能方程的两端同取对数,得到第二公式,所述第二公式为所述y2=lg(pR 2-pwf 2),所述x2=lgqg;
第三确定子模块,用于根据所述气井的地层静压和多组试井数据,确定所述第二公式对应的第二曲线;
第四确定子模块,用于将所述第二曲线的斜率的倒数确定为所述气井对应的第二系数;
第一获取子模块,用于获取所述第二曲线上的任一坐标点的坐标值中的第一数值和第二数值,所述第一数值为pR 2-pwf 2,所述第二数值为qg;
第五确定子模块,用于根据所述第一数值、所述第二数值和所述气井对应的第二系数,通过所述第一产能方程得到第一系数,将所得到的第一系数确定为所述气井对应的第一系数。
12.如权利要求8所述的装置,其特征在于,所述第三确定模块,包括:
第八确定子模块,用于当所述第三产能方程为lgqD=klgpD+b时,对于所述多个气井中的每个气井,根据所述气井的地层静压、多组试井数据和所述气井对应的第一无阻流量,确定所述气井对应的多组参数,所述多组参数中的每组参数包括pD和qD;
第九确定子模块,用于根据所述多个气井中每个气井对应的多组参数,确定多个坐标点,所述多个坐标点中每个坐标点的坐标值为(lgpD,lgqD);
拟合子模块,用于对所述多个坐标点进行线性拟合,得到第三曲线;
第十确定子模块,用于将所述第三曲线的斜率确定为目标第三系数中的第一子系数,将所述第三曲线的截距确定为目标第三系数中的第二子系数。
15.一种确定气井产量的装置,其特征在于,所述装置包括:
处理器;
用于存储处理器可执行指令的存储器;
其中,所述处理器被配置为执行权利要求1-7所述的任一项方法的步骤。
16.一种计算机可读存储介质,所述计算机可读存储介质上存储有指令,其特征在于,所述指令被处理器执行时实现权利要求1-7所述的任一项方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811168646.5A CN111005715B (zh) | 2018-10-08 | 2018-10-08 | 确定气井产量的方法、装置及存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811168646.5A CN111005715B (zh) | 2018-10-08 | 2018-10-08 | 确定气井产量的方法、装置及存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111005715A true CN111005715A (zh) | 2020-04-14 |
CN111005715B CN111005715B (zh) | 2023-02-28 |
Family
ID=70110713
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811168646.5A Active CN111005715B (zh) | 2018-10-08 | 2018-10-08 | 确定气井产量的方法、装置及存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111005715B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113047815A (zh) * | 2021-04-02 | 2021-06-29 | 西南石油大学 | 一种气井合理配产方法 |
CN113931621A (zh) * | 2020-07-14 | 2022-01-14 | 中国石油天然气股份有限公司 | 气井积液信息的确定方法、装置及存储介质 |
CN114060017A (zh) * | 2020-08-04 | 2022-02-18 | 中国石油天然气股份有限公司 | 一种低渗透油藏直井体积改造方法及改造系统 |
CN114382465A (zh) * | 2020-10-21 | 2022-04-22 | 中国石油天然气股份有限公司 | 确定气井的稳产年限的方法、装置、终端和存储介质 |
CN114517666A (zh) * | 2020-11-20 | 2022-05-20 | 中国石油天然气股份有限公司 | 气井的节流方法、装置、计算机设备及存储介质 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103046925A (zh) * | 2012-12-31 | 2013-04-17 | 中国石油天然气股份有限公司 | 基于二项式的凝析气藏绝对无阻流量的获取方法及系统 |
CN103046924A (zh) * | 2012-12-31 | 2013-04-17 | 中国石油天然气股份有限公司 | 基于指数式的凝析气藏绝对无阻流量的获取方法及系统 |
CN104504457A (zh) * | 2014-12-04 | 2015-04-08 | 中国石油大港油田勘探开发研究院 | 产水气井产能预测方法 |
CN104832166A (zh) * | 2015-03-20 | 2015-08-12 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | 一种页岩气水平井初期产能预测方法 |
CN104832156A (zh) * | 2015-05-05 | 2015-08-12 | 中国石油天然气股份有限公司 | 一种预估气井产量的方法 |
CN106869917A (zh) * | 2017-03-24 | 2017-06-20 | 中国石油集团川庆钻探工程有限公司 | 基于井底和井口流入动态曲线共同约束的产能评价方法 |
CN107301306A (zh) * | 2017-08-04 | 2017-10-27 | 中国石油化工股份有限公司华北油气分公司 | 用于致密砂岩气藏压裂水平井的动态无阻流量预测方法 |
CN108344853A (zh) * | 2018-01-24 | 2018-07-31 | 西南石油大学 | 平面非均质性定容干气气藏直井绝对无阻流量测试方法 |
CN108612525A (zh) * | 2018-04-19 | 2018-10-02 | 重庆科技学院 | 一种气藏动态储量计算方法 |
-
2018
- 2018-10-08 CN CN201811168646.5A patent/CN111005715B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103046925A (zh) * | 2012-12-31 | 2013-04-17 | 中国石油天然气股份有限公司 | 基于二项式的凝析气藏绝对无阻流量的获取方法及系统 |
CN103046924A (zh) * | 2012-12-31 | 2013-04-17 | 中国石油天然气股份有限公司 | 基于指数式的凝析气藏绝对无阻流量的获取方法及系统 |
CN104504457A (zh) * | 2014-12-04 | 2015-04-08 | 中国石油大港油田勘探开发研究院 | 产水气井产能预测方法 |
CN104832166A (zh) * | 2015-03-20 | 2015-08-12 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | 一种页岩气水平井初期产能预测方法 |
CN104832156A (zh) * | 2015-05-05 | 2015-08-12 | 中国石油天然气股份有限公司 | 一种预估气井产量的方法 |
CN106869917A (zh) * | 2017-03-24 | 2017-06-20 | 中国石油集团川庆钻探工程有限公司 | 基于井底和井口流入动态曲线共同约束的产能评价方法 |
CN107301306A (zh) * | 2017-08-04 | 2017-10-27 | 中国石油化工股份有限公司华北油气分公司 | 用于致密砂岩气藏压裂水平井的动态无阻流量预测方法 |
CN108344853A (zh) * | 2018-01-24 | 2018-07-31 | 西南石油大学 | 平面非均质性定容干气气藏直井绝对无阻流量测试方法 |
CN108612525A (zh) * | 2018-04-19 | 2018-10-02 | 重庆科技学院 | 一种气藏动态储量计算方法 |
Non-Patent Citations (2)
Title |
---|
杨元明等: "应用页岩气井地面测试数据分析气井产能", 《油气井测试》 * |
胡俊坤: "利用生产动态资料确定气井产能方程新方法", 《天然气地球科学》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113931621A (zh) * | 2020-07-14 | 2022-01-14 | 中国石油天然气股份有限公司 | 气井积液信息的确定方法、装置及存储介质 |
CN113931621B (zh) * | 2020-07-14 | 2023-08-22 | 中国石油天然气股份有限公司 | 气井积液信息的确定方法、装置及存储介质 |
CN114060017A (zh) * | 2020-08-04 | 2022-02-18 | 中国石油天然气股份有限公司 | 一种低渗透油藏直井体积改造方法及改造系统 |
CN114060017B (zh) * | 2020-08-04 | 2023-08-22 | 中国石油天然气股份有限公司 | 一种低渗透油藏直井体积改造方法及改造系统 |
CN114382465A (zh) * | 2020-10-21 | 2022-04-22 | 中国石油天然气股份有限公司 | 确定气井的稳产年限的方法、装置、终端和存储介质 |
CN114382465B (zh) * | 2020-10-21 | 2024-04-30 | 中国石油天然气股份有限公司 | 确定气井的稳产年限的方法、装置、终端和存储介质 |
CN114517666A (zh) * | 2020-11-20 | 2022-05-20 | 中国石油天然气股份有限公司 | 气井的节流方法、装置、计算机设备及存储介质 |
CN114517666B (zh) * | 2020-11-20 | 2023-08-22 | 中国石油天然气股份有限公司 | 气井的节流方法、装置、计算机设备及存储介质 |
CN113047815A (zh) * | 2021-04-02 | 2021-06-29 | 西南石油大学 | 一种气井合理配产方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111005715B (zh) | 2023-02-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107885533B (zh) | 管理组件代码的方法及装置 | |
CN111005715B (zh) | 确定气井产量的方法、装置及存储介质 | |
CN108965922B (zh) | 视频封面生成方法、装置和存储介质 | |
CN109102811B (zh) | 音频指纹的生成方法、装置及存储介质 | |
CN108764530B (zh) | 油井抽油机的工作参数配置方法及装置 | |
CN109281648B (zh) | 确定油藏的合理井网密度的方法和装置 | |
CN111255434B (zh) | 气井的试井方法、装置和计算机存储介质 | |
CN108733902B (zh) | 油藏渗透率的确定方法、装置及存储介质 | |
CN109281664B (zh) | 水平井产油量的预测方法和装置 | |
CN110865170B (zh) | 确定含砂量的方法、装置及存储介质 | |
CN111287718B (zh) | 确定射孔参数的方法、装置及存储介质 | |
CN109899051B (zh) | 油井设备的评价标准确定方法、装置及存储介质 | |
CN111063372B (zh) | 确定音高特征的方法、装置、设备及存储介质 | |
CN109388732B (zh) | 音乐地图的生成和显示方法、装置及存储介质 | |
CN109902844B (zh) | 注水系统的优化信息确定方法、装置及存储介质 | |
CN110580561B (zh) | 一种油井增油效果的分析方法、装置及存储介质 | |
CN113418848B (zh) | 岩心的渗透参数的确定方法、装置、设备及存储介质 | |
CN115263244B (zh) | 气水层的排水控制方法、装置及计算机存储介质 | |
CN112182986B (zh) | 确定面积波及信息的方法、装置、设备及存储介质 | |
CN111088969B (zh) | 注水井的分注方案确定方法、装置和存储介质 | |
CN110442975B (zh) | 发动机的温度场确定方法、装置、终端及存储介质 | |
CN110058306B (zh) | 确定三维速度体的方法、装置及计算机存储介质 | |
CN110633335B (zh) | 获取poi数据的方法、终端和可读存储介质 | |
CN110888992A (zh) | 多媒体数据处理方法、装置、计算机设备和可读存储介质 | |
CN111860543A (zh) | 图像处理方法、装置、设备及存储介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |