Problems and Solutions of Velocity Modeling in Natural Gas Exploration

Show more

1. Introduction

The west coast of

In this paper, three typical cases of velocity forecast and model building in Ying- Qiong basin are analyzed. And every case raises an unusual but pressing problem. There are three main problems: how to do velocity reprocessing and determine its accuracy, the system error existing between logging VSP velocity and velocity spectrum interpretation, and larger error occurs when predicting depth of new drilling well with the time-distance relationship of drilled wells nearby. Combined with problems above and the misunderstanding, method research is carried out and corresponding solution is put forward for every problem, which lays a solid foundation for accurate velocity prediction before drilling.

2. Misunderstanding in the Velocity Prediction Process of Conventional Natural Gas Exploration

As is known to all, velocity prediction is key and difficulty for the geophysical work. Conventionally, stacking velocity is obtained by velocity spectrum interpretation, and the accuracy of interpretation is judged by flattening effect of phase axis. Under the premise of formation with no large dip angle, the stacking velocity is equivalent to root mean square velocity. Meanwhile, when the seismic data is in good quality, velocity conversion formula (e.g. the DIX formula) can be used to convert root mean square velocity into interval velocity or the average velocity which can be used in pressure prediction and time-depth conversion [5] [6] [7] . But in practice, the assumptions that formation has no large dip angle and seismic data is in good quality are difficult to meet, and there also exists misunderstanding of velocity prediction [8] . So velocity prediction tends to have larger error.

1) Direct application of imaging velocity in time-depth conversion

This is the first misunderstanding of velocity prediction. Gathers data of the same velocity spectrum point are chosen from the deepwater area in the Qiong Dong Nan basin. Three different interpretations are made, including lesser velocity, medium velocity, and biggish velocity (Figure 1). Stacking velocity is obtained by velocity spectrum interpretation, and the accuracy of interpretation is judged by flattening effect of phase axis. We can see that the flattening effect of phase axes in three different solutions is basically consistent, but large errors (about 10% between low/high velocity and medium velocity) occur. So it is problematic to determine accuracy of velocity interpretation only using flattening effect of gathers data or imaging quality. Hereby, we come to the conclusion that the imaging velocity obtained through seismic processing usually cannot be directly used to make pressure prediction and time-depth conversion [9] .

2) The establishment of velocity volume can improve time-depth conversion accuracy

Many people believe that establishing velocity volume can improve the accuracy of the time-depth conversion. The idea is very wrong. The velocity volume we get through

Figure 1. Three different interpretation schemes for the same velocity spectrum point.

whatever methods and software is based on single point gathers and velocity interpretation, and takes reference of the lateral variation trends of layers and longitudinal change trends of VSP data on the well site, and is obtained by the method of interpolation extrapolation. In this way, velocity accuracy of single-point velocity gathers, and velocity volume accuracy still depends upon single point velocity interpretation precision. Only on the premise of higher single point velocity interpretation precision could the data be used to build velocity volume and extract information of single point, layer and volume [10] [11] [12] (Figure 2).

3. Typical Case Analysis

1) Target A in Dong Fang zone of Ying Ge Hai basin

Figure 2. Regional reservoir forming pattern map of target A.

Target A which is located in the central part of Ying Ge Hai basin and near Dong Fang 1-1 gas field belongs to medium deep favorable targets. The main development risk is high pressure. Before pressure prediction, we usually first carefully interpret the velocity spectrum on the design well site or nearby and then make pressure prediction combining with the typical pressure prediction model [13] . In this case, as is shown in the left part of Figure 3, there is an obvious influence of multiple waves and energy clusters on target layers are not concentrated. It leads to more uncertainties of velocity interpretation and cannot meet the requirements of target fine pressure prediction.

In order to solve this problem, the data is reprocessed (Figure 3 right). With more concentrated energy clusters on target layers and less multiple waves, uniqueness of velocity interpretation is better [14] [15] [16] . However, there is a certain velocity error before and after reprocessing. And comparing the two pressure in problem that how to do velocity reprocessing and determine its accuracy.

2) Target B in Bao Dao zone of Qiong Dong Nan basin

Target B is in the Bao Dao zone northeastern of Qiong Dong Nan basin, and located between Bao Dao sag and Song Tao bulge. Before drilling the target, we concern on time-depth conversion problem. Take a drilled well as an example: there is obvious phenomenon of scissors difference between the average velocity calculated by VSP and that obtained from velocity spectrum interpretation. For the tectonic location of target B is similar with this drilled well, the average velocity of Target B obtained from velocity spectrum interpretation is adjusted according to the correction of this well.

