CN117784278B - Prediction method and prediction system for dense sandstone gas dessert - Google Patents
Prediction method and prediction system for dense sandstone gas dessert Download PDFInfo
- Publication number
- CN117784278B CN117784278B CN202410206324.4A CN202410206324A CN117784278B CN 117784278 B CN117784278 B CN 117784278B CN 202410206324 A CN202410206324 A CN 202410206324A CN 117784278 B CN117784278 B CN 117784278B
- Authority
- CN
- China
- Prior art keywords
- gas
- decreasing
- yield
- cumulative
- intersection
- 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.)
- Active
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
Description
技术领域Technical Field
本发明涉及油气地质勘探技术领域,特别涉及一种致密砂岩气甜点的预测方法和预测系统。The invention relates to the technical field of oil and gas geological exploration, and in particular to a prediction method and a prediction system for tight sandstone gas sweet spots.
背景技术Background technique
随着非常规油气地质理论的不断发展与勘探开发技术的进步,致密砂岩气甜点预测成为非常规天然气地质学中的关键和核心,也是天然气成藏、富集与目标预测研究的薄弱环节。With the continuous development of unconventional oil and gas geological theory and the advancement of exploration and development technology, the prediction of tight sandstone gas sweet spots has become the key and core of unconventional natural gas geology, and is also the weak link in the research of natural gas accumulation, enrichment and target prediction.
国内多数气田中致密砂岩气勘探开发已进入“开发为主、勘探为辅”的局面,评价致密砂岩气甜点主要控制因素由六大要素逐渐过度到储层为主,现有技术中利用富集主控因素、多元线性回归、储层双甜点评价参数、基于机器学习的地震属性预测等技术计算致密砂岩气井无阻流量,从而评价致密砂岩气甜点。The exploration and development of tight sandstone gas in most gas fields in China has entered a situation of "development as the main factor and exploration as the auxiliary factor". The main controlling factors for evaluating tight sandstone gas sweet spots have gradually shifted from six major factors to reservoir-based. The existing technology uses enrichment main controlling factors, multivariate linear regression, reservoir double sweet spot evaluation parameters, and seismic attribute prediction based on machine learning to calculate the unimpeded flow rate of tight sandstone gas wells, thereby evaluating tight sandstone gas sweet spots.
上述现有技术存在的缺陷是:致密砂岩气储层具有“低孔渗-致密、开发时须经压裂改造”的特点,通常意义上的致密砂岩气井无阻流量不能反映出致密砂岩的富集程度,对致密砂岩气甜点评价不够精确。The defects of the above-mentioned existing technologies are: the tight sandstone gas reservoir has the characteristics of "low porosity and permeability - dense, and must be fracturing and transformed during development", the unimpeded flow rate of tight sandstone gas wells in the usual sense cannot reflect the enrichment degree of tight sandstone, and the evaluation of tight sandstone gas sweet spots is not accurate enough.
发明内容Summary of the invention
基于此,有必要针对上述技术问题,提供一种致密砂岩气甜点的预测方法和预测系统。Based on this, it is necessary to provide a prediction method and prediction system for tight sandstone gas sweet spots in response to the above technical problems.
本发明实施例提供一种致密砂岩气甜点的预测方法,包括:The embodiment of the present invention provides a method for predicting tight sandstone gas sweet spots, comprising:
绘制致密砂岩气井的累产气量、月产气量与生产年月的交会图,Draw a cross-plot of the cumulative gas production, monthly gas production and production year and month of tight sandstone gas wells.
根据交会图判识单井累产气量的递减趋势,通过递减趋势计算出致密砂岩气的单井累产气量;According to the intersection diagram, the decreasing trend of the cumulative gas production of a single well is identified, and the cumulative gas production of a single well of tight sandstone gas is calculated through the decreasing trend.
分别绘制单井累产气量与孔隙度Φ、单层砂体厚度H、富集系数EF的预测交会图以及单层砂体厚度H与含气饱和度S g 的预测交会图;Draw the predicted cross-plots of the cumulative gas production of a single well and porosity Φ , single-layer sand body thickness H , and enrichment factor EF , as well as the predicted cross-plots of the single-layer sand body thickness H and gas saturation Sg ;
分别绘制单井累产气量与孔隙度Φ、含气饱和度S g 和单层砂体厚度H、富集系数EF的预测交会图;Draw the predicted cross-plots of the cumulative gas production of a single well, porosity Φ , gas saturation Sg , single-layer sand body thickness H , and enrichment factor EF respectively;
根据所有的预测交会图得到致密砂岩气富集度的下限值,所述下限值为致密砂岩气甜点值。The lower limit value of the tight sandstone gas enrichment is obtained according to all the predicted cross-plots, and the lower limit value is the tight sandstone gas sweet spot value.
另外的,所述根据交会图判识单井累产气量的递减趋势,包括:In addition, the method of identifying the decreasing trend of the cumulative gas production of a single well according to the cross-plot includes:
致密砂岩气井的累产气量与月产气量在交会图上呈直线关系时,单井累产气量的递减趋势为指数递减;When the cumulative gas production and monthly gas production of tight sandstone gas wells present a linear relationship on the cross-plot, the decreasing trend of the cumulative gas production of a single well is exponential decrease;
致密砂岩气井的累产气量与月产气量在交会图上呈半对数直线关系时,单井累产气量的递减趋势为调和递减;When the cumulative gas production and monthly gas production of tight sandstone gas wells present a semi-logarithmic linear relationship on the cross-plot, the decreasing trend of the cumulative gas production of a single well is a harmonic decrease;
致密砂岩气井的累产气量与月产气量在交会图上既不呈直线关系也不呈半对数直线关系时,单井累产气量的递减趋势为双曲递减。When the cumulative gas production and monthly gas production of tight sandstone gas wells are neither in a linear relationship nor in a semi-logarithmic linear relationship on the intersection diagram, the decreasing trend of the cumulative gas production of a single well is a hyperbolic decrease.
另外的,所述通过递减趋势预测出致密砂岩气的单井累产气量包括:In addition, the single well cumulative gas production of tight sandstone gas predicted by the decreasing trend includes:
当递减趋势为调和递减,递减指数n=1,则有:When the decreasing trend is harmonic decreasing and the decreasing index n=1, then:
Q / Qi = (D / Di )n;Q / Qi = (D / Di ) n ;
当递减趋势为指数递减,递减指数,则有:When the decreasing trend is exponential, the decreasing exponential , then:
Q / Qi = (D / Di )n Q / Qi = (D / Di ) n
其中,Q为单井累产气量,Qi为递减阶段的初始产量,D为产量递减率,Di为初始瞬时递减率,n为递减指数,i为正整数;Wherein, Q is the cumulative gas production of a single well, Qi is the initial production in the decline stage, D is the production decline rate, Di is the initial instantaneous decline rate, n is the decline index, and i is a positive integer;
当递减趋势为双曲递减,1<递减指数n<∞,When the decreasing trend is hyperbolic, 1<decreasing exponent n<∞,
Q = Qi / {1 + (Di / n)×t}n Q = Qi /{1+( Di /n)×t} n
其中,Q为单井累产气量,Qi为递减阶段的初始产量,Di为初始瞬时递减率,t为递减阶段生产时间,n为递减指数,i为正整数。Among them, Q is the cumulative gas production of a single well, Qi is the initial production in the decline stage, Di is the initial instantaneous decline rate, t is the production time in the decline stage, n is the decline index, and i is a positive integer.
另外的,所述静态控制因素具体包括:In addition, the static control factors specifically include:
孔隙度Φ:利用测井组合参数自然伽马GR、中子测井CNL、声波测井AC和密度测井DEN,计算单层砂体的孔隙度Φ;Porosity Φ: The porosity Φ of a single sand body is calculated using the combined logging parameters of natural gamma GR, neutron logging CNL, acoustic logging AC and density logging DEN;
含气饱和度S g :通过Archie公式计算含水饱和度模型,并求取含气饱和度S g :Gas saturation S g : Calculate the water saturation model using Archie's formula and obtain the gas saturation S g :
Sw n = abRw/RtФm S w n = abR w /R t Ф m
Sg = 1 - Sw Sg = 1 - Sw
其中,Sw为含水饱和度,a、b为岩性系数,RW为地层水电阻率,Rt为地层电阻率,Ф为孔隙度,m为胶结系数,n为饱和度指数;Among them, Sw is water saturation, a and b are lithology coefficients, RW is formation water resistivity, Rt is formation resistivity, Ф is porosity, m is cementation coefficient, and n is saturation index;
单层砂体厚度H:根据单层砂体的孔隙度Φ和含气饱和度S g 确定单层砂体厚度。Single-layer sand body thickness H: The thickness of a single-layer sand body is determined based on the porosity Φ and gas saturation Sg of the single-layer sand body.
另外的,所述富集系数EF为:In addition, the enrichment factor EF is:
EF = H × Φ × S g EF = H × Φ × S g
其中,H为单层砂体厚度,Φ为孔隙度,S g 为含气饱和度。Among them, H is the thickness of a single sand layer, Φ is the porosity, and Sg is the gas saturation.
另外的,所述得到致密砂岩气富集的下限值包括:In addition, the lower limit values for obtaining tight sandstone gas enrichment include:
通过单井累产气量与单层砂体厚度H交会,建立致密砂岩气藏富集度单层砂体厚度评价标准;By intersecting the cumulative gas production of a single well with the thickness H of a single sand body, an evaluation standard for the thickness of a single sand body in terms of the richness of tight sandstone gas reservoirs is established.
通过单井累产气量与孔隙度Φ交会,建立致密砂岩气富集度孔隙度评价标准;Through the intersection of single well cumulative gas production and porosity Φ, the evaluation standard of tight sandstone gas enrichment and porosity is established;
通过单层砂体厚度H与含气饱和度S g 交会,建立致密砂岩气富集度含气饱和度评价标准;Through the intersection of single-layer sand body thickness H and gas saturation Sg , the evaluation standard of tight sandstone gas richness and gas saturation is established;
通过单井累产气量与富集系数EF交会,建立致密砂岩气富集度富集系数评价标准;Through the intersection of single well cumulative gas production and enrichment factor EF, the evaluation standard of tight sandstone gas enrichment and enrichment factor is established;
根据单层砂体厚度评价标准,孔隙度评价标准,含气饱和度评价标准和富集系数评价标准,确定致密砂岩气富集度的下限值。The lower limit of tight sandstone gas enrichment is determined based on the single-layer sand body thickness evaluation standard, porosity evaluation standard, gas saturation evaluation standard and enrichment coefficient evaluation standard.
另外的,一种致密砂岩气甜点的预测系统,包括:In addition, a prediction system for tight sandstone gas sweet spots includes:
绘制模块,用于绘制致密砂岩气井的累产气量、月产气量与生产年月的交会图;The drawing module is used to draw the intersection diagram of the cumulative gas production, monthly gas production and production year and month of tight sandstone gas wells;
判识模块,用于根据交会图判识单井累产气量的递减趋势,通过递减趋势计算出致密砂岩气的单井累产气量;The identification module is used to identify the decreasing trend of the cumulative gas production of a single well according to the intersection diagram, and calculate the cumulative gas production of a single well of tight sandstone gas through the decreasing trend;
联立模块,用于获取致密砂岩气的静态控制因素,所述静态控制因素包括:孔隙度Φ、含气饱和度S g 和单层砂体厚度H,联立所有的静态控制因素创建富集系数EF;A combined module is used to obtain the static control factors of tight sandstone gas, which include porosity Φ , gas saturation Sg and single-layer sand body thickness H. All static control factors are combined to create an enrichment factor EF ;
二次绘制模块,用于分别绘制单井累产气量与孔隙度Φ、单层砂体厚度H、富集系数EF的预测交会图以及单层砂体厚度H与含气饱和度S g 的预测交会图;The secondary drawing module is used to draw the predicted cross-plot of the cumulative gas production of a single well and the porosity Φ , the thickness of a single sand body H , and the enrichment factor EF , as well as the predicted cross-plot of the thickness of a single sand body H and the gas saturation Sg ;
推断模块,用于根据所有的预测交会图得到致密砂岩气富集度的下限值,所述下限值为致密砂岩气甜点值。The inference module is used to obtain the lower limit value of the tight sandstone gas enrichment according to all the predicted intersection diagrams, and the lower limit value is the tight sandstone gas sweet spot value.
本发明实施例提供的上述一种致密砂岩气甜点的预测方法和预测系统,与现有技术相比,其有益效果如下:Compared with the prior art, the above-mentioned prediction method and prediction system for tight sandstone gas sweet spots provided by the embodiment of the present invention have the following beneficial effects:
本发明通过绘制致密砂岩气井的累产气量、月产气量与生产年月的交会图,根据交会图判识单井累产气量的递减趋势;通过递减趋势计算出致密砂岩气的单井累产气量;获取致密砂岩气的静态控制因素,所述静态控制因素包括:孔隙度Φ、含气饱和度S g 和单层砂体厚度H,联立所有的静态控制因素创建富集系数EF;分别绘制单井累产气量与孔隙度Φ、单层砂体厚度H、富集系数EF的预测交会图以及单层砂体厚度H与含气饱和度S g 的预测交会图;根据预测交会图得到致密砂岩气富集的下限值,下限值即为致密砂岩气甜点值。相比于现有技术,从致密砂岩气储层的特点出发,分别绘制单井累产气量与孔隙度Φ、单层砂体厚度H、富集系数EF的预测交会图以及单层砂体厚度H与含气饱和度S g 的预测交会图;通过预测交会图得到致密砂岩气富集的下限值,下限值即为致密砂岩气甜点值。The present invention draws a cross-plot of the cumulative gas production, monthly gas production and production year and month of a tight sandstone gas well, and identifies the decreasing trend of the cumulative gas production of a single well according to the cross-plot; calculates the cumulative gas production of a single well of tight sandstone gas according to the decreasing trend; obtains static control factors of tight sandstone gas, the static control factors include: porosity Φ , gas saturation Sg and single-layer sand body thickness H , and creates an enrichment factor EF by combining all static control factors; draws predicted cross-plots of the cumulative gas production of a single well with porosity Φ , single-layer sand body thickness H , and enrichment factor EF, and predicted cross-plots of the single-layer sand body thickness H and gas saturation Sg ; and obtains a lower limit value of the enrichment of tight sandstone gas according to the predicted cross-plot, which is the sweet spot value of tight sandstone gas. Compared with the existing technology, based on the characteristics of tight sandstone gas reservoirs, the predicted intersection diagrams of single well cumulative gas production and porosity Φ , single-layer sand body thickness H , and enrichment factor EF , as well as the predicted intersection diagram of single-layer sand body thickness H and gas saturation Sg were drawn respectively ; the lower limit value of tight sandstone gas enrichment was obtained through the predicted intersection diagram, and the lower limit value was the tight sandstone gas sweet spot value.
上述方案能够替代致密砂岩气井无阻流量来反映致密砂岩的富集程度,对致密砂岩气甜点的探测更加准确。The above scheme can replace the open flow rate of tight sandstone gas wells to reflect the enrichment degree of tight sandstone, and the detection of tight sandstone gas sweet spots is more accurate.
附图说明BRIEF DESCRIPTION OF THE DRAWINGS
图1为一个实施例中提供的一种致密砂岩气甜点的预测方法的流程示意图;FIG1 is a schematic flow chart of a method for predicting tight sandstone gas sweet spots provided in one embodiment;
图2为一个实施例中提供的一种致密砂岩气甜点的预测方法的双曲递减趋势图;FIG2 is a hyperbolic decreasing trend diagram of a method for predicting tight sandstone gas sweet spots provided in one embodiment;
图3为一个实施例中提供的一种致密砂岩气甜点的预测方法的预测交会图。FIG. 3 is a prediction cross-plot of a method for predicting tight sandstone gas sweet spots provided in an embodiment.
具体实施方式Detailed ways
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。In order to make the purpose, technical solution and advantages of the present invention more clearly understood, the present invention is further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not intended to limit the present invention.
一个实施例中,提供的一种致密砂岩气甜点的预测方法,如图1所示,该方法包括:In one embodiment, a method for predicting tight sandstone gas sweet spots is provided, as shown in FIG1 , and the method includes:
1、绘制致密砂岩气井的累产气量、月产气量与生产年月的交会图。1. Draw a cross-plot of the cumulative gas production, monthly gas production and production year and month of tight sandstone gas wells.
2、根据交会图判识累产气量的递减趋势,通过递减趋势计算出致密砂岩气的单井累产气量。2. According to the intersection diagram, the decreasing trend of cumulative gas production is identified, and the cumulative gas production of a single well of tight sandstone gas is calculated based on the decreasing trend.
2.1、根据交会图判识致密砂岩气井单井的累产气量的递减趋势。2.1. Identify the decreasing trend of the cumulative gas production of a single tight sandstone gas well based on the intersection diagram.
(1)致密砂岩气井的累产气量与月产气量在交会图上呈直线关系时,单井累产气量的递减趋势为指数递减。(1) When the cumulative gas production and monthly gas production of tight sandstone gas wells present a linear relationship on the intersection diagram, the decreasing trend of the cumulative gas production of a single well is an exponential decrease.
(2)致密砂岩气井的累产气量与月产气量在交会图上呈半对数直线关系时,单井累产气量的递减趋势为调和递减。(2) When the cumulative gas production and monthly gas production of tight sandstone gas wells present a semi-logarithmic linear relationship on the cross-plot, the decreasing trend of the cumulative gas production of a single well is a harmonic decrease.
(3)上述二者特征均不明显,则单井累产气量的递减趋势为双曲递减。(3) If neither of the above two characteristics is obvious, the decreasing trend of the cumulative gas production of a single well is a hyperbolic decrease.
图2左上角为单井大吉6-7的交会图,图2右上角为单井大吉7-8的交会图,图2左下角为单井大吉5-4向4的交会图,图2右下角为单井1-7向6的交会图。The upper left corner of Figure 2 is the intersection diagram of single well Daji 6-7, the upper right corner of Figure 2 is the intersection diagram of single well Daji 7-8, the lower left corner of Figure 2 is the intersection diagram of single well Daji 5-4 to 4, and the lower right corner of Figure 2 is the intersection diagram of single well 1-7 to 6.
2.2、通过生产数据校正的单井产量递减趋势计算出致密砂岩气井单井累产气量值;2.2. Calculate the cumulative gas production value of a tight sandstone gas well by using the decline trend of single well production corrected by production data;
最常用的方法为Arps提出的产量递减规律方程式:The most commonly used method is the law of declining returns proposed by Arps:
Q / Qi = (D / Di )n Q / Qi = (D / Di ) n
式中,Q = 产气量,104 m3/d;Qi = 递减阶段的初始产量,104 m3/d;D = 产量递减率,1/月;Di = 初始瞬时递减率;n = 递减指数;i为正整数。其中,基于已有钻井试气成果,包括但不限于产气量(Q),递减阶段的初始产量(Q i )和试采时间(t)。进而通过产气量(Q)/试采时间(t)计算出产量递减率(D),通过递减阶段的初始产量(Q i )/试采时间(t)计算出初始瞬时递减率(D i )。Where, Q = gas production, 10 4 m 3 /d; Qi = initial production in the decline phase, 10 4 m 3 /d; D = production decline rate, 1/month; Di = initial instantaneous decline rate; n = decline index; i is a positive integer. Based on the existing drilling and gas test results, including but not limited to gas production ( Q ), initial production in the decline phase ( Qi ) and test production time (t). Then, the production decline rate (D) is calculated by gas production ( Q )/test production time (t), and the initial instantaneous decline rate ( Di ) is calculated by initial production in the decline phase ( Qi ) /test production time ( t ).
将递减趋势分为3种类型:指数递减、双曲递减和调和递减;The decreasing trend is divided into three types: exponential decreasing, hyperbolic decreasing and harmonic decreasing;
递减趋势为指数递减时,月产与累产在普通坐标系上呈直线关系,递减指数n →∞;When the decreasing trend is exponential, the monthly output and cumulative output are in a linear relationship on the ordinary coordinate system, and the decreasing exponent n →∞;
递减趋势为调和递减时,月产与累产呈半对数直线关系,递减指数n = 1;When the decreasing trend is harmonic decreasing, the monthly output and cumulative output are in a semi-logarithmic linear relationship, and the decreasing index n = 1;
递减趋势为双曲递减时,1<递减指数n<∞:When the decreasing trend is hyperbolic, 1<decreasing exponent n<∞:
Q = Qi / {1 + (Di / n)×t}n Q = Qi /{1+( Di /n)×t} n
式中,Q = 产气量,104 m3/d;Qi = 递减阶段的初始产量,104 m3/d;Di = 初始瞬时递减率,1/月;t为递减阶段生产时间;n = 递减指数;i为正整数。Where, Q = gas production, 10 4 m 3 /d; Qi = initial production in the decline phase, 10 4 m 3 /d; Di = initial instantaneous decline rate, 1/month; t is the production time in the decline phase; n = decline index; i is a positive integer.
表1 某区块致密砂岩气典型单井第7年累产气量预测与实际生产对比Table 1 Comparison between the predicted and actual cumulative gas production of a typical single well in a certain block of tight sandstone gas in the seventh year
应用双曲线递减规律计算第7年累产气量。如表1所示,通过大宁-吉县区块投产时间较长的典型井用以验证,发现双曲递减预测的累产气量与实际生产数据吻合度高,进一步验证本研究区上古生界致密气递减符合双曲递减规律,且该递减规律符合研究区致密气生产动态趋势。The hyperbolic decline law is applied to calculate the cumulative gas production in the seventh year. As shown in Table 1, through the typical wells with a long production time in the Daning-Jixian block, it is found that the cumulative gas production predicted by the hyperbolic decline is highly consistent with the actual production data, which further verifies that the decline of the tight gas in the Upper Paleozoic in this study area conforms to the hyperbolic decline law, and this decline law conforms to the dynamic trend of tight gas production in the study area.
3、获取致密砂岩气的静态控制因素,所述静态控制因素包括:孔隙度Φ、含气饱和度S g 和单层砂体厚度H,联立所有的静态控制因素创建富集系数EF。3. Obtain the static control factors of tight sandstone gas, which include porosity Φ, gas saturation Sg and single-layer sand body thickness H. Combine all the static control factors to create the enrichment factor EF.
3.1、优选致密砂岩气静态控制因素,如单层砂体厚度(H)、孔隙度(Φ)和含气饱和度(S g )等相关参数。3.1. Optimize the static control factors of tight sandstone gas, such as single-layer sand body thickness (H), porosity (Φ) and gas saturation ( S g ) and other related parameters.
(1)孔隙度Φ(1) Porosity Φ
测井资料计算储层孔隙度主要依赖中子、密度、声波3种孔隙度测井方法。在均质碎屑岩地层这3种孔隙度测井方法反映的地层孔隙度没有本质区别,考虑到致密砂岩在纵、横向上的非均质性较强,可能会影响单一因素反映的孔隙度值,因此选取自然伽马测井参数GR、中子测井参数CNL、声波测井参数AC和密度测井参数DEN的测井参数组合计算孔隙度。基于单井1000m3和3000m3产量的试气测试成果约束,分别读取大宁-吉县区块盒8+山2+本溪钻井气层单层砂体解释组合参数。根据统计分析结果,分别优选出1000m3和3000m3产量所对应的自然伽马测井参数GR、中子测井参数CNL、声波测井参数AC和密度测井参数DEN与孔隙度Φ交会关系,建立自然伽马测井参数GR、中子测井参数CNL、声波测井参数AC、密度测井参数DEN测井组合参数与孔隙度Φ的模型。The calculation of reservoir porosity by logging data mainly relies on three porosity logging methods: neutron, density, and sonic. There is no essential difference in the formation porosity reflected by these three porosity logging methods in homogeneous clastic rock formations. Considering that the tight sandstone has strong vertical and horizontal heterogeneity, which may affect the porosity value reflected by a single factor, the logging parameter combination of natural gamma logging parameter GR, neutron logging parameter CNL, sonic logging parameter AC, and density logging parameter DEN is selected to calculate the porosity. Based on the constraints of the gas test results of 1000m3 and 3000m3 production of a single well, the interpretation combination parameters of the single-layer sand body of the gas layer of He8+Shan2+Benxi drilling in Daning-Jixian block are read respectively. According to the statistical analysis results, the intersection relationships between the natural gamma logging parameters GR, neutron logging parameters CNL, acoustic logging parameters AC and density logging parameters DEN and the porosity Φ corresponding to the production of 1000m3 and 3000m3 were respectively selected, and a model of the combined logging parameters of the natural gamma logging parameters GR, neutron logging parameters CNL, acoustic logging parameters AC, density logging parameters DEN and the porosity Φ was established.
根据模型全面开展研究区单层砂体孔隙度的解释:测井资料计算储层孔隙度主要依赖中子、密度、声波3种孔隙度测井方法。在均质碎屑岩地层这3种孔隙度测井方法反映的地层孔隙度没有本质区别,考虑到致密砂岩在纵、横向上的非均质性较强,可能会影响单一因素反映的孔隙度值,因此选取基于自然伽马测井参数GR和中子测井参数CNL评价砂岩,之后,通过声波测井参数AC和密度测井参数DEN测井参数组合计算孔隙度。取二者计算值的平均值为该套砂岩的计算孔隙度值。The porosity of the single-layer sand body in the study area was fully interpreted according to the model: the calculation of reservoir porosity by logging data mainly relies on three porosity logging methods: neutron, density, and acoustic wave. There is no essential difference in the formation porosity reflected by these three porosity logging methods in homogeneous clastic rock formations. Considering that the tight sandstone has strong vertical and horizontal heterogeneity, it may affect the porosity value reflected by a single factor. Therefore, the sandstone is evaluated based on the natural gamma logging parameter GR and the neutron logging parameter CNL. After that, the porosity is calculated by the combination of the acoustic logging parameter AC and the density logging parameter DEN. The average value of the two calculated values is taken as the calculated porosity value of the sandstone.
利用Wyllie的声波时差方程计算泥岩孔隙度(Wyllie et al., 1956):The porosity of mudstone is calculated using Wyllie's acoustic time difference equation (Wyllie et al., 1956):
Φ AC = (Δt-Δt ma ) / (Δt f -Δt ma ) (1) Φ AC = ( Δt-Δt ma ) / ( Δt f - Δt ma ) (1)
式中,Δt是声波传播时间,埋深z处砂岩测井曲线数值,Δt f 是孔隙流体声波传播时间,Δt ma 是基质(或颗粒)声波传播时间,Φ AC代表埋深z处砂岩孔隙度。Where, Δt is the acoustic wave propagation time, Δt f is the acoustic wave propagation time of the pore fluid, Δt ma is the acoustic wave propagation time of the matrix (or particles), and Φ AC represents the porosity of the sandstone at the burial depth z.
密度测井计算泥岩孔隙度由以下公式确定(Dasgupta et al., 2016):The mudstone porosity calculated by density logging is determined by the following formula (Dasgupta et al., 2016):
Φ DEN = (ρ ma -ρ b ) / (ρ ma -ρ f ) (2) Φ DEN = ( ρ ma - ρ b ) / ( ρ ma - ρ f ) (2)
式中,ρ ma 为砂岩基质密度,ρ b 为密度测井中获得的砂岩体积密度,ρ f 为孔隙流体密度。Where ρma is the sandstone matrix density, ρb is the sandstone volume density obtained from density logging, and ρf is the pore fluid density.
(2)含气饱和度S g (2) Gas saturation S g
针对低渗透砂岩含气储层,通过岩石电阻率实验,分析了岩石Archie参数特征,在此基础上利用Archie公式建立了含水饱和度模型。利用高压压汞实验和核磁共振实验,结合研究区致密气实际成藏地质背景共同校验了含水饱和度的可靠性。For low permeability sandstone gas-bearing reservoirs, the rock Archie parameter characteristics were analyzed through rock resistivity experiments, and the water saturation model was established using the Archie formula. The reliability of water saturation was verified by high-pressure mercury injection experiments and nuclear magnetic resonance experiments combined with the actual geological background of tight gas accumulation in the study area.
岩石电阻率实验参考《SY T 5385-2007 岩石电阻率参数实验室测量及计算方法》标准。利用Archie公式,分别依次计算地层因素(F),电阻增大率(I),胶结指数(m)和饱和度指数(n)。对来自同一层组的多块岩样,在双对数坐标下,绘制地层因素(F)与孔隙度(Φ)的交会图、电阻增大率(I)与含水饱和度(S w )的交会图。用幂函数拟合地层因素(F)与孔隙度(Φ)曲线、阻增大率(I)与含水饱和度(S w )曲线,得到岩性系数(a和b)及胶结指数(m)、饱和度指数(n)。Archie公式为:The rock resistivity experiment refers to the standard "SY T 5385-2007 Laboratory Measurement and Calculation Method of Rock Resistivity Parameters". Using the Archie formula, the formation factor ( F ), resistivity increase rate ( I ), cementation index ( m ) and saturation index ( n ) are calculated in turn. For multiple rock samples from the same formation, the cross plot of formation factor ( F ) and porosity ( Φ ) and the cross plot of resistivity increase rate ( I ) and water saturation ( S w ) are plotted in double logarithmic coordinates. The formation factor ( F ) and porosity ( Φ ) curves, and the resistivity increase rate ( I ) and water saturation ( S w ) curves are fitted with a power function to obtain the lithology coefficients ( a and b ) and cementation index ( m ), saturation index ( n ). The Archie formula is:
其中,S w 为含水饱和度,a、b为岩性系数,R W 为地层水电阻率,R t 为地层电阻率,Ф为孔隙度,m为胶结系数,n为饱和度指数。Among them, Sw is water saturation, a and b are lithology coefficients , RW is formation water resistivity, Rt is formation resistivity, Ф is porosity , m is cementation coefficient, and n is saturation index.
(3)单层砂体厚度H(3) Single-layer sand body thickness H
根据单层砂体的孔隙度Φ和含气饱和度S g 确定单层砂体厚度,具体为:基于单井1000m3和3000m3产量约束下的气层解释,优选出大宁-吉县区块盒8+山2+本溪气层单层砂体解释组合参数取值为自然伽马测井参数GR、中子测井参数CNL、声波测井参数AC、密度测井参数DEN、孔隙度Φ、含气饱和度S g ,根据测井组合参数开展单层砂体厚度的解释。基于已有钻井的试气测试成果,也就是根据单井1000m3和3000m3产量的实测约束,分别读取大宁-吉县区块盒8+山2+本溪钻井气层单层砂体解释组合参数。根据统计分析结果,分别优选出1000m3和3000m3产量所对应的自然伽马测井参数GR、中子测井参数CNL、声波测井参数AC、密度测井参数DEN、孔隙度Φ和含气饱和度S g 测井参数的下限值,建立自然伽马测井参数GR、中子测井参数CNL、声波测井参数AC、密度测井参数DEN、孔隙度Φ、含气饱和度S g 的测井组合参数模型,开展研究区单层砂体气层厚度的解释,如实施例1。The thickness of the single-layer sand body is determined according to the porosity Φ and gas saturation Sg of the single-layer sand body. Specifically, based on the interpretation of the gas layer under the constraints of single-well production of 1000m3 and 3000m3 , the interpretation combination parameters of the single-layer sand body in the gas layer of He8+Shan2+Benxi in the Daning-Jixian block are selected as the natural gamma logging parameter GR, neutron logging parameter CNL, acoustic logging parameter AC, density logging parameter DEN, porosity Φ, and gas saturation Sg . The single-layer sand body thickness is interpreted according to the logging combination parameters. Based on the gas test results of existing wells, that is, according to the measured constraints of single-well production of 1000m3 and 3000m3 , the interpretation combination parameters of the single-layer sand body in the gas layer of He8+Shan2+Benxi in the Daning-Jixian block are read respectively. According to the statistical analysis results, the lower limits of the natural gamma logging parameter GR, neutron logging parameter CNL, acoustic logging parameter AC, density logging parameter DEN, porosity Φ and gas saturation Sg logging parameters corresponding to the production of 1000m3 and 3000m3 are respectively selected, and a logging combination parameter model of the natural gamma logging parameter GR, neutron logging parameter CNL, acoustic logging parameter AC, density logging parameter DEN, porosity Φ and gas saturation Sg is established to interpret the thickness of the single-layer sand body gas layer in the study area, as shown in Example 1.
实施例1Example 1
大宁-吉县区块盒8+山2+本溪气层单层砂体解释组合参数选取Selection of combined interpretation parameters for single-layer sand bodies in He8+Shan2+Benxi gas layers in Daning-Jixian block
单井日产量为1000m3时:声波测井参数AC≥195μs/m,电阻率值M2RX≥70Ω.m,中子测井参数CNL≤18.0%,密度测井参数DEN≤2.67g/cm3,自然伽马测井参数GR≤92API。When the daily production of a single well is 1000m 3 : the acoustic logging parameter AC≥195μs/m, the resistivity value M2RX≥70Ω.m, the neutron logging parameter CNL≤18.0%, the density logging parameter DEN≤2.67g/cm 3 , and the natural gamma logging parameter GR≤92API.
单井日产量为3000m3时:声波测井参数AC≥195μs/m,电阻率值M2RX≥70Ω.m,中子测井参数CNL≤16.0%,密度测井参数DEN≤2.64g/cm3,自然伽马测井参数GR≤90API。When the daily production of a single well is 3000m 3 : the acoustic logging parameter AC≥195μs/m, the resistivity value M2RX≥70Ω.m, the neutron logging parameter CNL≤16.0%, the density logging parameter DEN≤2.64g/cm 3 , and the natural gamma logging parameter GR≤90API.
通过90天单井日产气量对大宁-吉县区块上古生界二叠系山西组山2段电性值约束如下:The electrical property values of the Shanxi Formation Shan 2 Section of the Upper Paleozoic Permian in the Daning-Jixian block are constrained by the daily gas production of a single well for 90 days as follows:
单井日产量为1000m3时:含气饱和度S g ≥66%,测井计算孔隙度Φ≥4.1%,声波测井参数AC≥195μs/m,中子测井参数CNL≤19.5%。When the daily production of a single well is 1000m3 : gas saturation Sg ≥66%, porosity calculated by logging Φ ≥4.1 %, acoustic logging parameter AC≥195μs/m, neutron logging parameter CNL≤19.5%.
单井日产量为3000m3时:含气饱和度S g ≥75%,测井计算孔隙度Φ≥5.5%,声波测井参数AC≥195μs/m,中子测井参数CNL≤19.5%。When the daily production of a single well is 3000m3 : gas saturation Sg ≥75%, porosity calculated by logging Φ ≥5.5 %, acoustic logging parameter AC≥195μs/m, neutron logging parameter CNL≤19.5%.
3.2、联立上述静态参数,创建富集系数(EF),即EF = H×Φ×S g ,用于评价致密砂岩气甜点。3.2. Combining the above static parameters, the enrichment factor (EF) is established, that is, EF = H × Φ × S g , which is used to evaluate the tight sandstone gas sweet spot.
4、分别绘制单井累产气量与孔隙度Φ、含气饱和度S g 、富集系数EF的预测交会图以及单层砂体厚度H与含气饱和度S g 的交会图。4. Draw the predicted cross-plots of the cumulative gas production of a single well and porosity Φ , gas saturation S g , and enrichment factor EF , as well as the cross-plots of the thickness H of a single sand body and gas saturation S g .
根据预测交会图得到致密砂岩气富集的下限值,所述下限值即为致密砂岩气甜点值,如图3所示:According to the predicted cross-plot, the lower limit of tight sandstone gas enrichment is obtained, which is the sweet spot of tight sandstone gas, as shown in Figure 3:
4.1、通过单井累产气量与单层砂体厚度(H)交会,建立致密砂岩气藏富集度单层砂体厚度评价标准(如图3(a)所示);4.1. By intersecting the cumulative gas production of a single well with the thickness of a single sand body (H), an evaluation standard for the thickness of a single sand body in terms of the richness of a tight sandstone gas reservoir is established (as shown in Figure 3 (a)).
4.2、通过单层砂体厚度(H)与含气饱和度(S g )交会,建立致密砂岩气富集度含气饱和度评价标准(如图3(b)所示);4.2. By intersecting the single-layer sand body thickness (H) and gas saturation ( S g ), the evaluation standard of tight sandstone gas richness and gas saturation is established (as shown in Figure 3 (b));
4.3、通过单井累产气量与孔隙度(Φ)交会,建立致密砂岩气藏富集度孔隙度评价标准(如图3(c)所示);4.3. By intersecting the cumulative gas production of a single well with the porosity (Φ), the richness and porosity evaluation standard of tight sandstone gas reservoirs is established (as shown in Figure 3 (c));
4.4、通过单井累产气量与富集系数(EF)交会,建立致密砂岩气藏富集度富集系数评价标准(如图3(d)所示)。4.4. By intersecting the cumulative gas production of a single well with the enrichment factor (EF), an evaluation standard for the enrichment factor of tight sandstone gas reservoirs is established (as shown in Figure 3 (d)).
4.5、根据所有的预测交会图得到致密砂岩气富集度的下限值,所述下限值为致密砂岩气甜点值,如表2所示。表中,I类类别为7年累产气0.077×108m3,相当于7年内单井日产3000m3,II类类别为7年累产气0.026×108m3,相当于7年内单井日产1000m3,以上为中石油某区块工业规模气层下限。4.5. The lower limit of the tight sandstone gas enrichment is obtained according to all the prediction cross-plots. The lower limit is the sweet spot of tight sandstone gas, as shown in Table 2. In the table, the cumulative gas production of Class I is 0.077×10 8 m 3 in 7 years, which is equivalent to 3000m 3 per day per well in 7 years, and the cumulative gas production of Class II is 0.026×10 8 m 3 in 7 years, which is equivalent to 1000m 3 per day per well in 7 years. The above is the lower limit of industrial-scale gas layers in a certain block of PetroChina.
表2 鄂尔多斯盆地某区块致密砂岩气重点层位甜点评价标准Table 2 Evaluation criteria for sweet spots in key strata of tight sandstone gas in a block of the Ordos Basin
一个实施例中,提供的一种致密砂岩气甜点的预测系统,包括:In one embodiment, a prediction system for tight sandstone gas sweet spots is provided, comprising:
绘制模块,用于绘制致密砂岩气井的累产气量、月产气量与生产年月的交会图。The drawing module is used to draw the intersection diagram of the cumulative gas production, monthly gas production and production year and month of tight sandstone gas wells.
判识模块,用于根据交会图判识单井累产气量的递减趋势,通过递减趋势计算出致密砂岩气的单井累产气量。The identification module is used to identify the decreasing trend of the cumulative gas production of a single well according to the intersection diagram, and calculate the cumulative gas production of a single well of tight sandstone gas through the decreasing trend.
联立模块,用于获取致密砂岩气的静态控制因素,联立所有的静态控制因素创建富集系数EF。The combined module is used to obtain the static control factors of tight sandstone gas and to create the enrichment factor EF by combining all the static control factors.
二次绘制模块,用于分别绘制单井累产气量与孔隙度Φ、单层砂体厚度H、富集系数EF的预测交会图以及单层砂体厚度H与含气饱和度S g 的预测交会图。The secondary drawing module is used to draw the predicted cross-plot of the cumulative gas production of a single well and the porosity Φ , the single-layer sand body thickness H , and the enrichment factor EF , as well as the predicted cross-plot of the single-layer sand body thickness H and the gas saturation Sg .
推断模块,用于根据所有的预测交会图得到致密砂岩气富集度的下限值,所述下限值为致密砂岩气甜点值。The inference module is used to obtain the lower limit value of the tight sandstone gas enrichment according to all the predicted intersection diagrams, and the lower limit value is the tight sandstone gas sweet spot value.
以上所述实施例仅表达了本发明的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附权利要求为准。The above-mentioned embodiments only express several implementation methods of the present invention, and the descriptions thereof are relatively specific and detailed, but they cannot be understood as limiting the scope of the invention patent. It should be pointed out that, for ordinary technicians in this field, several variations and improvements can be made without departing from the concept of the present invention, and these all belong to the protection scope of the present invention. Therefore, the protection scope of the patent of the present invention shall be subject to the attached claims.
Claims (5)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410206324.4A CN117784278B (en) | 2024-02-26 | 2024-02-26 | Prediction method and prediction system for dense sandstone gas dessert |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410206324.4A CN117784278B (en) | 2024-02-26 | 2024-02-26 | Prediction method and prediction system for dense sandstone gas dessert |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117784278A CN117784278A (en) | 2024-03-29 |
CN117784278B true CN117784278B (en) | 2024-05-17 |
Family
ID=90402062
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202410206324.4A Active CN117784278B (en) | 2024-02-26 | 2024-02-26 | Prediction method and prediction system for dense sandstone gas dessert |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117784278B (en) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU2020103953A4 (en) * | 2020-12-08 | 2021-02-11 | Southwest Petroleum University | Method and system for predicting production of fractured horizontal well in shale gas reservoir |
CN112526107A (en) * | 2020-11-27 | 2021-03-19 | 中国地质大学(北京) | Method for recognizing and quantitatively characterizing desserts in fractured compact sandstone reservoir |
CN116402225A (en) * | 2023-04-13 | 2023-07-07 | 西南石油大学 | A Gas Production Prediction Method for Tight Sandstone Gas Reservoirs |
-
2024
- 2024-02-26 CN CN202410206324.4A patent/CN117784278B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112526107A (en) * | 2020-11-27 | 2021-03-19 | 中国地质大学(北京) | Method for recognizing and quantitatively characterizing desserts in fractured compact sandstone reservoir |
AU2020103953A4 (en) * | 2020-12-08 | 2021-02-11 | Southwest Petroleum University | Method and system for predicting production of fractured horizontal well in shale gas reservoir |
CN116402225A (en) * | 2023-04-13 | 2023-07-07 | 西南石油大学 | A Gas Production Prediction Method for Tight Sandstone Gas Reservoirs |
Non-Patent Citations (3)
Title |
---|
准噶尔盆地致密砾岩储层特征及"甜点"预测——以中佳地区佳木河组二段一砂组储层为例;杨川等;天然气地球科学;20230228(第2期);全文 * |
运城盆地油气地质与成藏条件研究;朱瑞静;工程科技Ⅰ辑;20230515(第5期);全文 * |
鄂尔多斯盆地七里村油田长7段致密储层孔隙结构特征;赵靖舟;第十七届全国古地理学及沉积学学术会议摘要集;20230831;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN117784278A (en) | 2024-03-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106468172B (en) | A logging interpretation method for low-resistivity reservoirs in ultra-low permeability sandstone reservoirs | |
CN110318745B (en) | Particle size lithology logging evaluation method under deposition microphase constraint | |
CN104564041B (en) | Hyposmosis clastic reservoir rock efficiency evaluation method based on exploitation permeability limits | |
CN111271057B (en) | A method for identifying gas layers by well logging | |
CN110688781B (en) | Well logging interpretation method for low-permeability heterogeneous gas reservoir | |
CN110344822B (en) | Determination method of water saturation in calcareous tight reservoir | |
CN107829731B (en) | Clay alteration volcanic porosity correction method | |
CN112965114B (en) | Dessert evaluation method for offshore deep natural gas reservoir | |
CN113503156B (en) | Quantitative determination method, system, medium, equipment and terminal for vertical heterogeneity of reservoir | |
CN106777649A (en) | A kind of fractured reservoir pore structure method for quantitatively evaluating | |
CN114428365A (en) | Gas-containing characteristic-based tight sandstone gas reservoir logging identification method and device | |
CN115542419A (en) | Comprehensive post-drilling evaluation method, system and medium for marine gas hydrate | |
CN117784278B (en) | Prediction method and prediction system for dense sandstone gas dessert | |
CN112381259B (en) | Dynamic prediction calculation method for productivity of tight conglomerate reservoir | |
CN118964792A (en) | A rapid method for dividing effective reservoirs in deep tight sandstone gas reservoirs | |
CN112346147A (en) | Reservoir evaluation method based on neutron density porosity difference | |
CN118131360A (en) | A method for predicting formation three pressures applicable to carbonate rocks | |
CN108805158A (en) | A kind of fine and close oily reservoir diagenetic phase division methods | |
CN112081582A (en) | Prediction method, system and device for dominant channel in water-drive oil reservoir development | |
Jin | An integrated geomechanics and petrophysics study of hydraulic fracturing in naturally fractured reservoirs | |
CN115434698B (en) | Stratum calcareous cement content calculation method based on photoelectric absorption section index | |
CN115012903B (en) | Logging evaluation method for judging shale bedding structure development | |
CN116050656A (en) | Prediction method of shale brittleness index based on principal component analysis and neural network | |
CN109917489B (en) | A New Method for Determining the Underground Confined Water Level | |
CN117709108B (en) | Model building method for oil layer porosity of grape-flowered oil field Fuyu |
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 |