CN110685651A - Yield splitting method and system for multilayer commingled production gas well - Google Patents

Yield splitting method and system for multilayer commingled production gas well Download PDF

Info

Publication number
CN110685651A
CN110685651A CN201910972695.2A CN201910972695A CN110685651A CN 110685651 A CN110685651 A CN 110685651A CN 201910972695 A CN201910972695 A CN 201910972695A CN 110685651 A CN110685651 A CN 110685651A
Authority
CN
China
Prior art keywords
gas
annual
contribution rate
yield
well
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
Application number
CN201910972695.2A
Other languages
Chinese (zh)
Other versions
CN110685651B (en
Inventor
庞进
刘洪�
卢灿阳
其他发明人请求不公开姓名
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chongqing University of Science and Technology
Original Assignee
Chongqing University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Chongqing University of Science and Technology filed Critical Chongqing University of Science and Technology
Priority to CN201910972695.2A priority Critical patent/CN110685651B/en
Publication of CN110685651A publication Critical patent/CN110685651A/en
Application granted granted Critical
Publication of CN110685651B publication Critical patent/CN110685651B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH DRILLING; MINING
    • E21BEARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B43/00Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
    • E21B43/14Obtaining from a multiple-zone well
    • EFIXED CONSTRUCTIONS
    • E21EARTH DRILLING; MINING
    • E21BEARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00Survey of boreholes or wells
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/02Agriculture; Fishing; Mining

Abstract

The invention discloses a yield splitting method and system for a multilayer commingled production gas well, and relates to the technical field of oil and gas field development. First, various beneficial factors and adverse factors influencing gas reservoir permeability are fully considered, wherein the beneficial factors comprise formation coefficient and gas saturation, and the adverse factors are mainly shale content. And then searching influence weights of all factors by constructing a comprehensive relation influencing the yield of the gas well, calculating the layered yield contribution rate of each gas-containing small layer by utilizing a multiple linear regression algorithm, and finally splitting by combining the specific test and production condition of the single well. By adopting the method and the system provided by the invention, the splitting result of the gas well yield can be accurately and effectively determined.

Description

Yield splitting method and system for multilayer commingled production gas well
Technical Field
The invention relates to the technical field of oil and gas field development, in particular to a method and a system for splitting the yield of a multilayer commingled gas production well.
Background
Multilayer commingling is an important way for the development of multilayer gas reservoirs, but the amount of each small layer of yield contribution is difficult to accurately determine in the gas well commingling process. The existing yield splitting method mainly comprises two types. The first type directly utilizes the gas production profile to split, and the method is based on the actual test output of each small layer during multi-layer commingled production, and the splitting result is more accurate and reliable, but because the test workload is larger and the cost is high, the gas well can not be normally produced in the test process, and the gas field can not be adopted in a large scale usually. And the second type is splitting by using logging information, the yield of the gas well and the formation coefficient (Kh) are in a linear relation according to a gas well productivity formula, and the formation coefficient obtained by logging interpretation calculation is in a direct proportional relation with the yield of each small layer.
Disclosure of Invention
The invention aims to provide a method and a system for splitting the yield of a multilayer commingled gas well, which can accurately and effectively determine the splitting result of the yield of the gas well.
In order to achieve the purpose, the invention provides the following scheme:
a multi-layer commingled production gas well yield splitting method comprises the following steps:
according to the annual profile test data and the well logging interpretation parameters of the gas reservoir components of the single well, establishing the comprehensive relation between the annual layered yield contribution rate of the gas reservoir components and the well logging interpretation parameters;
according to section test data produced by the gas layer component years of the single well, calculating the historical average yield contribution rate of the gas layer component years in layering;
according to the gas layer component annual stratification historical average yield contribution rate and the comprehensive relationship between the gas layer component annual stratification yield contribution rate and the logging interpretation parameters, fitting and determining the weight coefficient of each logging interpretation parameter in the comprehensive relationship by adopting a multiple linear regression algorithm, and further determining the final comprehensive relationship between the gas layer component annual stratification yield contribution rate and the logging interpretation parameters;
reversely deducing and determining a calculation relational expression of the annual and hierarchical yield contribution rate of the gas formation component according to the final comprehensive relationship between the annual and hierarchical yield contribution rate of the gas formation component and the well logging interpretation parameters;
and calculating a relational expression according to the annual profile test data of the gas layer components of the single well and the annual layered yield contribution rate of the gas layer components determined by reverse deduction, and splitting the yield of the single well.
Optionally, the establishing a comprehensive relationship between the contribution rate of the annual layered production of the gas formation component and the logging interpretation parameter according to the profile test data and the logging interpretation parameter produced by the individual well gas formation component year specifically includes:
screening well logging interpretation parameters; the well logging interpretation parameters comprise the shale content, the formation coefficient reflecting the seepage capability of the reservoir and the gas saturation reflecting the gas-bearing capability of the gas layer;
calculating parameters of each gas-containing small layer according to the well logging interpretation parameters;
and constructing a comprehensive relation between the annual stratified yield contribution rate of the gas layer components and the well logging interpretation parameters according to the gas-containing small layer parameters of each layer.
Optionally, the comprehensive relationship between the annual stratified production contribution rate of the gas formation components and the well logging interpretation parameter is as follows:
Fi=Xi·x+Yi·y+Zi·z
wherein: fiThe comprehensive coefficient of the ith layer; x + y + z is 1, and x, y and z are weight coefficients of each logging interpretation parameter;
Figure BDA0002232610820000021
kh is a formation coefficient, and n is the total number of commingled production layers;
Figure BDA0002232610820000022
Sgis the gas saturation;
Figure BDA0002232610820000023
Vshis the argillaceous content.
Optionally, the determining, according to the gas formation component annual stratification historical average yield contribution rate and the comprehensive relationship between the gas formation component annual stratification yield contribution rate and the logging interpretation parameters, a multiple linear regression algorithm is adopted to fit and determine the weight coefficient of each logging interpretation parameter in the comprehensive relationship, so as to determine the final comprehensive relationship between the gas formation component annual stratification yield contribution rate and the logging interpretation parameters, and specifically includes:
and fitting the weight coefficient of each logging interpretation parameter in the comprehensive relation by adopting a multiple linear regression algorithm according to the annual stratified historical average yield contribution rate of the gas layer component and the comprehensive relation between the annual stratified yield contribution rate of the gas layer component and the logging interpretation parameters, and selecting the weight coefficient set when the fitting result is the maximum comprehensive coefficient, thereby determining the final comprehensive relation between the annual stratified yield contribution rate of the gas layer component and the logging interpretation parameters.
Optionally, the method further includes: and splitting the single well yield according to the annual stratification historical average yield contribution rate of the gas formation component when the final comprehensive relationship between the annual stratification yield contribution rate of the gas formation component and the well logging interpretation parameters cannot be established.
Optionally, the calculating a relational expression according to the annual stratified production contribution rate of the gas formation component determined by the individual well gas formation component and the reverse deduction to split the individual well production specifically includes:
when output profile test data are not measured all the time after the single well is put into production, calculating the annual and stratified yield contribution rate of the gas formation component according to the gas formation component annual and stratified yield contribution rate calculation relational expression, and then splitting the single well yield according to the annual and stratified yield contribution rate of the gas formation component;
when only once production profile test data exists after the single well is put into production, calculating a relational expression according to the annual stratified production contribution rate of the gas formation components and the annual single well production before the production profile test date, and then splitting the single well production according to the annual stratified production contribution rate of the gas formation components;
when the single well is put into production and has multiple output profile test data, splitting the single well output according to three test conditions, wherein the three test conditions are as follows:
in the first test condition, when a single well starts to be tested before a first output profile, calculating a relational expression according to the gas formation component annual stratified output contribution rate and first output profile test data to calculate the gas formation component annual stratified output contribution rate, and then splitting the single well output according to the gas formation component annual stratified output contribution rate;
in a second practical situation, when a single well is among a plurality of times of output profile tests, obtaining process output profile test data by adopting a linear interpolation mode according to the adjacent two times of output profile test data, then calculating a relational expression according to the gas layer component annual hierarchical yield contribution rate and the process output profile test data, and finally splitting the single well yield according to the gas layer component annual hierarchical yield contribution rate;
in a third practical situation, when the single well is tested at the last production profile till now, the gas formation component annual stratified production contribution rate is calculated according to the gas formation component annual stratified production contribution rate calculation relational expression and the last production profile test data, and then the single well production is split according to the gas formation component annual stratified production contribution rate.
In order to achieve the above purpose, the invention also provides the following scheme:
a multi-layer commingled gas well yield split system comprising:
the comprehensive relationship establishing module is used for establishing a comprehensive relationship between the annual and layered yield contribution rate of the gas formation components and the logging interpretation parameters according to the annual profile test data and the logging interpretation parameters of the gas formation components of the single well;
the average yield contribution rate calculation module is used for calculating the historical average yield contribution rate of the gas formation component year layering according to profile test data produced by the single-well gas formation component year;
the comprehensive relationship determination module is used for fitting and determining the weight coefficient of each logging interpretation parameter in the comprehensive relationship by adopting a multiple linear regression algorithm according to the gas layer component annual stratification historical average yield contribution rate and the comprehensive relationship between the gas layer component annual stratification yield contribution rate and the logging interpretation parameters, and further determining the final comprehensive relationship between the gas layer component annual stratification yield contribution rate and the logging interpretation parameters;
the backstepping module is used for backstepping and determining a calculation relational expression of the annual and hierarchical yield contribution rate of the gas formation composition according to the final comprehensive relationship between the annual and hierarchical yield contribution rate of the gas formation composition and the well logging interpretation parameters;
the first splitting module is used for splitting the single well yield according to profile test data produced by the single well gas layer component years and a relational expression calculated by reversely deducing the contribution rate of the gas layer component year layered yield.
Optionally, the comprehensive relationship establishing module specifically includes:
the screening unit is used for screening the well logging interpretation parameters; the well logging interpretation parameters comprise the shale content, the formation coefficient reflecting the seepage capability of the reservoir and the gas saturation reflecting the gas-bearing capability of the gas layer;
the calculation unit is used for calculating parameters of each layer of gas-containing small layer according to the well logging interpretation parameters;
and the comprehensive relation establishing unit is used for establishing a comprehensive relation between the annual stratified yield contribution rate of the gas formation components and the well logging interpretation parameters according to the gas-containing small layer parameters of each layer.
Optionally, the comprehensive relationship determining module specifically includes:
and the comprehensive relationship determining unit is used for fitting the weight coefficient of each logging interpretation parameter in the comprehensive relationship by adopting a multiple linear regression algorithm according to the gas layer component annual stratification historical average yield contribution rate and the comprehensive relationship between the gas layer component annual stratification yield contribution rate and the logging interpretation parameters, selecting the weight coefficient set when the fitting result is the maximum comprehensive coefficient, and further determining the final comprehensive relationship between the gas layer component annual stratification yield contribution rate and the logging interpretation parameters.
Optionally, the multi-layer commingled gas well yield splitting system further includes: and the second splitting module is used for splitting the single-well yield according to the annual layering historical average yield contribution rate of the gas formation component when the comprehensive relation between the final annual layering yield contribution rate of the gas formation component and the well logging interpretation parameter cannot be established.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects:
the invention fully considers various beneficial factors and adverse factors influencing the gas reservoir permeability, wherein the beneficial factors comprise formation coefficient and gas saturation, and the adverse factor is mainly shale content. And (3) searching the influence weight of each factor by constructing a comprehensive coefficient (comprehensive relation) influencing the yield of the gas well, and calculating the layered yield contribution rate of each small layer by using a multivariate linear regression method. And finally splitting according to the specific test and production conditions of the single well. By adopting the method and the system provided by the invention, the splitting result of the gas well yield can be accurately and effectively determined.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings needed to be used in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings without inventive exercise.
FIG. 1 is a schematic flow chart of a multi-layer commingled gas well yield splitting method according to an embodiment of the invention;
FIG. 2 is a flow chart of a split flow for a single well specific test and production scenario in accordance with an embodiment of the present invention;
fig. 3 is a schematic structural diagram of a multi-layer commingled gas well yield splitting system in an embodiment of the invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The invention aims to provide a method and a system for splitting the yield of a multilayer commingled gas well, which can accurately and effectively determine the splitting result of the yield of the gas well.
In order to make the aforementioned objects, features and advantages of the present invention comprehensible, embodiments accompanied with figures are described in further detail below.
As shown in fig. 1, the method for splitting the yield of the multi-layer commingled producing gas well provided by the invention comprises the following steps:
step 101: and establishing a comprehensive relation between the annual layered yield contribution rate of the gas formation components and the logging interpretation parameters according to the annual profile test data and the logging interpretation parameters of the gas formation components of the single well.
The comprehensive relationship establishment process of the annual stratified yield contribution rate of the gas layer components and the well logging interpretation parameters is as follows:
(1) and screening well logging interpretation parameters.
The well logging interpretation parameters include favorable parameters and unfavorable parameters.
The favorable parameters selected are formation coefficient (Kh) reflecting the permeability of the reservoir, and gas saturation (S) reflecting the gas bearing capacity of the gas reservoirg) The unfavorable parameter is lithologic parameter-argillaceous content (V)sh)。
(2) And calculating parameters of each layer of gas-containing small layer according to the well logging interpretation parameters.
For the formation coefficient and the gas saturation, the larger the value of the gas saturation, the better the performance of the gas-bearing stratum and the larger the contribution of the layered production. Dividing the layering parameter value by the sum of the layering parameter values to obtain:
in the formula
Kh-formation coefficient, millidarxim mD.m;
x-parameter, dimensionless, decimal;
subscript i — ith layer;
n-total number of layers;
for the argillaceous content (V)sh) In other words, the smaller the value, the better the gas-bearing stratum performance and the greater the contribution to the stratified output. By dividing the reciprocal of the value of this parameter byThe sum of the reciprocal values of the layer stratification parameter values, namely:
Figure BDA0002232610820000071
in the formula
Vsh-argillaceous content, dimensionless, fractional;
y-parameter, dimensionless, decimal.
(3) And constructing a comprehensive relation between the annual stratified yield contribution rate of the gas layer components and the well logging interpretation parameters according to the gas-containing small layer parameters of each layer.
Setting different weight coefficients x, y and z for the three logging interpretation parameters respectively to obtain a comprehensive coefficient F:
Fi=Xi·x+Yi·y+Zi·z(4);
wherein: fiThe comprehensive coefficient of the ith layer; x + y + z is 1.
Step 102: and calculating the historical average yield contribution rate of the gas formation component year layering according to profile test data produced by the single-well gas formation component year. The calculation formula is as follows:
in the formula
Figure BDA0002232610820000073
-gas layer component annual stratification historical average yield contribution;
Vgthe annual stratification historical yield contribution rate of gas layer components, dimensionless, decimal;
z-parameter, dimensionless, decimal.
Step 103: and fitting and determining the weight coefficient of each logging interpretation parameter in the comprehensive relation by adopting a multiple linear regression algorithm according to the gas formation component annual stratification historical average yield contribution rate and the comprehensive relation between the gas formation component annual stratification yield contribution rate and the logging interpretation parameters, and further determining the final comprehensive relation between the gas formation component annual stratification yield contribution rate and the logging interpretation parameters. The method specifically comprises the following steps:
adopting a multiple linear regression algorithm, carrying out multi-parameter linear fitting on the formula (4) and the formula (5), observing fitting effects of different weight coefficients, then continuously adjusting the weight coefficients, selecting the weight coefficient set when the selected fitting result is the maximum comprehensive coefficient, and obtaining the final comprehensive relation between the annual stratified yield contribution rate of the gas layer components and the well logging interpretation parameters, wherein the formula comprises the following steps:
Figure BDA0002232610820000081
in the formula: a, b-linear fitting coefficient.
Step 104: and determining a calculation relational expression of the annual and hierarchical yield contribution rate of the gas formation component by reverse deduction according to the final comprehensive relationship between the annual and hierarchical yield contribution rate of the gas formation component and the well logging interpretation parameters.
The relation is as follows:
Figure BDA0002232610820000082
in the formula:
Figure BDA0002232610820000083
-gas layer component annual stratified yield contribution rate.
Step 105: and calculating a relational expression according to the annual profile test data of the gas layer components of the single well and the annual layered yield contribution rate of the gas layer components determined by reverse deduction, and splitting the yield of the single well.
Splitting is performed according to the specific test and production conditions of a single well, as shown in fig. 2, specifically including:
(1) and when the output profile test data is not measured all the time after the single well is put into production, calculating the annual and stratified output contribution rate of the gas formation component according to the established annual and stratified output contribution rate of the gas formation component, and splitting the single well output according to the annual and stratified output contribution rate of the gas formation component.
In a special case, if the comprehensive relationship cannot be established, the calculation relationship of the annual and hierarchical yield contribution rates of the gas layer components cannot be established, and at the moment, the single-well yield is split according to the statistical annual and hierarchical historical average yield contribution rate of the gas layer components. If the sum of the calculated or counted annual stratification historical average yield contribution rates of the gas layer components is not equal to 1, normalization processing is firstly carried out, and then splitting is carried out, wherein the detailed description is shown in FIG. 2.
(2) When only one-time production profile test data exists after the single well is put into production, the gas formation component annual stratified production contribution rate is calculated according to the gas formation component annual stratified production contribution rate calculation relational expression and the single well production annually before the production profile test date, and then the single well production is split according to the gas formation component annual stratified production contribution rate.
In special cases, if the comprehensive relationship cannot be established, the calculation relationship of the annual stratified yield contribution rate of the gas layer components cannot be established, splitting is carried out according to the statistical annual stratified historical average yield contribution rate of the gas layer components, and then splitting is carried out by using the test data of the output profile every year.
(3) When the test data of the production profile of a plurality of times exist after single well delivery, splitting is carried out according to the following three test conditions:
in the first test case, when the single well starts to the first time of the production profile test, the gas layer component annual stratified production contribution rate is calculated according to the gas layer component annual stratified production contribution rate calculation relational expression and the first time of the production profile test data, and then the single well production is split according to the gas layer component annual stratified production contribution rate.
In the second practical situation, when a single well is among a plurality of times of output profile tests, process output profile test data is obtained in a linear interpolation mode according to the adjacent two times of output profile test data, then a relational expression is calculated according to the gas layer component annual stratified output contribution rate and the process output profile test data, and finally the single well output is split according to the gas layer component annual stratified output contribution rate.
In a third practical situation, when the single well is tested at the last production profile till now, the gas formation component annual stratified production contribution rate is calculated according to the gas formation component annual stratified production contribution rate calculation relational expression and the last production profile test data, and then the single well production is split according to the gas formation component annual stratified production contribution rate.
In order to achieve the above object, the present invention further provides a multi-layer commingled producing gas well yield splitting system, as shown in fig. 3, including:
the comprehensive relationship establishing module 201 is configured to establish a comprehensive relationship between the annual and stratified production contribution rate of the gas formation component and the logging interpretation parameter according to the profile test data and the logging interpretation parameter produced by the individual well gas formation component.
The average production contribution rate calculating module 202 is configured to calculate an average production contribution rate of a gas formation component year stratification history according to profile test data produced by a single well gas formation component year.
And the comprehensive relationship determining module 203 is configured to fit and determine the weight coefficient of each logging interpretation parameter in the comprehensive relationship by using a multiple linear regression algorithm according to the gas formation component annual stratification historical average yield contribution rate and the comprehensive relationship between the gas formation component annual stratification yield contribution rate and the logging interpretation parameters, and further determine a final comprehensive relationship between the gas formation component annual stratification yield contribution rate and the logging interpretation parameters.
And a reverse-pushing module 204, configured to determine a calculation relation of the annual and hierarchical yield contribution rates of the gas formation components by reverse-pushing according to the final comprehensive relationship between the annual and hierarchical yield contribution rates of the gas formation components and the well logging interpretation parameters.
The first splitting module 205 is configured to split the single-well production according to the profile test data produced by the single-well gas formation component year and the gas formation component year layered production contribution rate calculation relational expression determined by reverse extrapolation.
The comprehensive relationship establishing module 201 specifically includes:
the screening unit is used for screening the well logging interpretation parameters; the well logging interpretation parameters comprise the shale content, the formation coefficient reflecting the seepage capability of the reservoir and the gas saturation reflecting the gas bearing capability of the gas reservoir.
And the calculating unit is used for calculating the parameters of the gas-containing small layer of each layer according to the well logging interpretation parameters.
And the comprehensive relation establishing unit is used for establishing a comprehensive relation between the annual stratified yield contribution rate of the gas formation components and the well logging interpretation parameters according to the gas-containing small layer parameters of each layer.
The comprehensive relationship determining module 203 specifically includes:
and the comprehensive relationship determining unit is used for fitting the weight coefficient of each logging interpretation parameter in the comprehensive relationship by adopting a multiple linear regression algorithm according to the gas layer component annual stratification historical average yield contribution rate and the comprehensive relationship between the gas layer component annual stratification yield contribution rate and the logging interpretation parameters, selecting the weight coefficient set when the fitting result is the maximum comprehensive coefficient, and further determining the final comprehensive relationship between the gas layer component annual stratification yield contribution rate and the logging interpretation parameters.
Preferably, the multi-layer commingled gas well yield splitting system further comprises: and the second splitting module is used for splitting the single-well yield according to the annual layering historical average yield contribution rate of the gas formation component when the comprehensive relation between the final annual layering yield contribution rate of the gas formation component and the well logging interpretation parameter cannot be established.
The invention provides a method and a system for splitting the yield of a multi-layer commingled production gas well, which are used for determining the yield contribution rate of each small layer of the multi-layer commingled production gas reservoir. The accuracy of the yield splitting method is obviously higher than that of the conventional method, and the yield splitting method is particularly more adaptive to gas reservoirs with higher mud content and water saturation; in addition, the well does not need to be shut down, the normal production of the gas field is not influenced, the test can be carried out on the basis of the presence or absence of a gas production profile test, the test cost is saved, and the requirements of complex multi-layer gas reservoir development can be met.
The embodiments in the present description are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments are referred to each other. For the system disclosed by the embodiment, the description is relatively simple because the system corresponds to the method disclosed by the embodiment, and the relevant points can be referred to the method part for description.
The principles and embodiments of the present invention have been described herein using specific examples, which are provided only to help understand the method and the core concept of the present invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, the specific embodiments and the application range may be changed. In view of the above, the present disclosure should not be construed as limiting the invention.