After target B is drilled, VSP logging velocity is projected onto the velocity spectrum (Figure 4 left), and it is found that phenomenon of scissors difference between the average velocity calculated by VSP and that obtained from velocity spectrum interpretation also appears (Figure 4 left). The unadjusted velocity has large error compared to

Figure 3. Velocity spectrum comparison before and after reprocessing.

Figure 4. Velocity comparison of target B before and after drilling.

measured velocity, and though the adjusted velocity is closer, there is still error. Hereby, this example raises the second main problem that the system error exists between logging VSP velocity and velocity spectrum interpretation.

3) Larger error occur when predicting depth of new drilling well with the time-dis- tance relation of drilled wells nearby

Recently HK29B well is drilled in

BD

Figure 5. Typical seismic profiles.

4. Solution and Method Discussion

1) Processing velocity with geological model for guidance

In order to answer the first main question, the method that velocity should be processed with geological model for guidance is proposed and a typical example in YGH basin is given to illustrate the idea (Figure 6). There is an obvious influence of multiple waves for medium-deep formation and effective energy clusters are difficult to concentrate. So the velocity spectrum should be reprocessed. As is shown in the upper picture of Figure 6, the processing outcome is alright when just using radon filtering without considering geological model. However, big difference happens when compared with regional velocity data. The drilled well reveals that a velocity inflection point appears around 1200 ms, but it doesn’t appear in this situation. As is shown in the lower picture of Figure 6, when geological model is taken for guidance, the inflection point appears around 1200 ms and the velocity spectrum changing trends are basically identical with that of logging velocity.

Figure 6. Effect diagram of velocity processed with geological model for guidance.

2) Establishing regional stratigraphic statics velocity correction templates

As for the second problem, we can see from the second case above that there is phenomenon of scissors difference between the average velocity calculated by VSP and that obtained from velocity spectrum interpretation in different structures of QDN basin. And the differences vary from one block to another. Take QDN basin as an example. To solve this kind of problem, regional stratigraphic statics velocity correction templates of every block in this area are established according to VSP velocity of drilled wells and velocity spectrum crossing the well points. Therefore, the velocity corrected by these templates will be closer to the actual velocity when making new velocity spectrum interpretation in this area in the future.

3) Using along-horizon velocity interpretation and establishing forward modeling in assistance

The velocity spectrum is restricted by its own accuracy, especially when the target layer velocity is abnormally low or high pressure exists. To deal with the third difficulty, two solutions are given. First, project the seismic interpretation horizon in the form of normal velocity into velocity spectrum interpretation work area, and pick up velocity along the horizon (Figure 7). In this way, resulting interval velocity has a certain geological meaning and vertical resolution of the velocity can be increased by interpreting horizon more intensive. Second, build a two-dimensional forward model to assist velocity interpretation. The forward model is based on typical seismic profiles which integrate seismic interpretation horizon with results of along-horizon velocity interpretation. It analyses the tectonic form and amplitude of the target layers in time domain, depth domain and velocity domain to guide the velocity interpretation and abnormal velocity judgment.

Figure 7. Velocity pick-up along the horizon.

4) Building the high precision velocity volume with velocity model corrected by information of drilled wells, horizons and geologic zonation

To solve the problem that regional velocity errors exist between new wells and drilled wells at the same structure around the diapir zone, the thinking of building velocity volume and correcting the velocity model is proposed in two steps (Figure 8). Firstly, seismic velocity model is corrected by logging VSP or acoustic velocity curve model and difference of interval velocity is used. Secondly, velocity model obtained in the first step is matched by geological zonation and time zonation, and ratio of interval velocity is adopted. During the process, three quality control points should be focused on. First, make sure the single well VSP or acoustic velocity reasonable, eliminate abnormal data and keep a vertical change tendency as smooth as possible. Second, pay attention to matching of geological layers and time horizons and their closed relationship. Third, compare velocity model error before and after correction, and find out whether the error is smaller and basically in a normal distribution.

5. Application Analysis

Velocity model of Dong Fang diapir area in Ying Ge Hai basin is established based on the method above. First, the lateral and vertical variation characteristics of the interval velocity are obvious. The velocity is relatively low at the center of diapir area as for the high pressure, which is consistent with geological knowledge. Then, along-horizon extracted interval velocity is similar to the average velocity lateral variation characteristics and time structure feature. Velocity is significantly less than that at the diapir peripheral, and the interval velocity at the location of DF

Figure 8. Flow chart of velocity model correction.

Figure 9. Error comparison histogram before and after velocity correction.

6. Conclusion and Understanding

