Yingqiong basin is a proven hydrocarbon-rich basin in South China Sea. There are a number of large exploration prospects in high temperature and over-pressured formations, especially in Yacheng Block of Qiongdongnan basin and Dongfang District of Yinggehai Basin. Owing to good exploration situation, we have already achieved proven geological reserves over 1000 × 10 8 m 3. In recent years, a few drilled HPHT wells have confirmed that pressure predicted by conventional method was wildly inaccurate. From the view of regional stress, the accuracy of the pressure prediction will be substantially improved. Accurate pressure prediction and three-dimensional pressure modeling which are based on three-dimensional lithology modeling are the cornerstone to achieve exploration breakthrough. In this paper, the use of the triple constraint trend lithology model broke through the traditional method of seismic lithology prediction only by means of impedance threshold value. Compared with actual data and prediction, it confirms that three-dimensional pressure modeling method is reasonable and effective, and has a wide prospect of application.
According to relevant references, more than two thirds of the world’s sedimentary basins develop at different levels of overpressure, among which hydrocarbon distribution of about 160 sedimentary basins is closely genetically related with overpressure body. Large and medium-sized gas fields have been found in these overpressure basins, which indicates that overpressure is closely related to hydrocarbon generation, migration, accumulation and distribution and so on. Yinggehai-Qiongdongnan Cenozoic sedimentary basins were rapidly subsiding basins, which were characterized by abnormal high temperature and high pressure. Abnormal high pressure is mainly due to the late rapid subsidence and under-compaction mudstone formed by rapid deposition, which relies on steady subsidence and weak tectonic activity to maintain [
Over the years, old pressure prediction method leads to reservoir geologic evaluation problems and potential safety hazard, and even seriously restricts the HPHT exploration breakthrough in this field. The error of forecast will be reduced from the perspective of regional stress magnitude and direction. It’s necessary to construct a three-dimensional pressure model mainly controlled by lithology model but the traditional forecasting method is only based on lithologic seismic impedance threshold value [
Structural model is the basis for the establishment of three-dimensional model, which reflects the formation of spatial framework and plays constrained role in phase model and property model. Structural model consists of interval model and fault model, which based mainly on structural surface owing to lack of faults. In the need of subsequent cap quality evaluation, mudstones can be divided into 4 layers based on structural interpretation of Huangliu formation. Flat grid spacing is set to 100 m * 100 m with 105 layers of vertical grid. Specific longitudinal layer is divided as follows: Four mudstones are successively divided into 20, 10, 10, 15 layers; the lower part of the reservoirs is divided into 50 layers with 2.2 m single average thickness.
During the first section of Huangliu formation (T30 - T31), the center of basin is conducive to the development of gravity flow. With transgression intensified, the scope of the deposition significantly expanded, the center of basin received a large reservoir system owing to surrounding active source system. Seismic facies have a transition from imbricate structure to chaotic characteristics, reflecting that the west of low delta directly connected to the submarine fan in the east. Main waterways have obviously undercut phenomenon in three-dimensional seismic profile, indicating rapid strong erosion with similar shape of submarine fan at home and abroad.
West of Dongfang a structure mainly developed delta front and fore-delta, up-leaning lithologic traps was formed by HST overlying mudstone and low fore-delta mudstone. During the first section of Huangliu formation, offshore beach and sand bars developed from Hainan provenance in the east. Inverted layered reservoir of III gas group in the thin edge of DFA well indicates hydrodynamic orientation. It’s predicted that higher permeability reservoirs may be in higher parts and thicker sand dam body.
On the basis of sedimentary facies and its digital boundary [
A number of favorable geological modeling spatial prediction information can be obtained from seismic attribute data. Numerous seismic attributes are through the same seismic trace, but from different perspectives of understanding and interpretation. Multiplicity of seismic attributes can cause abnormal response owing to various fluids, reservoir properties and depositional environment. Combined with the actual drilling condition and sedimentary facies, conventional seismic attributes can better reflect sedimentary reservoir changes in the target area.
Petro-physical studied on four wells shows that: P-impedance properties can better distinguish between sandstone and mudstone. Mudstone shows high impedance, mainly concentrates in more than 8,000,000 kg/m3 * m/s region, while sandstone shows low impedance, the main concentrates in an area of less than 8,000,000 kg/m3 * m/s. Therefore, wave impedance inversion can basically describe the variation between different wells. On the basis of the structural model, wave impedance inversion model will be obtained from wave impedance inversion results to grid model (
Three general ways of obtaining P-wave velocity: pre-stack simultaneous inversion, stack velocity from seismic data processing and then calculate longitudinal wave velocity through DIX formula, by using other seismic attributes with better correlation. Pre-stack simultaneous inversion method is not applicable in this area because of all from post-stack data. In this study, we use two methods to establish P-wave velocity model: One is obtained by DIX formula, which input discrete point of velocity, then sampled in the geological model and interpolated to form a layer velocity model [
Compared with these two methods, we can find that it’s difficult to explain partial low-speed or high-speed stratigraphy owing to taking no consideration of time-depth conversion and velocity prediction by pressure, as well as limitation of velocity resolution. Thus, changes in the processing velocity is relatively stable, in terms of velocity obtained by wave impedance can reflect high and low
Principle of lithology model is to establish lithology between wells and predict reservoir structure by using stochastic modeling method and a variety of constraints. The accuracy of lithology modeling on the one hand depends on geological accuracy, on the other hand depends on constraint conditions [
From the relationship between impedance of drilled wells and lithology, the impedance of the target area can well reflect lithology. Therefore, this prediction shows high accuracy. From the point of well logging and integrated lithology interpretation, lithology can be classified as two types (sandstone and mudstone) in order to make better use of impedance data. Specific procedure of the actual modeling process is as follows:
1) Roughen lithology of four drilled wells according to model grid in target area;
2) Sandstone content contour map can be made based on geo-statistical analysis of lithology, vertical ratio and variogram. It can serve as trend constraint of lithology model;
3) Statistical correspondence analysis between impedance characteristics and different lithology can be spatial constraints of lithology model;
4) Under the above constrained conditions, establish multiple lithology models by stochastic simulation algorithm and optimize model based on geological knowledge.
Establish lithology model of Dongfang district by the constraint of longitudinal, plane and three dimensional phase (
Four pore pressure models are established in this study. Pore pressure prediction model is set up by using velocity and seismic information in less well areas, while for relatively more wells areas pressure coefficient of drilled wells are more
commonly used. Four models can be used to predict the formation pore pressure.
Combined with acoustic travel time curves and MTD data, pore pressure prediction model is summarized to represent the relationship between P-wave velocity difference and pressure coefficient [
According to the analysis of prediction error, pressure predicted by wave impedance has high accuracy with actual pressure [
Different from pore pressure prediction models, solid model of pore pressure is obtained from interpolated trend control. Based on pressure coefficient of a single well, solid model of pore pressure is updated from time to time by newly drilled wells with seismic attributes as lateral trends. For Dongfang district, solid model of pore pressure is applied to data of good quality, more accurate velocity, pressure coefficient and interval velocity associated with good relativity.
1) Three-dimensional structural model, phase model, wave impedance model are the basis of velocity model and lithology model. Velocity model can be used to establish depth conversion and pressure prediction model. Lithology model based on longitudinal, plane and wave impedance constrained trends. Three- dimensional stress field model is established by regional stress analysis.
2) Pore pressure model is divided into pore pressure prediction model and pore pressure solid model, among which pore pressure prediction model has good predictive capabilities of formation pore pressure; pore pressure solid
model is much closer to the true formation pore pressure with constantly updated data.
In conclusion, more accurate forecast pore pressure before drilling can not only guide well construction design, but also greatly improve drilling safety and save drilling costs.
Liu, A.Q., Zhu, P.Y., Pan, G.C., Fan, C.W., Liu, B. and Wu, Y.P. (2017) Three-Dimensional Pressure Modeling of South China Sea in High Tem- perature High Pressure Field. Open Journal of Marine Science, 7, 271-280. https://doi.org/10.4236/ojms.2017.72019