A gas channeling discrimination method and system based on production well produced fluid phase state

Abstract

The application provides a gas channeling discrimination method and system based on production well output fluid phase state, comprising the following steps: field gas injection effect well analysis; sampling of gas injection effect well flow material; carrying out phase state analysis theoretical simulation and indoor analysis experiment on the obtained sample; comparing the phase state analysis theoretical simulation result and the indoor analysis experiment result to see whether the accuracy requirement is met, if not, repeating the phase state analysis theoretical simulation and the indoor analysis experiment until the accuracy requirement is met; forming a phase diagram envelope line group according to the phase state analysis theoretical simulation and the indoor analysis experiment result; adding a learning mechanism to fit and learn the actual sample parameter point, and establishing a gas channeling discrimination phase diagram graph for different gas oil ratios; and using the gas channeling discrimination phase diagram graph to discriminate gas channeling. Through the research on the phase state change law of the gas injection effect well output fluid, a phase diagram change graph for the whole gas injection process is formed to discriminate gas channeling, and support is provided for oilfield gas injection production regulation and control.

Description

TECHNICAL FIELD

[0001] The present application belongs to the technical field of oilfield development and phase equilibrium theory, and particularly relates to a gas channeling discrimination method and system based on the phase state of production well output fluid. BACKGROUND

[0002] With the increase of unconventional and low-permeability oilfields in the newly added reserves worldwide, and more and more oilfields entering the "two-high and one-low" stage of "high production degree, high water cut and low oil production rate", the gas injection development technology which has good adaptability to the above-mentioned situation has gradually become a commonly used technology for enhancing oil recovery worldwide. At present, the annual oil production obtained by using the gas injection enhanced oil recovery technology has accounted for one third of the annual oil production obtained by the entire enhanced oil recovery technology.

[0003] One of the key technologies of gas injection enhanced oil recovery is how to prolong the contact time of injected gas medium with underground crude oil, that is, to delay the breakthrough of injected gas medium to the production well. How to discriminate gas channeling and form a corresponding discrimination method is particularly important to provide support for field regulation and control.

[0004] Comprehensive main gas channeling discrimination methods can be mainly divided into two categories, one based on injection-production dynamic parameters and the other based on field dynamic monitoring. However, the current gas channeling discrimination method used in the production field still has the basic problem that it cannot accurately judge before gas channeling occurs, but only discriminates when gas channeling occurs. SUMMARY

[0005] In view of the above problems, the present application provides a gas channeling discrimination method and system based on the phase state of production well output fluid, which forms a phase diagram change chart for the entire gas injection process by studying the phase state change law of the injection effect well output fluid to discriminate gas channeling and provide support for oilfield gas injection production regulation and control.

[0006] The present application adopts the following technical solutions:

[0007] A gas channeling discrimination method based on the phase state of production well output fluid, comprising the following steps:

[0008] On-site injection effect well analysis;

[0009] Sampling of well stream matter of the injection effect well;

[0010] Carrying out phase state analysis theoretical simulation and indoor analysis experiment on the sampled samples;

[0011] Comparing the phase state analysis theoretical simulation result and the indoor analysis experiment result to see whether the accuracy requirement is met, if not, repeating the phase state analysis theoretical simulation and the indoor analysis experiment until the accuracy requirement is met;

[0012] Forming a phase diagram envelope line group according to the phase state analysis theoretical simulation and the indoor analysis experiment result;

[0013] A learning mechanism is added to fit and learn the actual sample parameter points, and to establish a gas channeling discrimination phase diagram map for different gas-oil ratios;

[0014] The gas channeling discrimination phase diagram map is used to discriminate gas channeling.

[0015] According to different phase diagram envelope line morphologies, combined with production gas-oil ratios and field production dynamic data, a basic database is formed, a core learning function is formed through a fitting learning method, and a gas channeling discrimination phase diagram map based on different gas-oil ratio phase diagram envelope line morphologies is established; the discrimination map is used to discriminate gas channeling, and as gas injection development proceeds, the phase diagram envelope line gradually transitions from an oil phase to a gas phase, the risk of gas channeling rises, and gas channeling time prediction is fitted from a production dynamic parameter database.

