Research on Quantitative Evaluation Methods and Influencing Factors of Natural Gas Reservoir Development

Current methods of evaluating the effects of gas reservoir exploitation primarily focus on evaluating production statistics of a well or field, and they often have a certain subjectivity. In order to further analyze key indexes that affect a gas reservoir during production, this paper considered the characteristic parameters of different types of gas reservoirs and established an index system of evaluating a gas reservoir “exploitation effect” from the three aspects of gas reservoir characteristics, rock characteristics, and economic benefits. *ere are numerous factors affecting the exploitation effect of a gas reservoir. *e relationships between different factors are complicated, and the evaluation result of each factor can also be vague. A fuzzy optimization theoretical model for gas reservoir evaluation and optimization was established on the basis of a quantitative study on the evaluation parameters in fuzzy mathematics. In order to clearly understand the influence of each factor on the gas reservoir exploitation effect, numerical simulation was used to perform a multifactor numerical simulation study of the quantitative sensitivity to determine the sequence of different influencing factors in terms of influencing the gas reservoir exploitation effect. *is study provides a quantitative scientific basis for the evaluation of gas reservoir exploitation effect.


Introduction
Due to complex geological conditions in China, low-quality reserves account for a higher proportion of the proven natural gas reserves. e exploitation with scale benefits is faced with long-term challenges. Accordingly, the evaluation of the gas field exploitation effect has become critical for proper field evaluation. However, in traditional theories and technologies of gas reservoir engineering, there are relatively few contents in this aspect. is has gradually become the technical bottleneck of gas reservoir exploitation. Compared with developed countries, in China, the quantitative analysis that couples natural gas exploitation technologies with economic situation is still lacking, and it is difficult to meet the demand for exploiting conventional and unconventional complex gas reservoirs under the new normal situation. is restricts further improvement of the natural gas exploitation level in China.
A comprehensive evaluation of the gas reservoir exploitation effect is a seemingly simple problem, but in fact, it is complicated. Because there are numerous factors influencing this effect, the interactions and contradictions among the factors increase the difficulty of analysis. It needs to be studied and determined by effective multi-index quantitative comprehensive evaluation methods. Guo and Evans [1] established a net present value model for the economic benefit evaluation of horizontal well exploitation in naturally fractured reservoirs with dissolved gas flooding. Samuel et al. [2] established a simple net present value model to evaluate the economic benefits of horizontal well exploitation by using the steady-state production equation of oil wells. Liu et al. [3] established the calculation formulas of horizontal well exploitation, vertical well exploitation, and their relative payback periods, which have a certain guiding role in the exploitation and deployment of horizontal wells. Wang et al. [4] used the Upper Shaximiao Formation (J 2 s 2 ) gas reservoir of the Middle Jurassic Series in the Xinchang Gas Field of the western Sichuan Depression as a representative. is study enhanced the geologic understanding of the reservoir by comparing the scattered distribution of the permeable reservoir sand body to the continuous distribution of favorable reservoirs. ree breakthroughs were determined in the production-increasing technology from multilayer fracturing of vertical wells to multilayer and multifracture fracturing of horizontal wells. In order to achieve efficient exploitation of a gas field, Wu et al. [5] performed a series of technical studies on the distribution and exploitation mode of organic reef reservoirs, including the fine characterization of small reefs and quantitative prediction of thin reservoirs, the optimization design of horizontal wells in band-like small reef gas reservoir, the real-time optimization and adjustment of horizontal well trajectory in ultradeep and thin reservoirs, and so on. Combined with the optimization exploitation technology policies that considered the influence of water invasion, Wang et al. [6] solved several key problems and developed a series of supporting exploitation technologies suitable for a gas reservoir. ese were applied to the exploitation practice of a gas field. Dong et al. [7] studied the characteristics of the pore structure of natural gas hydrate reservoir in the Qilian Mountain Permafrost. Ding et al. [8] studied the sand control in a natural gas hydrate reservoir in the South China Sea. Sun et al. [9] proposed a new approach to efficient and safe gas production from unsealed marine hydrate deposits. Wang et al. [10] established a new analysis method of cumulative gas production per unit pressure drop, and its effectiveness was verified by its application in the calculation of reserves of three gas reservoirs. Hu et al. [11] believed that the developed special technology for the efficient exploitation of deep-marine carbonate gas reservoirs provided strong support for Pet-roChina Southwest Oil & Gas Field to build a giant gas field of 10 billion cubic meters creating substantial growth of natural gas output. Wang et al. [12] developed a set of key technical systems suitable for the efficient exploitation of complex tight sandstone gas reservoirs in the Yan'an Gas Field located in the southeastern Ordos Basin. e Sulige Gas Field is located in the northern Ordos Basin, and it is China's largest tight sandstone gas field. After decades of research and integrated innovation, a full set of key technologies centered on reservoir characterization [13], full digital earthquake monitoring [14], big well group stereoscopic exploitation [15], "factory-like" fracturing transformation [16], downhole throttling [17], and drainage gas recovery has been developed for the exploration and development in this field [18]. ese technologies have had field-scale exploitation implementation and provide a good reference for the exploitation of tight sandstone gas reservoirs in China [19][20][21][22][23][24][25][26][27][28][29][30][31][32].
In summary, many scholars have performed the technical and economic analysis and classified evaluation for different types of gas reservoirs and developed several empirical modes in terms of exploitation indexes such as reasonable gas well output, the reasonable well spacing density of gas reservoir, and reasonable gas producing rates. However, due to the complexity of geological conditions, the determination of exploitation indexes is still a relatively difficult technical step in the process of gas reservoir exploitation. Moreover, the evaluation methods of gas reservoir exploitation effect mostly look for influences based upon statistics, which are subjective to a certain extent, and cannot directly guide the optimization of specific indexes. erefore, it is necessary to consider a variety of gas reservoir characteristic parameters and establish a comprehensive evaluation index system for gas reservoir development by combining statistics with the technical-economic coupling analysis method.