Claims (10)

1. A multi-layer commingled production gas well yield splitting method is characterized by comprising the following steps:
according to the annual profile test data and the well logging interpretation parameters of the gas reservoir components of the single well, establishing the comprehensive relation between the annual layered yield contribution rate of the gas reservoir components and the well logging interpretation parameters;
according to section test data produced by the gas layer component years of the single well, calculating the historical average yield contribution rate of the gas layer component years in layering;
according to the gas layer component annual stratification historical average yield contribution rate and the comprehensive relationship between the gas layer component annual stratification yield contribution rate and the logging interpretation parameters, fitting and determining the weight coefficient of each logging interpretation parameter in the comprehensive relationship by adopting a multiple linear regression algorithm, and further determining the final comprehensive relationship between the gas layer component annual stratification yield contribution rate and the logging interpretation parameters;
reversely deducing and determining a calculation relational expression of the annual and hierarchical yield contribution rate of the gas formation component according to the final comprehensive relationship between the annual and hierarchical yield contribution rate of the gas formation component and the well logging interpretation parameters;
and calculating a relational expression according to the annual profile test data of the gas layer components of the single well and the annual layered yield contribution rate of the gas layer components determined by reverse deduction, and splitting the yield of the single well.
2. The method for splitting the yield of the multi-layer commingled gas well as recited in claim 1, wherein the step of establishing a comprehensive relationship between the annual and stratified yield contribution rate of the gas formation component and the logging interpretation parameter according to the annual production profile test data and the logging interpretation parameter of the gas formation component of the single well specifically comprises the steps of:
screening well logging interpretation parameters; the well logging interpretation parameters comprise the shale content, the formation coefficient reflecting the seepage capability of the reservoir and the gas saturation reflecting the gas-bearing capability of the gas layer;
calculating parameters of each gas-containing small layer according to the well logging interpretation parameters;
and constructing a comprehensive relation between the annual stratified yield contribution rate of the gas layer components and the well logging interpretation parameters according to the gas-containing small layer parameters of each layer.
3. The method for splitting the yield of the multi-layer commingled gas well as recited in claim 1, wherein the comprehensive relationship between the annual and stratified yield contribution rate of the gas layer components and the logging interpretation parameters is as follows:
Fi=Xi·x+Yi·y+Zi·z
wherein: fiThe comprehensive coefficient of the ith layer; x + y + z is 1, and x, y and z are weight coefficients of each logging interpretation parameter;
Figure FDA0002232610810000021
kh is a formation coefficient, and n is the total number of commingled production layers;
Figure FDA0002232610810000022
Sgis the gas saturation;
Figure FDA0002232610810000023
Vshis the argillaceous content.
4. The method for splitting the yield of the multi-layer commingled gas production well according to claim 1, wherein the step of fitting and determining the weight coefficient of each logging interpretation parameter in the comprehensive relationship by adopting a multiple linear regression algorithm according to the annual and hierarchical historical average yield contribution rate of the gas layer components and the comprehensive relationship between the annual and hierarchical yield contribution rate of the gas layer components and the logging interpretation parameters further determines the final comprehensive relationship between the annual and hierarchical yield contribution rate of the gas layer components and the logging interpretation parameters specifically comprises the steps of:
and fitting the weight coefficient of each logging interpretation parameter in the comprehensive relation by adopting a multiple linear regression algorithm according to the annual stratified historical average yield contribution rate of the gas layer component and the comprehensive relation between the annual stratified yield contribution rate of the gas layer component and the logging interpretation parameters, and selecting the weight coefficient set when the fitting result is the maximum comprehensive coefficient, thereby determining the final comprehensive relation between the annual stratified yield contribution rate of the gas layer component and the logging interpretation parameters.
5. The method for split production of a multi-commingled producing gas well of claim 1 further comprising: and splitting the single well yield according to the annual stratification historical average yield contribution rate of the gas formation component when the final comprehensive relationship between the annual stratification yield contribution rate of the gas formation component and the well logging interpretation parameters cannot be established.
6. The method for splitting the yield of the multi-layer commingled gas well according to claim 1, wherein the splitting of the single well yield is performed according to the single well gas layer component annual yield profile test data and the gas layer component annual layered yield contribution rate calculation relational expression determined by reverse extrapolation, and specifically comprises the following steps:
when output profile test data are not measured all the time after the single well is put into production, calculating the annual and stratified yield contribution rate of the gas formation component according to the gas formation component annual and stratified yield contribution rate calculation relational expression, and then splitting the single well yield according to the annual and stratified yield contribution rate of the gas formation component;
when only once production profile test data exists after the single well is put into production, calculating a relational expression according to the annual stratified production contribution rate of the gas formation components and the annual single well production before the production profile test date, and then splitting the single well production according to the annual stratified production contribution rate of the gas formation components;
when the single well is put into production and has multiple output profile test data, splitting the single well output according to three test conditions, wherein the three test conditions are as follows:
in the first test condition, when a single well starts to be tested before a first output profile, calculating a relational expression according to the gas formation component annual stratified output contribution rate and first output profile test data to calculate the gas formation component annual stratified output contribution rate, and then splitting the single well output according to the gas formation component annual stratified output contribution rate;
in a second practical situation, when a single well is among a plurality of times of output profile tests, obtaining process output profile test data by adopting a linear interpolation mode according to the adjacent two times of output profile test data, then calculating a relational expression according to the gas layer component annual hierarchical yield contribution rate and the process output profile test data, and finally splitting the single well yield according to the gas layer component annual hierarchical yield contribution rate;
in a third practical situation, when the single well is tested at the last production profile till now, the gas formation component annual stratified production contribution rate is calculated according to the gas formation component annual stratified production contribution rate calculation relational expression and the last production profile test data, and then the single well production is split according to the gas formation component annual stratified production contribution rate.
7. A multi-layer commingled production gas well yield splitting system is characterized by comprising:
the comprehensive relationship establishing module is used for establishing a comprehensive relationship between the annual and layered yield contribution rate of the gas formation components and the logging interpretation parameters according to the annual profile test data and the logging interpretation parameters of the gas formation components of the single well;
the average yield contribution rate calculation module is used for calculating the historical average yield contribution rate of the gas formation component year layering according to profile test data produced by the single-well gas formation component year;
the comprehensive relationship determination module is used for fitting and determining the weight coefficient of each logging interpretation parameter in the comprehensive relationship by adopting a multiple linear regression algorithm according to the gas layer component annual stratification historical average yield contribution rate and the comprehensive relationship between the gas layer component annual stratification yield contribution rate and the logging interpretation parameters, and further determining the final comprehensive relationship between the gas layer component annual stratification yield contribution rate and the logging interpretation parameters;
the backstepping module is used for backstepping and determining a calculation relational expression of the annual and hierarchical yield contribution rate of the gas formation composition according to the final comprehensive relationship between the annual and hierarchical yield contribution rate of the gas formation composition and the well logging interpretation parameters;
the first splitting module is used for splitting the single well yield according to profile test data produced by the single well gas layer component years and a relational expression calculated by reversely deducing the contribution rate of the gas layer component year layered yield.
8. The system for split production of a multi-layer commingled gas well as recited in claim 7, wherein the integrated relationship establishing module specifically comprises:
the screening unit is used for screening the well logging interpretation parameters; the well logging interpretation parameters comprise the shale content, the formation coefficient reflecting the seepage capability of the reservoir and the gas saturation reflecting the gas-bearing capability of the gas layer;
the calculation unit is used for calculating parameters of each layer of gas-containing small layer according to the well logging interpretation parameters;
and the comprehensive relation establishing unit is used for establishing a comprehensive relation between the annual stratified yield contribution rate of the gas formation components and the well logging interpretation parameters according to the gas-containing small layer parameters of each layer.
9. The system for split production of a multi-layer commingled gas well as recited in claim 7, wherein the integrated relationship determination module specifically comprises:
and the comprehensive relationship determining unit is used for fitting the weight coefficient of each logging interpretation parameter in the comprehensive relationship by adopting a multiple linear regression algorithm according to the gas layer component annual stratification historical average yield contribution rate and the comprehensive relationship between the gas layer component annual stratification yield contribution rate and the logging interpretation parameters, selecting the weight coefficient set when the fitting result is the maximum comprehensive coefficient, and further determining the final comprehensive relationship between the gas layer component annual stratification yield contribution rate and the logging interpretation parameters.
10. The multi-layer commingled gas well yield split system of claim 7, further comprising: and the second splitting module is used for splitting the single-well yield according to the annual layering historical average yield contribution rate of the gas formation component when the comprehensive relation between the final annual layering yield contribution rate of the gas formation component and the well logging interpretation parameter cannot be established.
CN201910972695.2A 2019-10-14 2019-10-14 Yield splitting method and system for multilayer commingled production gas well Active CN110685651B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910972695.2A CN110685651B (en) 2019-10-14 2019-10-14 Yield splitting method and system for multilayer commingled production gas well

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910972695.2A CN110685651B (en) 2019-10-14 2019-10-14 Yield splitting method and system for multilayer commingled production gas well