[0016] Preferably, the effective well includes a well with a slowed rising or falling trend of gas production, a well with a slowed rising or falling trend of gas pressure, a well with a slowed rising or falling trend of water content of gas, or a gas well.

[0017] Preferably, the field gas injection effective well analysis parameter includes a bottom hole flowing pressure, an oil / casing pressure, an oil / gas component, and a gas-oil ratio.

[0018] Preferably, the sampling position is a bottom hole or a downhole separator, and the sampling frequency is 2-4 times per month.

[0019] Preferably, after sampling, preparation work for phase state analysis theoretical simulation and indoor analysis experiments is carried out, including the following steps:

[0020] The components of the sample and the molar content of each component are determined;

[0021] Each component is analyzed, and each component is checked and missing components are supplemented by using a reservoir fluid phase state parameter database software.

[0022] Preferably, the objects of the phase state analysis theoretical simulation and indoor analysis experiments include single degassing, constant mass expansion, differential separation, viscosity, separator, bubble and dew point, and constant volume depletion.

[0023] Preferably, the phase state analysis theoretical simulation step includes:

[0024] The sample is divided into components according to physical properties, polarity, and mass transfer similarity;

[0025] State equation selection:

[0026]

[0027] Wherein: p is the system pressure; T is the system temperature; v is the volume of the sampling system; R is the gas constant; ω is the molecular eccentric factor; a is the attraction coefficient; b is the volume coefficient; m, n are the related coefficients reflecting the influence of non-spherical asymmetry and the degree of molecular eccentricity, respectively;

[0028] Phase equilibrium calculation based on thermodynamic equilibrium fugacity:

[0029]

[0030] Wherein: x k , y k , z k are the mole fractions of the kth component in the liquid phase, the gas phase and the total body, respectively; V is the gas phase molar volume, mol; K k is the gas-liquid equilibrium constant of the kth component,

[0031] Selection of initial iteration values:

[0032]

[0033] Wherein: p ci is the critical pressure of the ith component; T ci is the critical temperature of the ith component; ω i is the eccentric factor of the ith component; K i is the gas-liquid equilibrium constant of the ith component.

[0034] Selection of iteration mode:

[0035] Wherein: R k is the iteration residual of the kth component; are the gas phase, liquid phase fugacity factors of the kth component, respectively.

[0036] Preferably, the accuracy requirement is specifically that the general density, volume, bubble point and dew point are all less than 5%, the volume coefficient is less than 10%, and the viscosity and miscibility pressure are less than 20%.

[0037] Preferably, the phase envelope and the isofluid line are obtained by carrying out phase state analysis theoretical simulation and indoor analysis experiment on the sample.

[0038] A gas channeling discrimination system based on the phase state of the production fluid of a production well, comprising

[0039] An analysis module for analyzing a gas injection responsive well in the field;

[0040] A sampling module for sampling the well stream of the gas injection responsive well;

[0041] An experiment module for carrying out phase state analysis theoretical simulation and indoor analysis experiment on the sample obtained by sampling;

[0042] a comparison module for comparing the phase state analysis theory simulation result and the indoor analysis experiment result to see if the accuracy requirement is met, if not, repeating the phase state analysis theory simulation and the indoor analysis experiment until the accuracy requirement is met;

[0043] a drawing module for forming a phase diagram envelope graph set according to the phase state analysis theory simulation and the indoor analysis experiment result;

[0044] a mapping module for adding a learning mechanism to perform fitting and learning on the actual sample parameter point and establishing a gas channeling discrimination phase diagram graph;

[0045] a discrimination module for performing gas channeling discrimination by using the gas channeling discrimination phase diagram graph.

[0046] Preferably, the sampling position of the sampling module is the well bottom or the downhole separator, and the sampling frequency is 2-4 times per month.

[0047] Preferably, the sampling module is used to perform the preparation work of the phase state analysis theory simulation and the indoor analysis experiment after sampling the well stream of the gas injection responsive well, including the following steps:

[0048] determining the components of the sample and the molar content of each component;

[0049] analyzing each component and checking and supplementing the missing components by using the oil and gas reservoir fluid phase state parameter database software.

[0050] Preferably, the object of the experiment module for performing the phase state analysis theory simulation and the indoor analysis experiment on the sample obtained by sampling includes single degassing, constant mass expansion, differential separation, viscosity, separator, bubble and dew point and constant volume depletion.