Gas Reservoir Evaluation Index System
In order to further analyze the important indexes affecting gas reservoirs, this paper separates the corresponding parameter values from the calculation formulas of geological reserves, dynamic reserves, gas well production, and economic return rates. e evaluation index system for the gas reservoir exploitation effect was established with consideration of the three aspects of gas reservoir characteristics, rock characteristics, and economic evaluation. A total of 12 indexes were used to evaluate the gas reservoir, as shown in Figure 1.

Fuzzy Optimization eory Model for Gas Reservoir Evaluation and Optimization.
In fuzzy decision problems, the decision-making objectives and constraint conditions typically have collinearity, and there are no obvious differences between objects. e theory of fuzzy mathematics can well solve fuzzy decision problems and has been widely used in practice.
X is a set composed of n schemes. Namely, it is a domain of discourse as follows: where x i represents a certain candidate evaluation scheme and i � 1, 2, . . ., n. en, fuzzy optimization is to determine the membership degree u A (x i ) of each element x i to the "optimal solution." e larger the membership degree, the closer the gas reservoir is to an "optimal gas reservoir." In other words, the problem of scheme optimization can be transformed into selecting a gas reservoir with the largest membership degree.
In general, the evaluation indexes can be divided into two categories: one consists of the-larger-the-better indexes and the other consists of the-smaller-the-better indexes. e membership degree of the-larger-the-better indexes can be calculated from the following: And the membership degree of the-smaller-the-better indexes can be calculated from the following: 2 Shock and Vibration where r ij represents the membership degree of the j-th evaluation index of gas reservoir i, x ij represents the value of the j-th evaluation index of gas reservoir i, x jmin represents the minimum value of the j-th evaluation index, and x jmax represents the maximum value of the j-th evaluation index. us, the matrix of the membership degree of each gas reservoir evaluation index R ij can be expressed as follows: e degree of membership of each index of the optimal gas reservoir G → can be obtained by evaluating the membership degree matrixes of the indexes. e maximum value among the membership degrees of each index for all gas reservoirs is the membership degree of the m-th index of the optimal gas reservoir G → . us, the optimal gas reservoir G → can be expressed as follows: G → � r 11 ∨r 12 ∨ · · · ∨r 1j , r 21 ∨r 22 ∨ · · · ∨r 2j , . . . , r i1 ∨r i2 · · · r ij T � g 1 , g 2 , . . . , g m T.
(6) Similarly, the minimum value among the membership degrees of each index for all gas reservoirs is the membership degree of the m-th index of the inferior gas reservoir B → . us, the inferior gas reservoir B → can be expressed as follows: where ∨ is the maximum operator and ∧ is the minimum operator.
Notably, different evaluation indexes have different influences on the evaluation results. e degree of influence of each index on the evaluation results can be expressed by their weight ω j , and then the weight vector of the m-th index W �→ can be set as follows: where ω i represents the weight of the j-th index and j � 1, 2, 3, . . ., m.

