The tight sandstone reservoir of the Cretaceous Bashijiqike Formation in Dabei Area of Kuqa Depression in Xinjiang Tarim Oilfield is buried below 6000 m; its porosity is generally between 3% and 8%, and its permeability is between 0.01 and 0.1 mD; it belongs to a deep tight sandstone reservoir  , where micro-fractures are developed, this enhances the heterogeneity of the reservoir; it is more difficult to evaluate the fluid property in the reservoir by well logging. Therefore, in addition to choose a downhole detecting instruments with higher detection precision, studies should be carried out on various logging interpretation methods.
Sulige Gas Field in China was used as example by Li Xia (2013), methods for identifying tight sandstone gas reservoirs were denoted based on conventional logging data, these methods included the contrast method of equivalent P-wave elasticity modulus, A-K intersection combination method, the intersection of P-wave and S-wave ratios with Poisson’s ratio and the method of Poisson’s ratio and compression coefficient ratio  . But in the target zone of 8th Member of Shihezi Formation, the average value of porosity is 8.3%, in Kuqa Area with lower porosity; its adaptability remains to be verified. Yang Shuangding (2005) introduced the subtraction method based on NMR logging data and the method of shift spectrum for identifying tight sandstone gas reservoirs in Ordos Basin  , but in the article, only the response of pure gas zones was discussed, for the effect of identifying water zone in the tight sandstone was difficult to be determined by using the method. Zhang Yongjun (2012) discussed in detail the methods of using array acoustic logging data for identifying physical foundation, parameter extraction and identifying the reservoir fluid in the tight sandstone gas reservoirs  . Guo Zhenhua (2010) applied the characteristic parameters of neutron, resistivity and density logging curves to establish the method for qualitatively evaluate the grade of open flow in the gas zone, certain effect was obtained in the 2nd Member of Shihezi Formation in Daniudi Gas Field  . Zhang Haitao (2012) introduced the method of using logging curve for identifying the types of diageneticfacies  , which expanded the application of logging data in geologic sphere. The methods described by the above scholars have included the current technical means for identifying gas reservoirs by using logging data, in which the methods for logging identification of water zone and gas-water coexistent zone are not discussed in depth, the detecting precision and applicability of the methods for identifying the tight sandstone reservoirs remain to be verified. Li Jun (2009) found that the error in calculating saturation is 15% - 30% by using an indequate m exponent because of the influences of low porosity  . Li Yuegang (2015) proposed a new method for calculation gas saturation by variable rock-eletro parameters  . But both of them do not carry out special researches on the methods for directly calculating the saturation in the reservoirs with low porosity and low permeability.
In allusion to the issues of deep burial, higher pressure, high temperature and hard for identification of gas and water layers in the tight sandstone gas reservoirs in Dabei Area, by using the logging data as the major means and on the basis of fine evaluation of reservoir fractures, a novel approach for calculating saturation is established in consideration of the influential factors of fractures and low porosity and etc., which can meet the requirements of evaluating the saturation and precision of tight sandstone gas reservoirs in Babei Area of Tarim Basin.
2. Fracture Evaluation and Reservoir Type Division in Tight Sandstone Reservoir
Whether the fractures in tight sandstone reservoirs are grown is the key indicator for determining the value of development in the reservoir, the qualitative and quantitative depiction of fractures is the prerequisite of consequent calculation of saturation in the reservoir. Micro-fractures are generally developed in the tight sandstone reservoirs in Kuqa Depression   , but actual logging interpretation indicates that a very poor effect is caused in the application of dual induction resistivity logging evaluation. Rock mass is basically not conductive, while the parallel network for conduction in fractures is the foundation that the dual induction logging can be used for identifying fractures in carbonate rock formations, in which the resistivity in the carbonate rock formations should be more than 2000 Ω∙m at least. Resistivity of the sandstone is much lower than that of carbonate, which is only 10 - 80 Ω∙m. In the sandstone formation, the resistivity in the fractural pores is not evidently different from that of rock skeleton, which is the major cause of resulting in poor effect in fracture evaluation using the dual induction logging.
Micro-resistivity imaging logging was one of the effective methods of fracture identification  , by which the status of fracture near the wellbore could be reflected, its radial extension and connection were difficult to be discriminated. When a low frequency Stoneley wave was used to reflect the mutual seepage of fluids in the fractures or in the pores and the fluid in the wellbore, the radial detecting range could be 1.5 m      . Apparent fracture porosity obtained from the micro-resistivity imaging logging is intersected with the permeability of Stoneley wave, by which the connectivity of fractures can be better evaluated. If there were no logging data available, the apparent fracture porosity obtained from the difference between density logging and acoustic logging could be used for reservoir evaluation, but poorer effect is obtained in the evaluation  . Based on the selection of acquisition of series logging data in Dabei Area, the micro-resistivity imaging logging is combined with the full-wave acoustic logging in the study for an effectiveness evaluation of reservoir fractures.
Figure 1 is the result of fractural evaluation by using the logging data from electric imaging logging and full-wave acoustic logging in Well Keshen801. It can be seen in the Figure that at 7053 - 7114 m the permeability from Stoneley wave (Track 5) and the reflecting coefficient (Track 8) are high values, which indicates that there exists a good connectivity in formation. At the same time, the parameters of apparent fracture porosity, fracture length and width and its density(Track 10) obtained from electric imaging logging data indicate that fractures are developed in the layer, the identification of acoustic logging shows that there exists good consistence in the fractures. In the Figure, the apparent fracture
Figure 1. The result of fractural evaluation by using the data of electric imaging logging and full-wave acoustic logging in Well Keshen 801.
porosity in Track 4 is obtained from the porosity difference calculated from the curves of density logging and sonic velocity logging  , it is commonly corresponded with the fracture in imaging logging, it also indicates that in deep tight sandstone formations, the combination of full-wave acoustic logging with micro-resistivity imaging logging should be the first choice for evaluating the effectiveness of fractures, while the conventional well logging can be a supplement.
For the evaluation of fracture development by taking the reservoir as a characteristic scale, the reservoir characteristic parameters from acoustic and electric loggings are output based on the reservoir. The tight sandstone reservoir types in Dabei Area are classified based on the data of mercury intrusion capillary curve and thin slice analysis. The reservoirs are divided into 4 types, Type I: (permeability) > 9%, K (porosity) > 1 mD; Type II: is 6% - 9%, K is 0.1 - 1 mD; Type III: is 3.5% - 6%, K is 0.055 - 0.11 mD; Dry Layer: < 3.5%, K < 0.055 mD. In Figure 1, the reservoir class in Track 7 is identified according to the
Figure 2. The cross-plot of permeability from Stoneley permeability and apparent fracture porosity in a wellblock of Dabei Area (according to 136 formations data from 8 wells, provided by Tarim Oilfield Company).
above criterion. Figure 2 is the cross plot of Stoneley wave permeability and apparent fracture porosity in a wellblock of Kuqa Area (here the apparent fracture porosity is obtained from the treatment of micro-resistivity imaging logging through core calibration). Based on Figure 2, by taking the apparent fracture porosity of 0.055% as the limit, the reservoirs are divided into fracture developed and fracture less developed reservoirs. By taking Stoneley wave calculated permeability of 0.2 mD as the limit, the reservoir connectivity is divided into a good or bad one. Type I, II and III reservoirs are mostly with good connectivity, the Dry Layer and part of Type III reservoirs are with poor connectivity, they are no productivity. In some of the Type II and III reservoirs, fractures are less developed and are lack of productivity.
3. A Novel Method for Calculating Matrix Saturation in Consideration of Low Porosity and Low Permeability
In a common resistivity logging, fluid type in reservoir can be qualitatively distinguished according to the characteristics of formation invasion, but this method is only suitable for the formations with conventional porosity and permeability. In the tight sandstone reservoirs in deep formation of Dabei Area, permeability is poor, at the same time, as it is influenced by fractures, the curve between the deep induction resistivity and shallow induction resistivity is not very different, the fluid is impossibly identified based on the invasion characteristics (it is shown as the curve in Track III of Figure 1). For a quantitative identification of reservoir fluid and for satisfying the requirement of accurate depiction of geological reserves, the establishment of saturation in tight sandstone reservoirs in Dabei Area is divided into 2 steps, a saturation model which takes account of low porosity and low permeability and is not affected by fractures established first and then takes account of the effect of fractures. The matrix porosity is the key factor for deciding the reserves in the reservoirs, therefore in this chapter, the method for calculating the matrix saturation of tight sandstone reservoirs in Dabei Area is emphatically discussed.
For sandstone formation, it is generally believable that Archie Theory and its improved method are used for saturation calculations      , but a non-Archie phenomenon is often induced in a low porosity and low permeability reservoirs  . As is shown in Figure 3, in a core test conducted in a wellblock of Babei Area, the resistivity increase coefficient is not in a single linear relation with water saturation and it is divided into 4 kinds of relations. Therefore the model for calculating the saturation of tight sandstone reservoirs should be reconsidered.
Figure 3. The relation between the resistivity increase coefficient and water saturation test (based on the data from 365 tests at 25˚C and 27 MPa, provided by Tarim Oilfield Company).
According to the Archie theory  , the relationship between the resistivity increase coefficient and water saturation is:
where in Equation (1), the exponent of saturation n is a given value, that is in the cross plot of resistivity increase coefficient and water saturation in a double logarithmic coordinate, the fitting relationship between the both is in a straight line, but there are roughly 4 relationships shown in Figure 3. The examination from Relation 1 to Relation 4 indicates that the porosity is raised gradually, which is consistent with the variation trend of resistivity increase coefficient in Archie Phenomenon by Montaron   . In allusion to the influence of arrangement of rock particles on porosity, a calculation equation based on probability statistics is introduced  .
As P only represents the conductive formation water in the rock pores, it is simply called a conductive pore water, thus , if , , then Equation (2) is turned into
If there exists a rock saturation, there is a rock resistivity
By combination of Equation (3) with Equation (4), a resistivity increase coefficient is obtained
According to Equation (2), , , a tendency plot of resistivity increase coefficient changing with different porosities can be derived (Figure 4). The comparison of Figure 4 and Figure 3 indicates that the shape of both of them is consistent, that is, the relationship between the resistivity increase coefficient and saturation can be established based on different porosities.
Figure 4. The relationship between resistivity increase coefficient and water saturation test ( ).
If a cementation index equation is introduced
It is combined with Equation (5), a new equation for saturation calculation is obtained:
In the equation, the effect of different porosities on the resistivity increase coefficient is taken into account, which is suitable for quantitative calculation of saturation in the reservoirs with low porosity.
4. The Method for Saturation Evaluation in Tight Sandstone Reservoirs in Dabei Area and Its Effect
In the conductive pore water, only the effects of low porosity and high irreducible water are taken into account without taking the fractures into account. By using the idea of dual porosity medium theory, the pore space volume is divided into 2 parts, which include the matrix pores and fractures. For the rock per unit volume, the volume that the rock takes is , then the volume that the matrix pores take is where because the matrix pores are smaller, it can be considered that there is no mud invasion, while because of good permeability in the fractures, it is easy for mud invasion. Based on the above rock model, the rock resistivity of matrix pores can be calculated through rock resistivity and fracture pore , the equation for calculation is :
1) The matrix saturation: in consideration of the impact of porosity and irreducible water on the matrix saturation, the equation based on conductive pore water is used.
where, the cementation index of the matrix pores can be obtained through a litho-electric experiment of small cores sampling from the area, as shown in Table 1. The of the conglomeratic sandstone is bigger than that of the other sandstones, it is not much different between the mid and fine sandstones, therefore in practical processing, the cementing indexes are established based on the lithologic property.
Table 1. The lithology-based statistics of matrix pore structure index
2) The water saturation in fractures: according to the permeability property of dual pore types, can be obtained from:
To examine the impact of pore space of fractures and etc in the reservoirs on the cementing index m, the rock core is simplified into a regular hexahedron, the length of its edge is a, a oblique crossing joint with the dip angle as and its width as breaks through the whole hexahedron, where pores and holes are developed internally. It is also simplified as a regular hexahedron with its length of edge as b, and the pores and holes are parallel with fractures (see the model shown in Figure 5). Based on the geometric relationship and Archie equation, is calculated as
Figure 5. The analytical model of pore-frature style.
According to Equation (11), make the edge length of rocks as 1, the relationship between the pore cementing index, width and dip angle of the fractures shown in Figure 6 is obtained, with the increase of its dip angle and width ( in the Figure is the fractural width) the bond index also increases.
In consideration of fractures with narrow width in the area, but their dip angles are generally larger, the bond index of fractures takes 1.2 - 1.5.
Figure 6. The relation between the pore structure index, its dip angle and width.
3) The total water saturation: the total water saturation in the formation is calculated with the following equation:
The result of well data process in a wellblock of Dabei Area indicates that when the apparent fractural pore is bigger than 0.055%, the saturation based on dual pores is 1.5% bigger than that of the one calculated with conductive pore water. Therefore the method for saturation calculation in the well block is: when the apparent fractural pore is smaller than 0.055%, the equation based on conductive pore water is used for saturation calculation. When the apparent fractural pore is bigger than 0.055%, the dual pore equation for saturation calculation is used (Figure 7). Figure 8 is the result of reservoir parameter logging process for Well Dabei102 in Dabei Area. As shown in Figure 8, the logging calculated porosity and permeability are close to the testing result in Track 6 and Track 7, at the same time, the saturation in Track 8 that is calculated in this paper is roughly equal to the result of mercury intrusion test. Table 2 is the relative error statistics between core saturaion and calculationg saturation from Well Dabei102. The error range is 0.01% - 26.04%, the average is 8.8%. Table 2 indicates that the method presented in this paper can be deployed for quantitative calculation of saturation in tight sandstone reservoirs in Dabei Area.
Figure 7. The flow diagram of the saturation calculating method and classification in Dabei Area.
Table 2. The relative error between core saturaion and calculating saturation.
Figure 8. The results of saturation calculation in Well Dabei102.
The accurate evaluation of fractures is an essential and important step in tight sandstone reservoir evaluation. The micro-resistivity logging and full wave acoustic logging can be used for a better fracture evaluation. The micro-resisti- vity logging imaging cab be used for obtaining the intersection of apparent fractural porosity and Stoneley wave permeability; the difference between porosities calculated by density and acoustic loggings can be used for obtaining the apparent fractural porosity for evaluation. If there are no data of micro-resistivity logging, the evaluation should be carried out by using the difference between porosities calculated by density and acoustic loggings, but the result is poorer. The tight sand reservoirs in Dabei Area can be divided into 4 types by stoneley permeability and apparent fractural porosity.
Influenced by physical property, the resistivity increase index is not a given value. The effects of low porosity and high irreducible water on the saturation of tight sandstone reservoirs are taken into account in the saturation model established by using conductive pore water. The RI-Sw relationship established in the study is consistent with the result of rock core testing, which indicates that the method is suitable for the evaluation of tight sandstone reservoirs.
According to the result of evaluation from Dabei Area, it indicates that when the apparent fractural porosity is less than 0.055%; the equation based on conductive pore water is used for calculation; when the apparent fractural porosity is more than 0.055%, the equation based on the theory of dual porous media is deployed. And the concrete operation flow can be seen in Figure 7. The relative error statistics between core saturaion and calculating saturation from WellDabei102 is from 0.01% to 26.04%, the average is 8.8%. The formation saturation calculated with the method is conformed to that of mercury intrusion in core testing, which further validates that the method is correct.
Supported by PetroChina Innovation Foundation (2015D-5006-0305).
is the resistivity increasing index, dimensionless;
is the formation true resistivity, Ω∙m;
is the full water saturated formation resistivity, Ω∙m;
is the water saturaton,%;
is the exponent of saturation, dimensionless;
is coefficient of saturation, dimensionless;
is the formation conductivity, 1/(Ω∙m);
is the formation water conductivity, 1/(Ω∙m);
is the volumetric ratio of unit rock fluid in the whole rock, %;
is the lower limit value of conduction in pore fluid in the unit rock, %;
is a proportional coefficient, dimensionless;
is formation water resistivity, Ω∙m;
is the volume of conductive pore water , %;
is the formation porosity, %;
is a proportional coefficient, same as , dimensionless;
is the coefficient in the formaiton resistivity Factor-Porosity relationship, dimensionless;
is the exponent in relationship: , dimensionless;
is the resistivity of matrix pores, Ω∙m;
is the porosity of matrix pores, %;
is the porosity of fracture, %;
is the mud filtrate resistivity, Ω∙m;
is the exponent in relationship: only for fracture, dimensionless;
is the exponent in relationship: only for matrix pores, dimensionless;
is the water saturaton in matrix pores,%;
is the water saturaton in fracture,%;
is the exponent of saturation for fracture, dimensionless;
is the resistivity of flushed zone, Ω∙m;
is the length of edge in the analytical model, m;
is the length of small pore in the analytical model, m;
is the incliation of fracture, radian;
is the width of fracture, m.