[0051] Preferably, the phase state analysis theory simulation step in the experiment module for performing the phase state analysis theory simulation and the indoor analysis experiment on the sample obtained by sampling includes:

[0052] performing pseudo-component division on the sample according to the physical property, polarity and mass transfer action similarity;

[0053] state equation selection:

[0054]

[0055] In the formula, p is the system pressure, T is the system temperature, v is the volume of the sampling system, R is the gas constant, ω is the molecular eccentric factor, a is the attraction coefficient, b is the volume coefficient, and m and n are respectively the related coefficients reflecting the influence of nonspherical asymmetry and molecular eccentricity degree;

[0056] phase equilibrium calculation based on thermodynamic equilibrium fugacity:

[0057]

[0058] wherein: x k y k z k are the mole fractions of the liquid phase, vapor phase and total, respectively, of the kth component; V is the vapor molar volume, mol; K k is the vapor-liquid equilibrium constant of the kth component,

[0059] The initial iteration value is selected as:

[0060]

[0061] wherein: p ci is the critical pressure of the ith component; T ci is the critical temperature of the ith component; ω i is the acentric factor of the ith component; K i is the vapor-liquid equilibrium constant of the ith component.

[0062] The iteration mode is selected as:

[0063] wherein: R k is the iteration residual of the kth component; are the fugacity factors of the kth component in the vapor and liquid phases, respectively.

[0064] Preferably, the precision requirement is specifically that the general density, volume, bubble point and dew point are all less than 5%, the volume coefficient is less than 10%, and the viscosity and miscibility pressure are less than 20%.

[0065] Preferably, the phase state analysis theoretical simulation and the indoor analysis experiment are carried out on the sample to obtain the phase envelope and the isofluid line.

[0066] The oil and gas reservoir fluid phase state parameter database software adopted in the present application has been disclosed in the software copyright, and the registration number is 2016SR339677.

[0067] The present application has the beneficial effects that a new gas channeling discrimination method based on the phase state change analysis of the production fluid of the gas injection responsive well is established, a gas channeling discrimination chart based on the phase diagram envelope is formed, the present application can realize the advanced early warning of the gas channeling, has high accuracy, and is simple and convenient to operate, and has excellent economic applicability.

[0068] Other features and advantages of the present application will be set forth in the following description, and in part will become apparent to those skilled in the art from the description, or can be learned by practice of the present application. The objects and other advantages of the present application will be realized and attained by the structure particularly pointed out in the written description and claims hereof as well as the appended drawings. BRIEF DESCRIPTION OF DRAWINGS

[0069] In order to more clearly illustrate the technical solutions in the embodiments of the present application or the prior art, the following will briefly introduce the drawings needed to be used in the embodiments or prior art description. Obviously, the drawings described below are only some embodiments of the present application, and the ordinary skilled in the art can obtain other drawings according to these drawings without any creative effort.

[0070] Figure 1 The gas channeling discrimination process schematic diagram of the present application is shown;

[0071] Figure 2 The gas channeling discrimination chart based on phase diagram envelope under different gas-oil ratios is shown;

[0072] Figure 3 The gas channeling discrimination chart based on phase diagram envelope of actual reservoir DHα is shown;

[0073] Figure 4 The gas channeling discrimination chart based on phase diagram envelope of actual reservoir DHβ is shown. DETAILED DESCRIPTION

[0074] In order to make the objects, technical solutions and advantages of the embodiments of the present application clearer, the following will clearly and completely describe the technical solutions in the embodiments of the present application with reference to the drawings in the embodiments of the present application. Obviously, the described embodiments are only some embodiments of the present application, rather than all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by the ordinary skilled in the art without any creative effort are within the protection scope of the present application.

[0075] Embodiment:

[0076] A gas channeling discrimination method based on production well output fluid phase state, taking actual reservoir (DH) hydrocarbon gas injection gravity drive development as an example, referring to Figure 1 and Figure 2 , comprises the following steps:

[0077] (1) On-site injection effect well analysis, selecting effect wells DHα and DHβ, both of which are wells whose gas effect has not been good in the early stage and has been good at present;

[0078] (2) Sampling of injection effect well stream material, sampling position is downhole separator, frequency is 2 times / month;