Shock and Vibration 3
Gas reservoir i can be expressed by the vector R i → as follows: . . , r nj T . (9) en, the difference between the i-th gas reservoir and the optimal gas reservoir can be described by the off-optimal distance as follows: Similarly, the difference between the i-th gas reservoir and the inferior gas reservoir can be described by the offinferior distance as follows: For gas reservoir x i , its membership degree to the optimal gas reservoir is u A (x i ), which is simplified as u i , and its membership degree to the inferior gas reservoir is u B (x i ), which is simplified as u c i . en, u i and u c i have the following relationship: In fuzzy theory, the membership degree can be viewed as a weight. e weighted off-optimal distance can be defined e objective function of an optimal criterion is defined as that the sum of the square of the weighted off-inferior distance and the square of the weighted off-optimal distance of each gas reservoir takes the minimum value, namely: where u i represents the membership degree of the i-th gas reservoir x i to the optimal fuzzy set, G → represents the membership degree matrix of the optimal gas reservoir, B → represents the membership degree matrix of the inferior gas reservoir, W �→ represents the index weight matrix, and R i → represents the membership degree matrix of the i-th gas reservoir. Take the derivation for the objective function (10) and assume the derivation to be equal to 0, namely: en, us, the fuzzy screening model can be obtained as follows: Equation (17) is the fuzzy optimization theory model for gas reservoir evaluation and optimization.

Gas Reservoir Evaluation and Optimization Model.
ere are numerous factors that affect the evaluation results of gas reservoirs, but it is difficult to evaluate each influencing factor. Based on the quantitative study on the evaluation parameters in fuzzy mathematics, a suitable comprehensive analysis method for gas reservoir evaluation and optimization is proposed in this paper.

Gas Reservoir Index Suitability Degree.
e fuzzy membership degree of the j-th evaluation index of the i-th gas reservoir (P ij ) to the optimal value, x ij ∈ (0, 1), can be calculated by the following: where P oj is the most suitable index value and P wj is the least suitable index value. From (17), it can be inferred that the closer the evaluation index of a gas reservoir is to the optimal ideal value, the smaller the fuzzy membership degree is.
In order to ensure that the optimal degree of membership can take a larger value when the evaluation parameters of the gas reservoir are closer to the optimal ideal values and meanwhile enlarge the differences among different evaluation indexes, an exponential function as follows is introduced: us, the fuzzy suitability degree of each evaluation index can be converted into a ij , a ij ∈ (0, 1). When the value of a ij is larger and closer to 1, this index is better. Hence, the fuzzy suitability matrix A of each gas reservoir can be expressed by using a ij as the evaluation index. Each element of the fuzzy suitability matrix A reflects the suitability degree of each index of a certain candidate gas reservoir and the differences of the suitability degree of the same evaluation index among different gas reservoirs.

Establishment of the Comprehensive Evaluation and
Screening Method. After establishing the suitability degree matrix A(a ij ) for m evaluation indexes of n gas reservoirs, the weight vector ω j1 for evaluating the degree of influence of each index considering the different influence degrees of different indexes on the evaluation results of gas reservoirs is required.
e product of element a ij in the fuzzy membership matrix A ij and the weight transposition matrix of each index ω T 1j , which is denoted as W i1 , can be used to reflect the weighted performance of the degree of suitability of the gas reservoir evaluation parameters, namely: e product of W ij and the weighted matrix ω 1j , which is denoted as M ij , can reflect the differences in the different evaluation parameters of different gas reservoirs, namely: Hence, equations (21) and (22) can be used to obtain the evaluation and optimization sequence of all gas reservoirs. e comprehensive evaluation results of gas reservoir evaluation and optimization can be ranked according to the values of V i ′ as follows:

Weight Allocation Method.
For the gas reservoir comprehensive evaluation model, it is necessary to obtain the influence weight of each evaluation index in the gas reservoir. e weight of each evaluation index can be obtained by the fuzzy analytic hierarchy process. Fuzzy analytic hierarchy process (FAHP) is a method to calculate the weights of evaluation indexes by constructing fuzzy judgment, which can reflect the real situation well. e fuzzy judgment matrix R can be established as shown in the following: R � C a 1 a 2 · · · a n a 1 r 11 r 12 · · · r 1n a 2 r 21 r 22 · · · r 2n ⋮ ⋮ ⋮ ⋮ ⋮ a n r n1 r n2 · · · r nn . (23) e meaning of element r ij in the fuzzy judgment matrix is that when the evaluation index a i and the evaluation index a j are compared relative to C (which can be expressed as the recovery factor), the evaluation index a i, and the evaluation index a j follow the rule that a i is more important (or less important) than a j ″ . e degrees of importance of the two evaluation indexes relative to C are quantitatively described in Table 1.
Similarly, the fuzzy membership degrees of the evaluation index a j and the evaluation index a i can be obtained. r ij is defined as the fuzzy membership degree scale of the evaluation index a i to the evaluation index a j , and r ji is defined as the fuzzy membership degree scale of the evaluation index a j to the evaluation index a i . en, r ij and r ji have the following relationship: Further, the fuzzy judgment matrix expressed by the fuzzy membership degree scale can be obtained as follows: R � r 11 r 12 · · · r 1m r 21 r 22 · · · r 2m ⋮ ⋮ ⋮ ⋮ r m1 r m2 · · · r mm R has the following properties: (1) When i � j, r ij � 0.5 (i, j � 1, 2, 3, . . ., m) When the fuzzy judgment matrix R is a consistent fuzzy matrix, the weight ω i can be calculated by the following: In equation (26), a represents the difference between the emphasis degrees for different indexes, which is related to the number of evaluation indexes and the degree of difference. When a � (m − 1)/2, the weights have the largest difference degree; that is, a is equal to (m − 1)/2 when we attach the most important to the importance of evaluation indexes. is value of a is also the used value in our practical applications.
Shock and Vibration 5

