To fully extract and mine the multi-scale features of reservoirs and geologic structures in time/depth and space dimensions, a new 3D multi-scale volumetric curvature (MSVC) methodology is presented in this paper. W...To fully extract and mine the multi-scale features of reservoirs and geologic structures in time/depth and space dimensions, a new 3D multi-scale volumetric curvature (MSVC) methodology is presented in this paper. We also propose a fast algorithm for computing 3D volumetric curvature. In comparison to conventional volumetric curvature attributes, its main improvements and key algorithms introduce multi-frequency components expansion in time-frequency domain and the corresponding multi-scale adaptive differential operator in the wavenumber domain, into the volumetric curvature calculation. This methodology can simultaneously depict seismic multi-scale features in both time and space. Additionally, we use data fusion of volumetric curvatures at various scales to take full advantage of the geologic features and anomalies extracted by curvature measurements at different scales. The 3D MSVC can highlight geologic anomalies and reduce noise at the same time. Thus, it improves the interpretation efficiency of curvature attributes analysis. The 3D MSVC is applied to both land and marine 3D seismic data. The results demonstrate that it can indicate the spatial distribution of reservoirs, detect faults and fracture zones, and identify their multi-scale properties.展开更多
Beginning with the method of whole path iterative ray-tracing and according to the positive definiteness of the coefficient matrix of the systems of linear equations, a symmetry olock tridiagonal matrix was decomposed...Beginning with the method of whole path iterative ray-tracing and according to the positive definiteness of the coefficient matrix of the systems of linear equations, a symmetry olock tridiagonal matrix was decomposed into the product of block bidiagonal triangular matrix and its transpose by means of Cholesky decomposition. Then an algorithm for solving systems of block bidiagonal triangular linear equations was given, which is not necessary to treat with the zero elements out of banded systems. A fast algorithm for solving the systems of symmetry block tridiagonal linear equations was deduced, which can quicken the speed of ray-tracing. Finally, the simulation based on this algorithm for ray-tracing in three dimensional media was carried out. Meanwhile, the segmentally-iterative ray-tracing method and banded method for solving the systems of block tridiagonal linear equations were compared in the same model mentioned above. The convergence condition was assumed that the L-2 norm summation for mk, 1 and mk. 2 in the whole ray path was limited in 10-6. And the calculating speeds of these methods were compared. The results show that the calculating speed of this algorithm is faster than that of conventional method and the calculated results are accurate enough. In addition, its precision can be controlled according to the requirement of ray-tracing.展开更多
The characteristics of seismogenic structures are an important basis for delineating the potential seismic source areas and determining the annual occurrence rate of earthquakes. The potential seismic source area does...The characteristics of seismogenic structures are an important basis for delineating the potential seismic source areas and determining the annual occurrence rate of earthquakes. The potential seismic source area does not only have the intension that “this area has the possibility for destructive earthquakes to occur in the future" but also means that earthquakes of high magnitude interval have the characteristics of similar recurrence. When determining the seismic activity parameters of a statistical unit, some active tectonic blocks in the unit may have different background earthquakes. In order to better reflect the heterogeneity in space of seismic activities, it is necessary to divide the potential seismic source areas into three orders. By analyzing the recurrence characteristics of earthquakes of high magnitude interval in the potential source area and calculating the occurrence probability of earthquakes of high magnitude interval in the potential seismic source area in the time window for prediction, the average annual occurrence rate of earthquakes can be obtained by the method of probability equivalent conversion in the time window for prediction. This would be helpful for considering the recurrence characteristics of strong earthquakes in potential source areas within the framework of seismic risk analysis of China. Besides, the insufficient frequency of characteristic earthquakes of the next high magnitude interval in the potential source area and the heterogeneity of strong earthquakes on seismogenic structures are analyzed to see their application in seismic risk analysis.展开更多
In the process of accurate interpretation of multi-wave seismic data,we wanted to solve the problem of multi-wave information recognition.Based on techniques of elastic wave forwarding,targeting the geological model o...In the process of accurate interpretation of multi-wave seismic data,we wanted to solve the problem of multi-wave information recognition.Based on techniques of elastic wave forwarding,targeting the geological model of a reservoir of an oil field exploration area,we used a high-order staggered-grid difference technology to simulate many shots of seismic records of nonzero offset shots,implemented multi-wave seismic data processing to acquire the CMP of P waves and converted waves,NMO traces of CCP pre stacks,including AVA information and superposition profiles.Based on the AVA calculation of the model,the layer parameters of the model and the forwarding wave field relations of the P-S wave,we also compared and studied the correspondence between P waves and converted waves.The results of our analysis show that the results from simulation and from the AVO analysis are consistent.Significant wave field differences between P waves and converted waves in the same reservoir were found,which are helpful in recognizing and interpreting the multi-wave information in this area.We made use of the multi-wave data to provide the important guidelines for reservoir prediction.展开更多
基金supported by the National Natural Science Foundation of China (No. 41004054) Research Fund for the Doctoral Program of Higher Education of China (No. 20105122120002)Natural Science Key Project, Sichuan Provincial Department of Education (No. 092A011)
文摘To fully extract and mine the multi-scale features of reservoirs and geologic structures in time/depth and space dimensions, a new 3D multi-scale volumetric curvature (MSVC) methodology is presented in this paper. We also propose a fast algorithm for computing 3D volumetric curvature. In comparison to conventional volumetric curvature attributes, its main improvements and key algorithms introduce multi-frequency components expansion in time-frequency domain and the corresponding multi-scale adaptive differential operator in the wavenumber domain, into the volumetric curvature calculation. This methodology can simultaneously depict seismic multi-scale features in both time and space. Additionally, we use data fusion of volumetric curvatures at various scales to take full advantage of the geologic features and anomalies extracted by curvature measurements at different scales. The 3D MSVC can highlight geologic anomalies and reduce noise at the same time. Thus, it improves the interpretation efficiency of curvature attributes analysis. The 3D MSVC is applied to both land and marine 3D seismic data. The results demonstrate that it can indicate the spatial distribution of reservoirs, detect faults and fracture zones, and identify their multi-scale properties.
基金Project(40674071) supported by the National Natural Science Foundation of ChinaProject(KFAS2002-2003) supported by the Korea Foundation for Advanced Studies
文摘Beginning with the method of whole path iterative ray-tracing and according to the positive definiteness of the coefficient matrix of the systems of linear equations, a symmetry olock tridiagonal matrix was decomposed into the product of block bidiagonal triangular matrix and its transpose by means of Cholesky decomposition. Then an algorithm for solving systems of block bidiagonal triangular linear equations was given, which is not necessary to treat with the zero elements out of banded systems. A fast algorithm for solving the systems of symmetry block tridiagonal linear equations was deduced, which can quicken the speed of ray-tracing. Finally, the simulation based on this algorithm for ray-tracing in three dimensional media was carried out. Meanwhile, the segmentally-iterative ray-tracing method and banded method for solving the systems of block tridiagonal linear equations were compared in the same model mentioned above. The convergence condition was assumed that the L-2 norm summation for mk, 1 and mk. 2 in the whole ray path was limited in 10-6. And the calculating speeds of these methods were compared. The results show that the calculating speed of this algorithm is faster than that of conventional method and the calculated results are accurate enough. In addition, its precision can be controlled according to the requirement of ray-tracing.
文摘The characteristics of seismogenic structures are an important basis for delineating the potential seismic source areas and determining the annual occurrence rate of earthquakes. The potential seismic source area does not only have the intension that “this area has the possibility for destructive earthquakes to occur in the future" but also means that earthquakes of high magnitude interval have the characteristics of similar recurrence. When determining the seismic activity parameters of a statistical unit, some active tectonic blocks in the unit may have different background earthquakes. In order to better reflect the heterogeneity in space of seismic activities, it is necessary to divide the potential seismic source areas into three orders. By analyzing the recurrence characteristics of earthquakes of high magnitude interval in the potential source area and calculating the occurrence probability of earthquakes of high magnitude interval in the potential seismic source area in the time window for prediction, the average annual occurrence rate of earthquakes can be obtained by the method of probability equivalent conversion in the time window for prediction. This would be helpful for considering the recurrence characteristics of strong earthquakes in potential source areas within the framework of seismic risk analysis of China. Besides, the insufficient frequency of characteristic earthquakes of the next high magnitude interval in the potential source area and the heterogeneity of strong earthquakes on seismogenic structures are analyzed to see their application in seismic risk analysis.
基金the Doctor Research Fund for Universities of China (No.20070616004)the National High Technology Research and Development Program of China (No.2007AA060505)
文摘In the process of accurate interpretation of multi-wave seismic data,we wanted to solve the problem of multi-wave information recognition.Based on techniques of elastic wave forwarding,targeting the geological model of a reservoir of an oil field exploration area,we used a high-order staggered-grid difference technology to simulate many shots of seismic records of nonzero offset shots,implemented multi-wave seismic data processing to acquire the CMP of P waves and converted waves,NMO traces of CCP pre stacks,including AVA information and superposition profiles.Based on the AVA calculation of the model,the layer parameters of the model and the forwarding wave field relations of the P-S wave,we also compared and studied the correspondence between P waves and converted waves.The results of our analysis show that the results from simulation and from the AVO analysis are consistent.Significant wave field differences between P waves and converted waves in the same reservoir were found,which are helpful in recognizing and interpreting the multi-wave information in this area.We made use of the multi-wave data to provide the important guidelines for reservoir prediction.