[0079] (3) The injection effect well stream material sampling in step (2) is downhole separator sampling, which is used for phase state theoretical simulation and indoor analysis experiment after laboratory recombination verification;

[0080] (4) Preparing work for developing phase state theoretical simulation and indoor analysis experiment, including:

[0081] 1) Laboratory test well stream components and determine the molar content of each component;

[0082] 2) Application component analysis and independent research and development software "Oil and gas reservoir fluid phase state parameter database software" to check and supplement each component and missing components in well stream;

[0083] (5) Develop phase analysis theory simulation and laboratory analysis experiment; Phase analysis theory simulation includes single degassing, CCE (isochoric expansion), CME (constant mass expansion), DL (decline), multiple degassing, viscosity, separator and bubble dew point, etc.

[0084] (6) The phase analysis theory simulation step in step (5) includes:

[0085] 1) The original crude oil sample is divided into six pseudo-components C1+N2, C2+C3+CO2, C4+C5, C6-C9, C 10 -C 29 , C 30+ , see Table 1 according to physical properties, polarity and mass transfer similarity;

[0086] 2) The formula is selected into PR state equation;

[0087]

[0088] 3) Based on the thermodynamic equilibrium fugacity phase equilibrium calculation, the full component model is shown in the formula:

[0089]

[0090] 4) Initial iteration value selection, select Wilson equation, see formula:

[0091]

[0092] 5) Newton iteration is selected for iteration method, and the iteration residual is shown in the formula:

[0093]

[0094] (7) According to the results of laboratory analysis experiment in step (5), the phase envelope and isoflux line are obtained;

[0095] (8) Compare the results obtained in steps (6) and (7) to see if the composite precision requirement is met, if not, repeat steps (5)-(8) until the precision requirement is met; Among them, the general density, volume, bubble point and dew point are less than 5%, the volume coefficient is less than 10%, the viscosity and miscible pressure are less than 20%;

[0096] (9) Forming phase envelope line diagram set, as shown in Figure 1 and Figure 2 ;

[0097] (10) Adding learning mechanism to fit and learn actual sample parameter points; to form a chart, but only the phase diagram obtained from existing indoor experiments, the number may not be enough, so through the learning mechanism to fit the existing experimental data points, predict some parts that have not been experimented, form a complete chart; after the later experiment is supplemented, replace it; according to different phase envelope line forms, combined with production gas-oil ratio and field production dynamic data, form a basic database, through fitting and learning method, form the core learning function, establish the gas channeling discrimination phase diagram based on different gas-oil ratio phase envelope line forms (see Figure 2 , Figure 3 );

[0098] (11) After the formation of the chart, use the gas channeling discrimination phase diagram to discriminate gas channeling, the steps are as follows: after the formation of the chart, the steps are as follows: the chart includes phase envelope line and production gas-oil ratio; different phase envelope lines represent different production periods; with the passage of production time, the temperature and pressure range of the phase envelope line wrapping area presents the characteristics of first increasing and then decreasing; with the passage of production time, the critical point gradually moves from the oil phase to the gas phase, the critical pressure first increases and then decreases, and the critical temperature gradually decreases; when one of the conditions that the critical point migration speed decreases, the critical pressure increases close to the inflection point, and the temperature and pressure range of the phase envelope line wrapping area approaches stability appears and the production gas-oil ratio continues to rise, it can be judged that gas channeling will occur; the temperature and pressure range, critical value and gas-oil ratio in the phase envelope line chart corresponding to different reservoir temperature and pressure conditions and crude oil properties are different, but the law is consistent.

[0099] Exemplarily, take DHα and DHβ wells as examples, based on the above method, form the gas channeling discrimination phase diagram based on the phase envelope line, according to the gas channeling discrimination phase diagram, it can be clearly judged that the two wells occurred gas channeling in 2018, see Figure 3 for DHα well, see Figure 4 for DHβ well.

[0100] Table 1 Sampling different component division method (before gas injection)

[0101]

[0102] A gas channeling discrimination system based on the phase state of production well output fluid, comprising

[0103] an analysis module for analyzing field gas injection effective wells;