Publications (2)

Publication Number Publication Date
CN110685651A true CN110685651A (en) 2020-01-14
CN110685651B CN110685651B (en) 2021-11-30

Family

ID=69112410

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910972695.2A Active CN110685651B (en) 2019-10-14 2019-10-14 Yield splitting method and system for multilayer commingled production gas well

Country Status (1)

Country Link
CN (1) CN110685651B (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111382523A (en) * 2020-03-17 2020-07-07 中国石油化工股份有限公司 New method for splitting oil well yield
CN111810119A (en) * 2020-07-21 2020-10-23 重庆科技学院 Method for calculating productivity of gas well of high-pressure carbonate rock with water-gas reservoir
CN112031752A (en) * 2020-09-11 2020-12-04 重庆科技学院 Method for calculating stratified formation pressure of multilayer commingled gas well based on flowing pressure test
CN112949974A (en) * 2021-01-11 2021-06-11 中国石油天然气股份有限公司 Method for evaluating contribution rate of layered yield of composite deposition gas reservoir
CN113006760A (en) * 2021-03-17 2021-06-22 西南石油大学 Compact sandstone gas multi-section fracturing horizontal well yield splitting method
CN113762559A (en) * 2020-06-05 2021-12-07 中国石油化工股份有限公司 Intelligent yield splitting method and system for multilayer commingled production gas well

Citations (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1693893A (en) * 2004-04-30 2005-11-09 施卢默格海外有限公司 Method for determining properties of formation fluids
CN101042048A (en) * 2006-03-24 2007-09-26 中国石油天然气股份有限公司 Complicated fault block fluvial facies reservoir oil water well using situation split system
US20090194333A1 (en) * 2007-10-19 2009-08-06 Macdonald Duncan Ranging methods for developing wellbores in subsurface formations
WO2009158160A2 (en) * 2008-06-25 2009-12-30 Schlumberger Canada Limited Method for estimating formation permeability using time lapse measurements
WO2011008405A2 (en) * 2009-06-30 2011-01-20 Bp Corporation North America Inc. Isotopic identification of production by individual formations in commingled gas wells
US20110297370A1 (en) * 2010-06-07 2011-12-08 Conocophillips Company Hydrocarbon production allocation methods and systems
CN102747991A (en) * 2011-04-22 2012-10-24 长江大学 Method for determining single-layer output of commingled producing well
CN104653156A (en) * 2015-01-22 2015-05-27 中国地质大学(武汉) Fragmental rock oil and gas reservoir development unit partitioning method
CN104933246A (en) * 2015-06-15 2015-09-23 中国海洋石油总公司 Method for quantitatively representing contribution of gas field reserves parameters to reserves change
CN105247546A (en) * 2013-06-10 2016-01-13 埃克森美孚上游研究公司 Determining well parameters for optimization of well performance
US20160061020A1 (en) * 2014-08-22 2016-03-03 Chevron U.S.A. Inc. Flooding analysis tool and method thereof
US20160356125A1 (en) * 2015-06-02 2016-12-08 Baker Hughes Incorporated System and method for real-time monitoring and estimation of well system production performance
CN106499370A (en) * 2016-08-22 2017-03-15 中国石油天然气股份有限公司 Well group synchronization dispensing divides each interval Liquid output computational methods of the separate-zone producer that adopts and device
CN106600443A (en) * 2016-12-09 2017-04-26 北京源博科技有限公司 Water saturation-based dynamic oil well yield splitting method
CN106777651A (en) * 2016-12-09 2017-05-31 北京源博科技有限公司 The oil-water well production split method of balanced flood principle
CN106988740A (en) * 2017-06-12 2017-07-28 重庆科技学院 Method based on early yield data prediction shale gas well recoverable reserves
CN107869348A (en) * 2017-10-27 2018-04-03 西北大学 A kind of method of thick-layer sandstone oil reservoir producing well production split
CN109858075A (en) * 2018-12-13 2019-06-07 朱亚婷 A kind of oil well delamination production split method based on unsteady flow

Patent Citations (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1693893A (en) * 2004-04-30 2005-11-09 施卢默格海外有限公司 Method for determining properties of formation fluids
CN101042048A (en) * 2006-03-24 2007-09-26 中国石油天然气股份有限公司 Complicated fault block fluvial facies reservoir oil water well using situation split system
US20090194333A1 (en) * 2007-10-19 2009-08-06 Macdonald Duncan Ranging methods for developing wellbores in subsurface formations
WO2009158160A2 (en) * 2008-06-25 2009-12-30 Schlumberger Canada Limited Method for estimating formation permeability using time lapse measurements
WO2011008405A2 (en) * 2009-06-30 2011-01-20 Bp Corporation North America Inc. Isotopic identification of production by individual formations in commingled gas wells
US20110297370A1 (en) * 2010-06-07 2011-12-08 Conocophillips Company Hydrocarbon production allocation methods and systems
CN102747991A (en) * 2011-04-22 2012-10-24 长江大学 Method for determining single-layer output of commingled producing well
CN105247546A (en) * 2013-06-10 2016-01-13 埃克森美孚上游研究公司 Determining well parameters for optimization of well performance
US20160061020A1 (en) * 2014-08-22 2016-03-03 Chevron U.S.A. Inc. Flooding analysis tool and method thereof
CN104653156A (en) * 2015-01-22 2015-05-27 中国地质大学(武汉) Fragmental rock oil and gas reservoir development unit partitioning method
US20160356125A1 (en) * 2015-06-02 2016-12-08 Baker Hughes Incorporated System and method for real-time monitoring and estimation of well system production performance
CN104933246A (en) * 2015-06-15 2015-09-23 中国海洋石油总公司 Method for quantitatively representing contribution of gas field reserves parameters to reserves change
CN106499370A (en) * 2016-08-22 2017-03-15 中国石油天然气股份有限公司 Well group synchronization dispensing divides each interval Liquid output computational methods of the separate-zone producer that adopts and device
CN106600443A (en) * 2016-12-09 2017-04-26 北京源博科技有限公司 Water saturation-based dynamic oil well yield splitting method
CN106777651A (en) * 2016-12-09 2017-05-31 北京源博科技有限公司 The oil-water well production split method of balanced flood principle
CN106988740A (en) * 2017-06-12 2017-07-28 重庆科技学院 Method based on early yield data prediction shale gas well recoverable reserves
CN107869348A (en) * 2017-10-27 2018-04-03 西北大学 A kind of method of thick-layer sandstone oil reservoir producing well production split
CN109858075A (en) * 2018-12-13 2019-06-07 朱亚婷 A kind of oil well delamination production split method based on unsteady flow

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
PANG, J .ET AL: "Optimization Deployment for Well Patterns of three-Dimensional Injection-production in the Buried-Hill Reservoirs", 《 INTERNATIONAL CONFERENCE ON FUTURE ENERGY, ENVIRONMENT, AND MATERIALS (FEEM)》 *
刘香山等: "多层合采水驱油藏产量劈分方法适应性研究", 《科技信息》 *
庞进等: "基于流态划分的页岩气井产量预测可靠性分析", 《特种油气藏》 *
曾庆恒等: "气田合采层间干扰分析及参数优化", 《油气田地面工程》 *
林孟雄等: "苏里格气田苏东区块多层产量劈分新方法", 《新疆地质》 *
苏伟明等: "稠油油藏多层合采产量劈分方法研究——以辽河油田D29块为例", 《录井工程》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111382523A (en) * 2020-03-17 2020-07-07 中国石油化工股份有限公司 New method for splitting oil well yield
CN113762559A (en) * 2020-06-05 2021-12-07 中国石油化工股份有限公司 Intelligent yield splitting method and system for multilayer commingled production gas well
CN111810119A (en) * 2020-07-21 2020-10-23 重庆科技学院 Method for calculating productivity of gas well of high-pressure carbonate rock with water-gas reservoir
CN111810119B (en) * 2020-07-21 2022-06-28 重庆科技学院 Method for calculating productivity of gas well of high-pressure carbonate rock having water gas reservoir
CN112031752A (en) * 2020-09-11 2020-12-04 重庆科技学院 Method for calculating stratified formation pressure of multilayer commingled gas well based on flowing pressure test
CN112031752B (en) * 2020-09-11 2022-07-01 重庆科技学院 Method for calculating stratified formation pressure of multilayer commingled gas well based on flowing pressure test
CN112949974A (en) * 2021-01-11 2021-06-11 中国石油天然气股份有限公司 Method for evaluating contribution rate of layered yield of composite deposition gas reservoir
CN112949974B (en) * 2021-01-11 2024-03-05 中国石油天然气股份有限公司 Method for evaluating layered yield contribution rate of composite sedimentary gas reservoir
CN113006760A (en) * 2021-03-17 2021-06-22 西南石油大学 Compact sandstone gas multi-section fracturing horizontal well yield splitting method
CN113006760B (en) * 2021-03-17 2022-03-11 西南石油大学 Compact sandstone gas multi-section fracturing horizontal well yield splitting method

Also Published As

Publication number Publication date
CN110685651B (en) 2021-11-30

Similar Documents

Publication Publication Date Title
CN110685651B (en) Yield splitting method and system for multilayer commingled production gas well
CN105569646B (en) A kind of oil/gas well technologically recoverable reserves prediction technique
CN101560879B (en) Well test analysis control method for low permeability gas reservoir and method thereof
CN111425193B (en) Reservoir compressibility evaluation method based on clustering analysis logging rock physical facies division
CN106150477A (en) A kind of method determining single well controlled reserves
CN106600443A (en) Water saturation-based dynamic oil well yield splitting method
CN107038516B (en) Quantitative evaluation method for water-flooding development effect of medium-permeability complex fault block oil reservoir
CN109882163A (en) A kind of PRODUCTION FORECASTING METHODS for compact oil reservoir pressure break horizontal well
CN109284545A (en) The method for solving of structural realism based on optimal conditions importance sampling technique
CN110878690A (en) Method for determining shale gas well production decreasing curve based on carbon isotope evolution
CN112926893A (en) Horizontal well profile control effect evaluation method based on fuzzy comprehensive evaluation and hierarchical analysis
CN111476448B (en) Multi-layer combined production vertical well dynamic yield splitting method
CN113445988A (en) Method for evaluating productivity of gas well of low-permeability carbonate rock gas reservoir
CN113762614A (en) Shale gas well productivity prediction method based on static oil and gas field data
CN111209646B (en) Gas production amount splitting method and device for compact low-permeability sandstone gas reservoir multi-layer commingled production well
CN116976519A (en) Shale oil reservoir single well recoverable reserve prediction method and system
CN108665090B (en) Urban power grid saturation load prediction method based on principal component analysis and Verhulst model
CN110486009A (en) A kind of automatic reverse method of parameter and system on infinity stratum
CN109356577A (en) Tight gas reservoir reserves measuring method based on gas-bearing formation Drilling ratio
CN114462690A (en) Deep shale gas fracturing effect prediction method based on mutation theory
CN111222227B (en) Mine curve determining method and system
CN114066666A (en) Method for analyzing connectivity among wells through injection-production profile monitoring data
CN113094816A (en) Method for constructing comprehensive working condition vibration spectrum and long-life test spectrum of armored vehicle
CN110765660A (en) Method for quickly determining movable reserves of low-permeability gas under different gas production speeds
Guo et al. Application of improved multiple linear regression method in oilfield output forecasting

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
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20200114

Assignee: Guangxi ronghua Ship Technology Co.,Ltd.

Assignor: Chongqing University of Science & Technology

Contract record no.: X2023980053987

Denomination of invention: A Production Splitting Method and System for Multi layer Combined Production Gas Wells

Granted publication date: 20211130

License type: Common License

Record date: 20231227

EE01 Entry into force of recordation of patent licensing contract
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20200114

Assignee: GUANGZHOU JUMING AUTOMATIC CONTROL TECHNOLOGY CO.,LTD.

Assignor: Chongqing University of Science & Technology

Contract record no.: X2023980054328

Denomination of invention: A Production Splitting Method and System for Multi layer Combined Production Gas Wells

Granted publication date: 20211130

License type: Common License

Record date: 20240103

Application publication date: 20200114

Assignee: GUANGZHOU TIANYI TECHNOLOGY Co.,Ltd.

Assignor: Chongqing University of Science & Technology

Contract record no.: X2023980054325

Denomination of invention: A Production Splitting Method and System for Multi layer Combined Production Gas Wells

Granted publication date: 20211130

License type: Common License

Record date: 20240103

Application publication date: 20200114

Assignee: Guangzhou Xiangyun Technology Co.,Ltd.

Assignor: Chongqing University of Science & Technology

Contract record no.: X2023980054323

Denomination of invention: A Production Splitting Method and System for Multi layer Combined Production Gas Wells

Granted publication date: 20211130

License type: Common License

Record date: 20240103

Application publication date: 20200114

Assignee: Laishi (Guangzhou) Digital Technology Co.,Ltd.

Assignor: Chongqing University of Science & Technology

Contract record no.: X2023980054322

Denomination of invention: A Production Splitting Method and System for Multi layer Combined Production Gas Wells

Granted publication date: 20211130

License type: Common License

Record date: 20240103

EE01 Entry into force of recordation of patent licensing contract
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20200114

Assignee: Yuao Holdings Co.,Ltd.

Assignor: Chongqing University of Science & Technology

Contract record no.: X2024980000642

Denomination of invention: A Production Splitting Method and System for Multi layer Combined Production Gas Wells

Granted publication date: 20211130

License type: Common License

Record date: 20240119

Application publication date: 20200114

Assignee: Zhongke Liben (Chongqing) Technology Development Co.,Ltd.

Assignor: Chongqing University of Science & Technology

Contract record no.: X2024980000644

Denomination of invention: A Production Splitting Method and System for Multi layer Combined Production Gas Wells

Granted publication date: 20211130

License type: Common License

Record date: 20240119

EE01 Entry into force of recordation of patent licensing contract
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20200114

Assignee: FOSHAN TENGPU INDUSTRIAL DESIGN CO.,LTD.

Assignor: Chongqing University of Science & Technology

Contract record no.: X2024980003021

Denomination of invention: A Production Splitting Method and System for Multi layer Combined Production Gas Wells

Granted publication date: 20211130

License type: Common License

Record date: 20240322

Application publication date: 20200114

Assignee: FOSHAN YIQING TECHNOLOGY Co.,Ltd.

Assignor: Chongqing University of Science & Technology

Contract record no.: X2024980003019

Denomination of invention: A Production Splitting Method and System for Multi layer Combined Production Gas Wells

Granted publication date: 20211130

License type: Common License

Record date: 20240322