It’s very difficult to predict the velocity of favorable targets in a basin with low overall prospecting degree like Ying-Qiong basin. Combined with problems and the misunderstanding, in this paper method research is carried out and corresponding solution is put forward for every problem, which lays a solid foundation for accurate velocity prediction before drilling. First, processing velocity with geological model for guidance solves the problem how to do velocity reprocessing and determine its accuracy. Then, establishing regional stratigraphic statics velocity correction templates diminish the system error existing between logging VSP velocity and velocity spectrum interpretation. Next, using along-horizon velocity interpretation and establishing forward modeling in assistance overcome the difficulty that velocity spectrum is restricted by its own accuracy. Finally, high precision velocity volume is established with velocity model corrected by information of drilled wells, horizons and geologic zonation.

In short, as long as we know more about misunderstandings, overcome the difficulties of velocity analysis and solve the important problems in the process of gas velocity prediction, ideal velocity prediction results can be obtained.

Fund Project

National Science and technology major special projects supported by 2016ZX05024- 005.

References

[1] Chen, X.-R. (1991) An Analysis of Exploration Prospects of YingDong Structure Belt in Yinggehai Basin. Natural Gas Industry, 19, 31-34.

[2] He, J.-X., Xian, Z.-Y., Yang, X.-B. and Mao, Y.-X. (2001) Petroleum Geology and Current Exploration Trend in YingDong Slope Area, Yinggehai Basin, South China Sea. China Offshore Oil and Gas (Geology), 15, 242-248.

[3] He, J.-X., Xia, B., Zhang, S.-L. and Sun, D.-S. (2006) Migration-Accumulation of Gas under the Abnaormal High Geotemperature and Superpressure Circumstance in Yinggehai Basin. Marine Geology & Quaternary Geology, 26, 81-89.

[4] Li, L., Song, H.-B. and Yang, J.-H. (2006) A Preliminary Study of Seafloor Gas Seepage in Central Sag Zone of Yinggehai Basin. Progress in Geophysics, 21, 1244-1247.

[5] Chen, Y., Wei, Y. and Ge, Y. (2004) Multiple Reflection Analysis and Rejection in L Depression of East China Sea. China Offshore Oil and Gas (Geology), 16, 373-376.

[6] Hao, F., Dong, W.-H., Zou, H.-Y. and Yang, X.-S. (2003) Overpressure Fluid Flow and Rapid Accumulation of Natural Gas in Yinggehai Basin. Acta Petrolei Sinica, 24, 7-12.

[7] Cai, G. and Qu, Z.-Y. (2005) Study of Seismic Data Velocity and Mapping Method in Complicated Structure Area and Its Application. Natural Gas Geoscience, 16, 246-249.

[8] Liu, A.-Q., Tong, C.-X. and Li, L. (2008) Research of Velocity Analysis Method and Analysis of Velocity Influence Factors in the East Slope of Yinggehai Basin. Progress in Geophysics, 23, 1909-1917.

[9] Liu, J.-P., Yang, Y.-Q. and Li, X.-Z. (1999) The Method for Extracting Velocity Under the Complex Topographic Condition and its Analysis. Geophysical and Geochemical Exploration, 23, 259-264.

[10] Liu, Q.-F. (2003) Analysis and Application of the Velocity in Imaging of Seismic Data and Mapping of Structure. Journal of Jianghan Petroleum University of Staff and Workers, 17, 52-54.

[11] Liu, T.-Y., Li, Y.-S. and Shi, X.-M. (1995) Geopressure Prediction and Evaluation Methods and its Application in Yinggehai Basin. China Offshore Oil and Gas (Geology), 9, 333-339.

[12] He, J.-X., Li, M.-X. and Huang, B.-J. (2000) The Analysis of the Oil and Gas Exploration Prospect and the Distribution of Outflow of Oil and Gas in the Northern Slope of the Yinggehai Basin. Natural Gas Geoscience, 11, 1-9.

[13] Meng, Q.-S., Chu, X.-F., Guo, X.-J., Fan, Y.-Q. and Jia, Y.-G. (2007) The Application of High Resolution Seismic Data Processing Technique in Multi-Channel Shallow Offshore Engineering Seismic Surveys. Progress in Geophysics, No. 3, 1006-1010.

[14] Tian, S.-C., Chen, Y.-J. and Shi, F.-C. (2004) The Role of Abnormal Pressure Compartment in Oil and Gas Accumulation. Earth Science Frontiers, 11, 283-284.

[15] Xie, R.-J., Zhu, G.-S. and Qi, J.-F. (2003) Application of Sonic Logging in Average Velocity Field. Journal of Southwest Petroleum Institute, 25, 9-12.

[16] Wang, L.-M. and Li, Q.-C. (2006) Improving the Accuracy of Converted-Wave Velocity Analysis with Selected-Correlation Method. Progress in Geophysics, 21, 1213-1220.