[0104] a sampling module for sampling the well stream of the gas injection affected well; the sampling position is the well bottom or the downhole separator, and the sampling frequency is 2-4 times per month; after the sampling of the well stream of the gas injection affected well, the preparation work for the phase state analysis theoretical simulation and the indoor analysis experiment includes the following steps: determining the components of the sample and the molar content of each component; analyzing each component and checking each component by using the oil and gas reservoir fluid phase state parameter database software and supplementing the missing components;

[0105] an experiment module for carrying out the phase state analysis theoretical simulation and the indoor analysis experiment on the sample obtained by sampling; the objects of the experiment include: single degassing, constant mass expansion, differential separation, viscosity, separator, bubble dew point and constant volume depletion;

[0106] a comparison module for comparing the phase state analysis theoretical simulation result and the indoor analysis experiment result to see whether the accuracy requirement is met, if not, repeating the phase state analysis theoretical simulation and the indoor analysis experiment until the accuracy requirement is met; the accuracy requirement is specifically: the general density, volume, bubble point and dew point are all less than 5%, the volume coefficient is less than 10%, the viscosity and the miscible pressure are less than 20%;

[0107] a drawing module for forming a phase diagram envelope line graph set according to the phase state analysis theoretical simulation result and the indoor analysis experiment result;

[0108] a mapping module for adding a learning mechanism to fit and learn the actual sample parameter points and establish a gas channeling discrimination phase diagram graph;

[0109] a discrimination module for discriminating gas channeling by using the gas channeling discrimination phase diagram graph.

[0110] Further, the experiment module is used for carrying out the phase state analysis theoretical simulation and the indoor analysis experiment on the sample obtained by sampling, and the phase state analysis theoretical simulation step in the indoor analysis experiment includes:

[0111] grouping the sample according to the physical properties, polarity and mass transfer action similarity;

[0112] state equation selection:

[0113]

[0114] In the formula: p is the system pressure; T is the system temperature; v is the volume of the sampling system; R is the gas constant; ω is the molecular eccentric factor; a is the attraction coefficient; b is the volume coefficient; m and n are respectively the related coefficients reflecting the influence of non-spherical asymmetry and the degree of molecular eccentricity;

[0115] phase equilibrium calculation based on thermodynamic equilibrium fugacity:

[0116]

[0117] wherein: x k y k z k are the liquid, vapor and overall mole fractions of the kth component, respectively; V is the vapor molar volume, mol; K k is the vapor-liquid equilibrium constant of the kth component,

[0118] Initial iteration value selection:

[0119]

[0120] wherein: p ci is the critical pressure of the ith component; T ci is the critical temperature of the ith component; ω i is the acentric factor of the ith component; K i is the vapor-liquid equilibrium constant of the ith component,

[0121] Iteration mode selection:

[0122] wherein: R k is the iteration residual of the kth component; are the vapor and liquid fugacity factors of the k components, respectively.

[0123] Although the present application has been described in detail with reference to the foregoing embodiments, it should be understood by those skilled in the art that the technical solutions recorded in the foregoing embodiments can still be modified, or some of the technical features can be replaced by equivalents; and such modifications or replacements do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present application.

Claims

1. A gas channeling discrimination method based on the phase behavior of the production well effluent, characterized in that, It comprises the following steps: On-site injection effect well analysis, the effect well comprising: no gas production rising or falling trend slowing well, no gas pressure rising or falling trend slowing well, no gas water cut down or up trend slowing well or gas well; Sampling the injection effect well fluid, after sampling, carrying out phase analysis theoretical simulation and indoor analysis experiment preparation, determining the composition and molar content of each component of the sample, analyzing each component and checking and supplementing the missing components through oil and gas reservoir fluid phase state parameter database software; Carrying out phase analysis theoretical simulation and indoor analysis experiment on the sampled sample; Comparing the phase analysis theoretical simulation result and the indoor analysis experiment result to see whether it meets the accuracy requirement, If not, repeating the phase analysis theoretical simulation and indoor analysis experiment until the accuracy requirement is met; Forming a phase diagram envelope line group according to the phase analysis theoretical simulation and indoor analysis experiment result; Adding a learning mechanism to fit and learn the actual sample parameter points and establish a gas channeling discrimination phase diagram graph for different gas oil ratios; Using the gas channeling discrimination phase diagram graph to discriminate gas channeling; The adding of the learning mechanism to fit and learn the actual sample parameter points and establish a gas channeling discrimination phase diagram graph for different gas oil ratios comprises: According to different phase diagram envelope line morphologies, combining production gas oil ratio and field production dynamic data to form a basic database, through fitting and learning method, forming a core learning function, and establishing a gas channeling discrimination phase diagram graph based on different gas oil ratio phase diagram envelope line morphologies; The using of the gas channeling discrimination phase diagram graph to discriminate gas channeling comprises: The gas channeling discrimination phase diagram graph comprises two parts of phase envelope line and production gas oil ratio, and different phase envelope lines represent different production periods; With the passage of production time, the temperature and pressure range of the phase envelope line wrapped region presents the characteristics of first increasing and then decreasing; With the passage of production time, the critical point gradually moves from the oil phase to the gas phase, the critical pressure first increases and then decreases, and the critical temperature gradually decreases; When the critical point moving speed decreases or the critical pressure increasing speed approaches the inflection point or the temperature and pressure range of the phase envelope line wrapped region approaches stability, and the production gas oil ratio continuously rises, it is judged that gas channeling will occur soon.