Influencing Factors of Gas Reservoir Exploitation Effect
In order to clearly understand the influences of each influencing factor on the gas reservoir exploitation effect, the numerical simulation model was established by using the numerical simulation technology to preliminarily study the influencing rules.

Skin
Factor. e skin factor affects the conductivity between the well and the grid in the numerical model. e larger the skin factor, the worse the exploitation effect. is is because the skin factor affects the conductivity between the wellbore and the grid in the numerical model. As shown in Figures 2-4, the skin factor indeed has a certain influence on the exploitation effect.

Stress Sensitivity.
In the stress-sensitive formation, with a decrease in the formation pressure, the opened fractures under high pressure gradually close and the permeability also simultaneously decreases. is overall affects the flow of the fluids, especially for a gas reservoir with induced fractures. Figures 5-7 show the simulation results. Stress sensitivity leads to the shortening of stable production periods and the reduction of final accumulated gas yield, which greatly affects the exploitation effect of gas reservoirs. Hence, for stress-sensitive gas reservoirs, horizontal wells should be used as far as possible to reduce the bottom-hole flow pressure; meanwhile, a reasonable gas production rate should be maintained to avoid the formation of larger pressure drawdown distribution.

Water Saturation.
Water saturation mainly affects the calculation of effective permeability of a gas well, thus affecting the calculation of conductivity. For different watercut models, it can be seen from Figures 8-10 that water saturation has a great influence on stable gas production and accumulative gas production. erefore, wells should be deployed on the high part of the structure as far as possible so that the perforation segments can be restricted in the regions with high gas saturation as far as possible, thus reducing the risk of water invasion.

Permeability.
Permeability is an important parameter in Darcy's formula, which greatly affects the calculation of conductivity in a numerical simulation. In actual production, permeability has a great influence on the gas well productivity and the gas reservoir exploitation effect. For different permeability models, it can be seen from Figures 11-13 that permeability has a great influence on the stable gas production and accumulative gas production of a gas reservoir. It is considered that permeability is an important index that has a greater influence on the exploitation effect.

Water Invasion.
Water invasion greatly influences gas reservoir production over time. Water invasion can result in a sharp decrease in the gas phase permeability, leading to a decrease in the gas well production capacity. Conversely, water invasion can result in wellbore effusion, and the gas-liquid two-phase flow can increase wellbore friction, thus affecting the gas production ability. Concurrently, the "water encasing gas" phenomenon may appear due to water invasion, making it difficult to exploit this part of the gas reserves. It has been proven in actual production that water invasion has a great influence on the gas reservoir exploitation effect. Figures 14-16 show the numerical simulation results. Once water invasion occurs, the stable production period of the gas reservoir is sharply shortened or even there completely absent. e accumulative gas yield at the end of the forecast period is also greatly reduced. erefore, water invasion is an important index affecting the gas reservoir exploitation effect.
Based on the above analysis and according to the actual production of gas reservoirs, the above influencing factors are ranked as follows: permeability > water saturation >water invasion > stress sensitivity > skin factor. e ranking of influencing factors obtained by numerical simulation provides a reference basis for the weight allocation of influencing factors.

Conclusions
In order to further analyze key indexes that affect a gas reservoir during production, this paper considered the characteristic parameters of different types of gas reservoirs and established an index system of evaluating a gas reservoir "exploitation effect" from the three aspects of gas reservoir characteristics, rock characteristics, and economic benefits.
ere are numerous factors affecting the exploitation effect of a gas reservoir. e relationships between different factors are complicated, and the evaluation result of each factor can also be vague. A fuzzy optimization theoretical model for gas reservoir evaluation and optimization was established on the basis of a quantitative study on the evaluation parameters in fuzzy mathematics.
In order to clearly understand the influence of each factor on the gas reservoir exploitation effect, numerical simulation was used to perform a multifactor numerical simulation study of the quantitative sensitivity to determine the sequence of different influencing factors in terms of influencing the gas reservoir exploitation effect.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that there are no conflicts of interest. Shock and Vibration 13