Facies-Controlled Modeling for Permeability of Tight Gas ReservoirBasedonHydrodynamicandGeophysicsCharacteristics

Due to poor physical properties and strong heterogeneity of Daniudi tight gasfield (China), traditional methods are not ideal for predicting reservoir permeability. Based on geoscience data mining algorithms and modeling techniques, this parameter is predicted and characterized from a new perspective. (e high precision hydraulic unit information was outputted by the BP neural network and the interwell database was predicted by sequential indicator simulation based on the logging and core data, the exponential and power law relationship functions were selected to participate in the construction of the permeability model. (e results showed that the facies-controlled modeling system for permeability based on the hydrodynamic characteristics of tight gas reservoirs could effectively improve the accuracy of reservoir prediction, and the logging information in the longitudinal direction and the facies information in the plane were combined by the hydraulic unit.


Introduction
In recent years, in case of the gradual increase in the scale of exploration and development of tight oil and gas reserves, the reservoir characterization methods have become more abundant [1][2][3]. Permeability is an important parameter for the evaluation of this type of reservoir, and the acquisition of its accurate value directly affects the results of the evaluation of heterogeneity and the design of scheme [4][5][6]. If the data points are unified for regression calculation of the permeability-porosity model, the weakly producing pores and invalid reservoirs could easily be classified into the same category [7,8]. e disadvantages of this method are particularly obvious in the evaluation process of strongly heterogeneous reservoirs [9][10][11]. In recent years, various data mining and deep learning algorithms have attracted attention from reservoir engineers due to their excellent calculation ability [12,13], and this type of algorithms have great potential for improving accuracy [14]. e neural network technology is widely used in the data mining field for the prediction of hydraulic unit (HU)/permeability on the single-well scale, especially BP network. BP network has strong nonlinear mapping capabilities and good effect in solving problems of complex reservoir classification.
Tight gas reservoir has the strong heterogeneity in the distribution of throats-pores environment, and the pores' geometrical characteristic could not only constraint physical properties of reservoir geological body but also directly affect the recovery of tight gasfield [15,16]. e target area has poor microstructure, complex logging response, and poor correlation between overall porosity and permeability. erefore, the traditional methods were abandoned, and the geoscience data mining algorithms were applied to solve the prediction problem of permeability, and the permeability distribution was characterized through the optimized porepermeability function based on the platform of modeling software [17][18][19][20][21][22][23]. e result could effectively eliminate the interference of weak signal-to-noise ratio and improve the prediction accuracy of the model. Generally, the HU/permeability prediction of well trajectory and interwell space were usually two independent tasks; in this study, the prediction work on the vertical and plane scale was innovatively combined; the BP neural network and faciescontrolled simulation system could directly provide an important reference for the prediction and evaluation of related parameters of the same type of reservoirs.

Geological Setting.
e target area is located in the northeast of the Ordos Basin, China, and its construction location is in the northern section of the Yishan Slope ( Figure 1). e slope range is about 6-9 m/km, and the dip angle ranges from 0.3°to 0.6°. is is a monoclinic structure high in the northeast and low in the southwest with sporadic development of nose-like uplifts. e geological characteristics of the H 3 gas-bearing reservoir were studied based on the sedimentary and stratigraphic data of the Xiaohaotu-Da16 well group, Daniudi gasfield. e H 3 reservoir was divided into two layers (H 1 3 and H 2 3 ) according to the rock electricity and cycle differences. e sedimentary environment belongs to river sedimentation, and the lithology of the reservoir was composed of gravel sandstone, medium-fine sandstone intercalated with siltstone. e degree of particle sorting was medium-poor.

Hydrodynamic Parameters and Permeability
Characterization. According to the geological settings, various reservoir characterization parameters have different adaptability. e target area has the characteristics of low permeability, compactness, large seepage changes, and complex pore structure. In the process of constructing the traditional porosity-permeability logarithm linear model, the trend line has segmented characteristics, which will affect the accuracy of the overall permeability prediction model inevitably ( Figure 2). And then, the flow zone indicator (I FZ ) and reservoir quality index (I RQ ) were selected to participate in the prediction and characterization research [24]. Amafule et al. [25] conceptualize the pore throat into a series of capillaries and applied the hydrodynamic method to divide the geometric pore units on the basis of permeability changes: where I FZ reflects the characteristics of rock pores and pore throat and mineral structure, μm and I RQ characterizes the seepage capacity of the reservoir, μm [26].

BP Neural Network.
e BP neural network is composed of an input layer, a hidden layer, and an output layer. e layers are connected by weights. e network training and learning process includes error back propagation and forward propagation of information.
e external information is collected and passed to the hidden layer (each neuron node) through the input layer and further passed to the output layer under the action of the excitation function to complete a forward propagation. e system will enter the back propagation process if the output is inconsistent with expectations. e general principle is to reduce the error between the output and the expectation, the hidden layer and the input layer are conveyed error information through a mode in which the error gradient is gradually reduced, and the weights are also adjusted at the same time. e following formula is the error function expression: where (d k − O k ) corresponds to the difference between the expected and output results. In the back propagation process, the adjustment calculation expression of the weight w is as follows: where the proportional coefficient is represented by a constant, which is the learning rate of the network, and the negative sign represents the gradual decrease of the gradient.

e Hydraulic Unit
Identification of the Core Samples. I FZ of core samples was calculated based on the parameters of porosity and permeability, and then the cluster analysis of 2 Geofluids the sum of squared deviations method was applied to perform calculations since the sum of squares of deviations of the same class is small and the sum of squares of deviations between classes is large, and this method could show a good effect on the classification of hydraulic units [27][28][29][30]. Finally, the hydraulic units in the target area were divided into 4 types. From Figure 3, each type of hydraulic unit corresponds to a fitted trend line with a slope of 1, and the intercept of the trend line is the average I FZ value of this type of hydraulic unit. e average I FZ of HU#1∼HU#4 were 1.3497, 1.1889, 0.9967, and 0.8367, respectively. Because the reservoir is tight and the range of porosity-permeability parameters is small, the difference between the average values of I FZ is also relatively small, but it does not affect the area division of different hydraulic units on the chart of standard porosity-reservoir quality index.

Neural Network Pattern Recognition Technology for
Single-Well. Loggings are the comprehensive parameters that can reflect the geophysical characteristics of the reservoir, and the parameters were introduced into the training process of the artificial neural network. e entire recognition network was composed of input, output, and hidden layer neuron systems under the artificial neural network technology. Logging parameters (RT, AC, SP, GR, and DEN) of the target area were used as input elements, and four types of hydraulic unit types were used as output elements to perform BP network pattern recognition and prediction. e schematic diagram is shown in Figure 4.
According to the empirical formula and the test deduction results, the number of hidden layer neurons was finally determined to be 15, and the sum of the test and validation data weights was set to be 30%. Among them, several key parameters were set as follows: the network learning rate is 0.1; the allowable error is 0.001, and the maximum number of training is 1000 (see Table 1). e network training was stopped execution when the training effect reaches the set requirements.
From the results, the positive judgment rate of test set on the confusion matrix was as high as 95.50%, and the positive judgment rate of the overall set was 98.20% ( Figure 5). e training network was saved, and the sample points were randomly selected for the test on the effect of network prediction. Test results are shown in Figure 6. It can be seen that the statistical frequency of HU prediction on the main diagonal occupies a large proportion, and the overall prediction effect was close to 85% with the environment of insufficient data points. Among them, the prediction accuracy of HU#1, HU#2, HU#3, and HU#4 were 50%, 75%, 81.82%, and 90.48%, respectively. And then, the uncored section of the single well and the uncored well were predicted by the network which provides sufficient hydraulic unit information of single well for the next spatial modeling of the hydraulic unit.
en, the spatial model was further constructed by the sequential indicator simulation method based on the predicted results which were treated as 4 types of discrete variables [31]. e search path values were estimated in each layer and unit, and the layers were simulated under the fitting environment of the variogram. e spatial geological body of H 3 was established finally.

Permeability Model
After determining the spatial distribution model of hydraulic units in the target area, the regression calculation of porosity-permeability parameters was performed according to each type of hydraulic unit. us, each regression model represents the permeability prediction model of each type of hydraulic unit, and then the permeability of each model grid is calculated by the prediction model through the software

Establishment of the Permeability Prediction Model.
Analyzing the past prediction results of physical property parameters, it can be found that it is easy to produce pessimistic results by processing the data to predict the permeability according to the traditional linear method. e main reason is that the part of high-permeability and low-permeability showed obvious heterogeneity changes, which ultimately affects the regression model [32,33]. erefore, this study abandoned the traditional data processing methods and matched the classification results of the hydraulic unit to the permeability model; the porosity-permeability regression model based on each type of unit is as follows: where k is permeability (10 −3 μm 2 ) and ∅ is porosity. e optimal exponential and power law relationship functions were used in each type of hydraulic unit system to predict the permeability based on the modeling software platform. And then, the 3D-geological body of predicted permeability was constructed.  Figure 7(a)) was consistent with the provenance direction. e 1 st layer of the geological body was randomly selected to draw the relationship chart of the porosity-permeability-facies, from Figure 8, and the porosity and permeability were interrelated in three dimensions. Random point selection was performed in the interpolation area between the model wells, and the distribution law of the permeability value was consistent with the geological setting. According to the plane distribution of the facies, it can be seen that the high permeability area matched the types of microfacies with good physical properties such as channel bar and point bar.

Validation of the Permeability
e distribution characteristics of the prediction model on the plane were reasonable except for a few extreme values that were not controlled by the well points.  (Figure 9). However, it is worth noting that when the curve passes through the origin, it jumps directly to a    high value interval above 0.8000. e jump-like fluctuation of the semivariance value at the origin reflects the heterogeneity which was also consistent with the actual strong heterogeneity of the target reservoir on the vertical direction. At the same time, the wells D52 and D30 were selected to randomly eliminate, and the previous modeling ideas were repeated to reconstruct the model and output the predicted permeability profiles on the trajectories of these two wells ( Figure 10). Comparing the predicted permeability data with the core test data of the two wells, it can be seen that the repeated action did not interfere with the accuracy of the model, and the predicted permeability information of the two wells were basically consistent with the core data. is comparison also verified the reliability of the model information database.

Results and Discussion
e heterogeneity of tight gas reservoirs is usually serious, and conventional prediction methods often cannot accurately reflect the permeability characteristics of the reservoir. In the study of Izadi et al., HFU model was applied to a lowpermeability reservoir based on a series core data, and the model identified more flow units and also predicted moreaccurate permeability/porosity relationships [34]. Lin et al. conducted a tight reservoir permeability prediction based on an artificial neural network for a small sample environment, and the prediction results were good [35]. Unlike the above neural network model, the permeability was used as an indirect output parameter in this study, and the BP network with higher accuracy was trained to predict the types of hydraulic units. e three-dimension of the hydraulic unit was established by the sequential indicator simulation technology, and on the basis of the geological body, the permeability mode grids were fully calculated based on the porosity-permeability relationship of each unit. e model was controlled by high precision hydraulic unit on the vertical and sedimentary facies on the plane. e final permeability prediction system was more reliable, which could provide a powerful force for next numerical simulation work and related evaluation work.

Conclusions
e hydrodynamic and geophysical characteristics of tight gas reservoirs were fully considered by this study, and its application effect was better. e conclusion is as follows: (1) e flow zone indicator and the cluster method of the sum of square deviations were selected to classify the reservoir hydraulic units into 4 types. (2) e best network was trained; in the results of statistical frequency of HU prediction, the overall prediction effect was close to 85%, and the prediction accuracy of HU#4 reached 90.48%. (3) By calculating the variogram of the permeability on the plane and vertical, the reliability of the predicted geological body was verified; the rationality of the model space was also reflected by the analysis results of the porosity-permeability-facies and the repeated action. (4) e facies-controlled modeling method for permeability based on hydrodynamic and geophysics characteristics can improve prediction accuracy. It provides an effective method for spatial permeability prediction using sedimentary facies and logging data.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e author declares that there are no conflicts of interest.