2. The gas breakthrough discrimination method based on the phase state of the production well fluid output according to claim 1, characterized in that, The on-site injection effect well analysis parameters comprise bottom hole flowing pressure, tubing / casing pressure, oil / gas components and gas oil ratio.

3. The gas breakthrough discrimination method based on the phase state of the production well fluid output according to claim 1, characterized in that, The sampling position is the bottom hole or downhole separator, and the sampling frequency is 2-4 times / month.

4. The gas coning discrimination method based on the phase behavior of the production well fluid in claim 1, wherein, The objects of the phase analysis theoretical simulation and indoor analysis experiment comprise single degassing, constant mass expansion, differential separation, viscosity, separator, bubble dew point and constant volume depletion.

5. The gas coning discrimination method based on the phase behavior of the production well fluid in claim 1, wherein, The phase analysis theoretical simulation step comprises: Dividing the sample into components according to physical properties, polarity and mass transfer similarity; State equation selection: ， In the formula: p is the system pressure; T is the system temperature; v is the volume of the sampling system; R is the gas constant; ω is the molecular eccentric factor; a is the attraction coefficient; b is the volume coefficient; m and n are respectively the related coefficients reflecting the influence of non-spherical asymmetry and molecular eccentricity; Phase equilibrium calculation based on thermodynamic equilibrium fugacity: ， where: x k ,y k ,z k are the liquid, vapor and overall mole fractions, respectively, of the kth component; V is the vapor molar volume, mol; K k is the vapor-liquid equilibrium constant for the kth component, Initial iteration value selection: ， where: p ci is the critical pressure of the ith component; T ci is the critical temperature of the ith component; ω i is the acentric factor of the ith component; K i is the gas-liquid equilibrium constant of the ith component, Iterative approach: , where: R k is the kth component iteration residual; φ k V ,φ k L are the k component gas and liquid fugacity factors, respectively.

6. The gas coning discrimination method based on the phase behavior of the production well fluid of claim 1, wherein, The accuracy requirement is specifically: less than 5% for general density, volume, bubble point and dew point, less than 10% for volume coefficient, less than 20% for viscosity and miscibility pressure.

7. The gas breakthrough discrimination method based on the phase state of the production well output fluid according to any one of claims 1-6, characterized in that, Theoretical simulation and laboratory analysis experiment of phase state analysis are carried out on the sample.

