Bieng Zih Hsieh

Lithology identification of aquifers from geophysical well logs and fuzzy logic analysis: Shui-Lin Area, Taiwan

The purpose of this study is to construct a fuzzy lithology system from well logs to identify formation lithology of a groundwater aquifer system in order to better apply conventional well logging interpretation in hydro-geologic studies because well log responses of aquifers are sometimes different from those of conventional oil and gas reservoirs. The input variables for this system are the gamma-ray log reading, the separation between the spherically focused resistivity and the deep very-enhanced resistivity curves, and the borehole compensated sonic log reading. The output variable is groundwater formation lithology. All linguistic variables are based on five linguistic terms with a trapezoidal membership function. In this study, 50 data sets are clustered into 40 training sets and 10 testing sets for constructing the fuzzy lithology system and validating the ability of system prediction, respectively. The rule-based database containing 12 fuzzy lithology rules is developed from the training data sets, and the rule strength is weighted. A Madani inference system and the bisector of area defuzzification method are used for fuzzy inference and defuzzification. The success of training performance and the prediction ability were both 90%, with the calculated correlation of training and testing equal to 0.925 and 0.928, respectively. Well logs and core data from a clastic aquifer (depths 100-198m) in the Shui-Lin area of west-central Taiwan are used for testing the systems construction. Comparison of results from core analysis, well logging and the fuzzy lithology system indicates that even though the well logging method can easily define a permeable sand formation, distinguishing between silts and sands and determining grain size variation in sands is more subjective. These shortcomings can be improved by a fuzzy lithology system that is able to yield more objective decisions than some conventional methods of log interpretation.

Identification of stratigraphic formation interfaces using wavelet and Fourier transforms

The purpose of this study was to identify the formation interfaces from geophysical well log data using the wavelet transform, and a combination of the wavelet transform and the Fourier transform methods. In the wavelet transform method, the identification of formation interfaces is based on the wavelet coefficients from the wavelet transform of spontaneous potential (SP) log and gamma ray (GR) log data. In the combination of the wavelet transform and the Fourier transform methods, the wavelet transform, spectrum analysis, and logarithmic transform of well logs were applied to the SP and GR log data successively to obtain clear signals for identifying the stratigraphic formation interface. In this study, a set of ideal log data was first created and analyzed to test the validity of the developed procedures. In analyzing the SP and GR logs from a field, both the wavelet transform method and conventional well log analysis showed similar results. The results from a combination of the wavelet transform and the Fourier transform methods, however, were better than those from the wavelet transform method and the conventional well log analysis.

Propagation of Radius of Investigation from Producing Well

Abstract The purpose of this study is to estimate the pressure disturbance area, or the propagation of the radius of investigation, from a producing well in an infinite reservoir by using both analytical and numerical methods. A linear (or radius) coefficient in the relationship between the square of the dimensionless radius of investigation and the dimensionless time is studied and derived. The radius coefficient in the equation is a constant, and varies with different criteria of radius of investigation defined, i.e., the amount of pressure change from the initial formation pressure at the pressure front of the pressure disturbance area. As the dimensionless pressure defined at the pressure front changed from 0.1095 to 10−9, the radius coefficient varied from 4 to 71.15. The radius coefficient is independent of the level of the flow rate for a well producing at a constant flow rate. For a well producing with variable flow rates, the radius coefficient is not a constant for the case of larger pressure drops defined at the pressure front. The skin factor does not affect the result of the calculated radius of investigation. The wellbore storage volume will affect the propagation of the radius of investigation only at an early time, depending on the wellbore storage volume.

Estimation of groundwater aquifer formation-strength parameters from geophysical well logs: The southwestern coastal area of Yun-Lin, Taiwan

Abstract The purpose of this study is to estimate groundwater aquifer formation-strength parameters including shear modulus, bulk modulus, Poissons ratio, and Youngs modulus by using geophysical well logs. A new dispersed-shale index equation was developed by using the natural gamma-ray log and the compensated formation density log to solve a confusing problem of the compaction factor setting in the calculation of sonic porosity for an unconsolidated groundwater aquifer. A useful Poissons ratio estimation method was employed to estimate groundwater aquifer formation-strength parameters when shear-wave transit time data is lacking in groundwater wells. Hydrogeologic parameters are characterized in estimating formation-strength parameters. Five wells in the southwestern coastal area of Yun-Lin, Taiwan, were logged, and four shallow aquifers were identified from log-derived hydrogeologic characteristics less than 200 m in depth. The formation-strength parameters for aquifers between 310 and 500 m in depth were calculated in two wells because complete formation density and compressional-wave transit time data were available. The results of the aquifers formation-strength parameters demonstrate that both shear modulus, ranging from 0.15 to 0.42 * 106 psi, and Youngs modulus, ranging from 0.40 to 1.07 * 106 psi, increase with depth, whereas bulk compressibility, ranging from 1.2 to 2.6 * 10−6 psi−1, decreases with increasing depth.

Determination of the constant coefficient in pressure propagation equation

Abstract Many studies on the pressure propagation in porous media show that the relationship between the square of the radius of investigation and the producing time of a well is linear. The linear relationship is the pressure propagation equation (or pressure equation). However, these studies do not arrive at the same constant coefficient representing the linear relationship. The purpose of this study is to find the constant coefficient for the pressure propagation equation. The constant coefficient from this study is derived from analyzing pressure transient behavior and observing the boundary effect time (the starting time of wellbore pressure affected by the reservoir boundary). The results show that the constant coefficient in the pressure equation, i.e., the linear constant between the square of dimensionless radius of investigation and the dimensionless producing time is 17.82.

A prediction of pressure behavior with production data only

Abstract The purpose of this study is to derive and calculate a pressure response function in order to predict the pressure behavior of a producing or testing well without knowledge of reservoir parameters. The equation of pressure response function is theoretically derived from an analytical solution of the diffusivity flow equation. The pressure response function for a specific reservoir can be derived from the equation of the pressure response function with known production data, such as wellbore pressure and flow rate information. Once calculated, the pressure response function can be used to calculate future wellbore pressures for given flow rates. In these calculations or predictions, no reservoir and fluid parameters are necessary. The derived equations and the calculation methodology have been verified in this study by using reservoir simulation data. The ability to predict future wellbore pressure behavior from a pressure response function is also demonstrated.

Estimation of Water Saturation in Rocks with Conductive Matrix

Abstract The purposes of this study are to develop a water saturation model (MCRM model) for a formation with conductive rock matrix and to calculate the empirical parameters used in the model by using the optimization method. The empirical parameters in the new model (MCRM model) include the conductive formation parameter (b′), tortuosity factor (a), cementation exponent (m), and saturation exponent (n). A field data is used to illustrate the procedures to derive the empirical parameters by the optimization method. The results of the empirical parameters from three models, the MCRM model, the Givens model, and the Archie equation (model), are compared. The conductive formation parameter (b′) for the field data analyzed is less then one, indicating that the formation contains conductive rock matrix.