8. A gas channeling discrimination system based on the phase behavior of the production well effluent, characterized by, The method comprises the following steps: The analysis module is used for analyzing the gas injection effective well, and the effective well comprises: a well with no rising or falling trend of gas production, a well with no rising or falling trend of gas pressure, a well with no rising or falling trend of water content, or a gas well; The sampling module is used for sampling the well stream of the gas injection effective well, and the preparation work of the theoretical simulation and laboratory analysis experiment of phase state analysis after sampling is carried out, the components and the molar content of each component of the sample are determined, each component is analyzed and verified by the oil and gas reservoir fluid phase state parameter database software, and the missing components are supplemented; The experimental module is used for carrying out the theoretical simulation and laboratory analysis experiment of phase state analysis on the sample obtained by sampling; The comparison module is used for comparing the results of the theoretical simulation and laboratory analysis experiment of phase state analysis to see whether the accuracy requirement is met, and if not, repeating the theoretical simulation and laboratory analysis experiment of phase state analysis until the accuracy requirement is met; The drawing module is used for forming a phase diagram envelope line graph set according to the results of the theoretical simulation and laboratory analysis experiment of phase state analysis; The mapping module is used for adding a learning mechanism to fit and learn the actual sample parameter points, and establishing a gas channeling discrimination phase diagram graph, which comprises the following steps: forming a basic database according to different phase diagram envelope line morphologies, combining production gas oil ratio and field production dynamic data, forming a core learning function through fitting and learning methods, and establishing a gas channeling discrimination phase diagram graph based on different gas oil ratio phase diagram envelope line morphologies; The discrimination module is used for discriminating gas channeling by using the gas channeling discrimination phase diagram graph; the gas channeling discrimination phase diagram graph comprises two parts of phase envelope line and production gas oil ratio, and different phase envelope lines represent different production periods; with the lapse of production time, the temperature and pressure range of the phase envelope line wrapped region presents the characteristics of first increasing and then decreasing; with the lapse of production time, the critical point gradually moves from the oil phase to the gas phase, the critical pressure first increases and then decreases, and the critical temperature gradually decreases; when the critical point moving speed decreases or the critical pressure increasing speed approaches the inflection point or the temperature and pressure range of the phase envelope line wrapped region approaches stability, and the production gas oil ratio continuously rises, it is judged that gas channeling will occur.

9. The gas coning discrimination system based on production well fluid phase behavior of claim 8, wherein, The sampling position of the sampling module is the bottom of the well or the downhole separator, and the sampling frequency is 2-4 times per month.

10. The gas coning discrimination system based on production well fluid phase behavior of claim 8, wherein, The sampling module is used for carrying out the preparation work of the theoretical simulation and laboratory analysis experiment of phase state analysis after sampling the well stream of the gas injection effective well, and comprises the following steps: Determine the components and the molar content of each component of the sample; Analyze each component and verify each component by the oil and gas reservoir fluid phase state parameter database software, and supplement the missing components.

11. The gas coning discrimination system based on production well fluid phase behavior of claim 8, wherein, The object of the experimental module for carrying out the theoretical simulation and laboratory analysis experiment of phase state analysis on the sample obtained by sampling comprises: single degassing, constant mass expansion, differential separation, viscosity, separator, bubble point and dew point, and constant volume depletion.

12. The gas coning discrimination system based on production well fluid phase behavior of claim 8, wherein, The step of the theoretical simulation of phase state analysis in the experimental module for carrying out the theoretical simulation and laboratory analysis experiment of phase state analysis on the sample obtained by sampling comprises: The sample is classified according to the physical properties, polarity and mass transfer similarity; State equation selection: ， where: p P is the system pressure; T T is the system temperature; v V is the volume of the sampled system; R R is the gas constant; ω k is the molecular eccentricity factor; a G is the gravitational factor; b V is the volume factor; m , n are the correlation coefficients reflecting the influence of the non-spherical asymmetry and the degree of molecular eccentricity, respectively; Phase equilibrium calculation based on thermodynamic equilibrium fugacity: ， wherein: x k , y k , z k are the liquid, vapor and overall mole fractions of the k th component, respectively; V is the vapor mole volume, mol; K k is the vapor-liquid equilibrium constant of the k th component, Initial iteration value selection: ， wherein: p ci the critical pressure of the i th component; T ci the critical temperature of the i th component; ω i the eccentric factor of the i th component; Ki i the gas-liquid equilibrium constant of theth component, Iterative approach: , In the formula: R k For the first k The iterative residuals of each component; φ k V , φ k L They are respectively k The gas-phase and liquid-phase fugacity factors of each component.

13. The gas coning discrimination system based on production well fluid phase behavior of claim 8, wherein, The accuracy requirement is that the general density, volume, bubble point and dew point are all less than 5%, the volume coefficient is less than 10%, and the viscosity and miscibility pressure are less than 20%.

14. The gas coning discrimination system based on the phase behavior of the production well effluent of any of claims 8-13, wherein, Phase state analysis theoretical simulation and indoor analysis experiment are carried out on the sample to obtain phase envelope and isoflux line.

Images