Based on a detailed analysis of differences between seismic data and well logs, we discuss the problem of matching seismic traces and well logs and present a new matching method based on event search in instantaneous ...Based on a detailed analysis of differences between seismic data and well logs, we discuss the problem of matching seismic traces and well logs and present a new matching method based on event search in instantaneous phase which greatly improves seismic resolution. The method is based on flattening events in instantaneous phase to compare the seismic traces to the well log traces with the same phase. We calculate the coefficients using the singular value decomposition method to extrapolate the well logs. As a result, the events in the seismic profile are continuous and match well with well logs. We apply this method to the Mao-2 well in Daqing Oilfield with good results.展开更多
To test the performance of the Chinese whole-surface electromagnetic prospecting (SEP) system, system integrations, instrument performances, and large-scale production viabilities in Liaoning province and Inner Mong...To test the performance of the Chinese whole-surface electromagnetic prospecting (SEP) system, system integrations, instrument performances, and large-scale production viabilities in Liaoning province and Inner Mongolia were measured via extensive field tests. Resultant electric fields, magnetic fields, apparent resistivities, impedance phases, and inversion profiles compared favorably with results of commercial equipment from other countries. The inversion results agreed well with the geologic information from boreholes. Field tests showed that the SEP system is stable, reliable, lightweight, and easy to operate, making it suitable and ready for real-field exploration.展开更多
The conventional method of seismic data acquisition geometry design is based on the assumption of horizontal subsurface reflectors, which often is not suitable for complex structure. We start from a controlled illumin...The conventional method of seismic data acquisition geometry design is based on the assumption of horizontal subsurface reflectors, which often is not suitable for complex structure. We start from a controlled illumination analysis and put forward a method of seismic survey geometry design for target-oriented imaging. The method needs a velocity model obtained by a preliminary seismic interpretation. The one-way Fourier finite-difference wave propagator is used to extrapolate plane wave sources on the target layer to the surface. By analyzing the wave energy distribution at the surface extrapolated from the target layer, the shot or receiver locations needed for target layer imaging can be determined. Numerical tests using the SEG-EAGE salt model suggest that this method is useful for confirming the special seismic acquisition geometry layout for target-oriented imaging.展开更多
This paper discusses the principle and procedures of the second-generation wavelet transform and its application to the denoising of seismic data. Based on lifting steps, it is a flexible wavelet construction method u...This paper discusses the principle and procedures of the second-generation wavelet transform and its application to the denoising of seismic data. Based on lifting steps, it is a flexible wavelet construction method using linear and nonlinear spatial prediction and operators to implement the wavelet transform and to make it reversible. The lifting scheme transform -includes three steps: split, predict, and update. Deslauriers-Dubuc (4, 2) wavelet transforms are used to process both synthetic and real data in our second-generation wavelet transform. The processing results show that random noise is effectively suppressed and the signal to noise ratio improves remarkably. The lifting wavelet transform is an efficient algorithm.展开更多
Air-gun arrays are used in marine-seismic exploration. Far-field wavelets in subsurface media represent the stacking of single air-gun ideal wavelets. We derived single air-gun ideal wavelets using near-field wavelets...Air-gun arrays are used in marine-seismic exploration. Far-field wavelets in subsurface media represent the stacking of single air-gun ideal wavelets. We derived single air-gun ideal wavelets using near-field wavelets recorded from near-field geophones and then synthesized them into far-field wavelets. This is critical for processing wavelets in marine- seismic exploration. For this purpose, several algorithms are currently used to decompose and synthesize wavelets in the time domain. If the traveltime of single air-gun wavelets is not an integral multiple of the sampling interval, the complex and error-prone resampling of the seismic signals using the time-domain method is necessary. Based on the relation between the frequency-domain phase and the time-domain time delay, we propose a method that first transforms the real near-field wavelet to the frequency domain via Fourier transforms; then, it decomposes it and composes the wavelet spectrum in the frequency domain, and then back transforms it to the time domain. Thus, the resampling problem is avoided and single air-gun wavelets and far-field wavelets can be reliably derived. The effect of ghost reflections is also considered, while decomposing the wavelet and removing the ghost reflections. Modeling and real data processing were used to demonstrate the feasibility of the proposed method.展开更多
In multi-component seismic exploration, the horizontal and vertical components both contain P- and SV-waves. The P- and SV-wavefields in a seismic record can be separated by their horizontal and vertical displacements...In multi-component seismic exploration, the horizontal and vertical components both contain P- and SV-waves. The P- and SV-wavefields in a seismic record can be separated by their horizontal and vertical displacements when upgoing P- and SV-waves arrive at the sea floor. If the sea floor P wave velocity, S wave velocity, and density are known, the separation can be achieved in ther-p domain. The separated wavefields are then transformed to the time domain. A method of separating P- and SV-wavefields is presented in this paper and used to effectively separate P- and SV-wavefields in synthetic and real data. The application to real data shows that this method is feasible and effective. It also can be used for free surface data.展开更多
The Ordos basin is rich in oil and gas resources. The structures in the basin are generally rather gentle. Most of the oil and gas pools already found are belong to the types of lithologic, lithostratigraphic and pale...The Ordos basin is rich in oil and gas resources. The structures in the basin are generally rather gentle. Most of the oil and gas pools already found are belong to the types of lithologic, lithostratigraphic and paleo-geomorphic with thin reservoirs. The geology of oil and gas and the generation and distribution of oil and gas resources in the Ordos basin is analyzed in detail in the paper. It is also pointed out in the paper that what direction should be taken in oil and gas exploration in the basin, and that the important role of the seismic prospecting technologies in the discovery of large oil and gas fields in the basin could not be replaced.展开更多
3D inversion of borehole-surface electrical data for complex geo-electrical models is still a challenging problem in geophysical exploration. We have developed a program for 3D inversion to borehole-surface electrical...3D inversion of borehole-surface electrical data for complex geo-electrical models is still a challenging problem in geophysical exploration. We have developed a program for 3D inversion to borehole-surface electrical data based on the quasi-analytical approximation (QA) and re-weighted regularized conjugate gradient method (RRCG) algorithms using Visual Fortran 6.5. Application of the QA approximation to forward modeling and Frechet derivative computations speeds up the calculation dramatically. The trial calculation for synthetic data of theoretical model showed that the program is fast and highly precise.展开更多
This article provides the application of the high-order, staggered-grid, finite-difference scheme to model elastic wave propagation in 3-D isotropic media. Here, we use second-order, tempo- ral- and high-order spatial...This article provides the application of the high-order, staggered-grid, finite-difference scheme to model elastic wave propagation in 3-D isotropic media. Here, we use second-order, tempo- ral- and high-order spatial finite-difference formulations with a staggered grid for discretization of the 3-D elastic wave equations of motion. The set of absorbing boundary conditions based on paraxial approximations of 3-D elastic wave equations are applied to the numerical boundaries. The trial re- sults for the salt model show that the numerical dispersion is decreased to a minimum extent, the accuracy high and diffracted waves abundant. It also shows that this method can be used for modeling wave propagation in complex media with the lateral variation of velocity.展开更多
The overburden failure causes the changes of the constant electric field, and the resistivity is the main parameter in these changes. The experimental simulation about tbe response relation between the overburden fail...The overburden failure causes the changes of the constant electric field, and the resistivity is the main parameter in these changes. The experimental simulation about tbe response relation between the overburden failure and its electrical parameter changes is made by building the similar material physics model of mining. The experiment results are used to analyze and test the in-situ detection. The research indicates that the resistivity changes with the electric characteristic of the rock in cracked zone and caving zone caused by overburden failure, the response characteristics of resistivity vary with the failure degrees at different overburden failure zone and that they are corresponding. The resistivity method used in monitoring the overburden failure can determine the height and the affecting scopes of the cracked zone and caving zone. This can provide reliable techaological guarantee for mining design and safe production.展开更多
This paper will present a detailed analysis of the deformation mechanism and stability assessment of the slope through field investigations, numerical modeling and measurements. Field investigation indicated that thre...This paper will present a detailed analysis of the deformation mechanism and stability assessment of the slope through field investigations, numerical modeling and measurements. Field investigation indicated that three thin coal seams encountered large mined-out area at one side and free surface of hill slope at the other side, which lead to the caving of roof strata movement, ground movement and crown crack along the preferred orientations of joints. The three-dimensional numeri- cal modeling study on the case demonstrated that the plasticity failure occurred gradually along with the extension of mined-out area in depth. When the depth of mining reached the verge defined by the seismic prospecting method, a large mount of tension failure occurred on the crown of the slope. The factor of safety was 1.36 calculated by the shear strength reduction technique, which indicated the slope was in stable state. The measurement showed that the residual deformation occurred before 1998 and became stable subsequently, which indicated that the residual deformation almost finished and the slope is in stable state.展开更多
We carried out surveying on the shallow structure and faulted-stratum of the Haihe fault in Tianjin using acoustic surveying methods such as the single-channel seismic technique. The result shows that the method can o...We carried out surveying on the shallow structure and faulted-stratum of the Haihe fault in Tianjin using acoustic surveying methods such as the single-channel seismic technique. The result shows that the method can obtain satisfactory results in wide and deep river courses. It also shows that in the Tanggu area of Tianjin, the upper fault point of Haihe fault is about 30m beneath the river bed and the corresponding latest active period is Qp^3- Qh^1 , which is consistent with the former borehole survey result. In the offshore area of the Bohai Sea, Haihe fault shows as a NWW-NEE strike dense fault zone and its upper fault point is less than 30m beneath the seabed. It shows that the active characteristics of Haihe fault in the Bohai Sea correspond to the Tanggu area.展开更多
Borehole gravity has been used in mineral exploration recently with the advent of slim-hole gravimeters. It is logical to perform inversion to utilize the information in the newly acquired data. The inversions were ca...Borehole gravity has been used in mineral exploration recently with the advent of slim-hole gravimeters. It is logical to perform inversion to utilize the information in the newly acquired data. The inversions were carried out by using cokriging,which is a geostatistical method of estimation that minimizes the error variance by applying cross-correlation between several variables. In this study the estimated densities are derived by using boreholes gravity and known densities along the borehole. This method does not need iterative process and computes efficiently. The selection of examples demonstrates that this method has the ability to include physical property from borehole measurements in the inversion. The synthetic examples demonstrate the density variation along a borehole can be well determined without depth constraints in the inversion. The resolution of the recovered model can be further improved by including the densities along the borehole for inversion. However,this capability decreases dramatically with the increasing of distance between the anomalous body and the borehole.展开更多
By shallow seismic prospecting, the Cenozoic Group in the sea area near the Yangtze River Mouth can be divided into five seismic sequences. They correspond to the Quaternary, Pliocene, Upper Miocene, Lower Miocene and...By shallow seismic prospecting, the Cenozoic Group in the sea area near the Yangtze River Mouth can be divided into five seismic sequences. They correspond to the Quaternary, Pliocene, Upper Miocene, Lower Miocene and Eocene respectively. The Quaternary System covers all the detecting area. The Tertiary System overlaps and thins out from NE to SW. The sedimentary basement mainly consists of volcanic rock (J 3) and acidic rock (r 3 5). Paleogene or Late Cretaceous basins are not found there. The faults that have been detected are all normal faults. They can be divided into three groups (NE, NW, near EW) by their trend. The NE and NW trending faults are predominant, and agree with aeromagnetic anomaly. Their length and displacement are larger than that of the EW trending faults. The activity of the NE trending faults is different in different segments. The SW segment is a Quaternary fault, the middle segment is a Neogene fault, The NE is Paleogene. But the segment of the NW trending fault is not obvious. The average vertical displacement rate is about 0 015mm/a.展开更多
基金supported by National Natural Science Foundation of China (Grant No. 40674072)National Hi-techResearch and Development Program of China (863 Program) (Grant No. 2006AA09A102-08)+1 种基金National Basic ResearchProgram of China (the 973 Program. Grant No. 007CB209603) the Opening fund of State Key Laboratory of Geological Processes and Mineral Resources, China University of Geosciences (GPMR200633)
文摘Based on a detailed analysis of differences between seismic data and well logs, we discuss the problem of matching seismic traces and well logs and present a new matching method based on event search in instantaneous phase which greatly improves seismic resolution. The method is based on flattening events in instantaneous phase to compare the seismic traces to the well log traces with the same phase. We calculate the coefficients using the singular value decomposition method to extrapolate the well logs. As a result, the events in the seismic profile are continuous and match well with well logs. We apply this method to the Mao-2 well in Daqing Oilfield with good results.
基金supported by Sino Probe Project(Sino Probe-09-02,201011079)Development Project of National Key Scientific Equipment(No.ZDYZ2012-1-05)the Strategic Priority Research Program of Chinese Academy of Sciences(No.XDA14050100)
文摘To test the performance of the Chinese whole-surface electromagnetic prospecting (SEP) system, system integrations, instrument performances, and large-scale production viabilities in Liaoning province and Inner Mongolia were measured via extensive field tests. Resultant electric fields, magnetic fields, apparent resistivities, impedance phases, and inversion profiles compared favorably with results of commercial equipment from other countries. The inversion results agreed well with the geologic information from boreholes. Field tests showed that the SEP system is stable, reliable, lightweight, and easy to operate, making it suitable and ready for real-field exploration.
文摘The conventional method of seismic data acquisition geometry design is based on the assumption of horizontal subsurface reflectors, which often is not suitable for complex structure. We start from a controlled illumination analysis and put forward a method of seismic survey geometry design for target-oriented imaging. The method needs a velocity model obtained by a preliminary seismic interpretation. The one-way Fourier finite-difference wave propagator is used to extrapolate plane wave sources on the target layer to the surface. By analyzing the wave energy distribution at the surface extrapolated from the target layer, the shot or receiver locations needed for target layer imaging can be determined. Numerical tests using the SEG-EAGE salt model suggest that this method is useful for confirming the special seismic acquisition geometry layout for target-oriented imaging.
文摘This paper discusses the principle and procedures of the second-generation wavelet transform and its application to the denoising of seismic data. Based on lifting steps, it is a flexible wavelet construction method using linear and nonlinear spatial prediction and operators to implement the wavelet transform and to make it reversible. The lifting scheme transform -includes three steps: split, predict, and update. Deslauriers-Dubuc (4, 2) wavelet transforms are used to process both synthetic and real data in our second-generation wavelet transform. The processing results show that random noise is effectively suppressed and the signal to noise ratio improves remarkably. The lifting wavelet transform is an efficient algorithm.
基金supported by the Geosciences and Technology Academy of China University of Petroleum(East China)
文摘Air-gun arrays are used in marine-seismic exploration. Far-field wavelets in subsurface media represent the stacking of single air-gun ideal wavelets. We derived single air-gun ideal wavelets using near-field wavelets recorded from near-field geophones and then synthesized them into far-field wavelets. This is critical for processing wavelets in marine- seismic exploration. For this purpose, several algorithms are currently used to decompose and synthesize wavelets in the time domain. If the traveltime of single air-gun wavelets is not an integral multiple of the sampling interval, the complex and error-prone resampling of the seismic signals using the time-domain method is necessary. Based on the relation between the frequency-domain phase and the time-domain time delay, we propose a method that first transforms the real near-field wavelet to the frequency domain via Fourier transforms; then, it decomposes it and composes the wavelet spectrum in the frequency domain, and then back transforms it to the time domain. Thus, the resampling problem is avoided and single air-gun wavelets and far-field wavelets can be reliably derived. The effect of ghost reflections is also considered, while decomposing the wavelet and removing the ghost reflections. Modeling and real data processing were used to demonstrate the feasibility of the proposed method.
基金This research is sponsored by National Natural Science Foundation of China (No. 40272041) and Innovative Foundation of CNPC (N0. 04E702).
文摘In multi-component seismic exploration, the horizontal and vertical components both contain P- and SV-waves. The P- and SV-wavefields in a seismic record can be separated by their horizontal and vertical displacements when upgoing P- and SV-waves arrive at the sea floor. If the sea floor P wave velocity, S wave velocity, and density are known, the separation can be achieved in ther-p domain. The separated wavefields are then transformed to the time domain. A method of separating P- and SV-wavefields is presented in this paper and used to effectively separate P- and SV-wavefields in synthetic and real data. The application to real data shows that this method is feasible and effective. It also can be used for free surface data.
文摘The Ordos basin is rich in oil and gas resources. The structures in the basin are generally rather gentle. Most of the oil and gas pools already found are belong to the types of lithologic, lithostratigraphic and paleo-geomorphic with thin reservoirs. The geology of oil and gas and the generation and distribution of oil and gas resources in the Ordos basin is analyzed in detail in the paper. It is also pointed out in the paper that what direction should be taken in oil and gas exploration in the basin, and that the important role of the seismic prospecting technologies in the discovery of large oil and gas fields in the basin could not be replaced.
文摘3D inversion of borehole-surface electrical data for complex geo-electrical models is still a challenging problem in geophysical exploration. We have developed a program for 3D inversion to borehole-surface electrical data based on the quasi-analytical approximation (QA) and re-weighted regularized conjugate gradient method (RRCG) algorithms using Visual Fortran 6.5. Application of the QA approximation to forward modeling and Frechet derivative computations speeds up the calculation dramatically. The trial calculation for synthetic data of theoretical model showed that the program is fast and highly precise.
文摘This article provides the application of the high-order, staggered-grid, finite-difference scheme to model elastic wave propagation in 3-D isotropic media. Here, we use second-order, tempo- ral- and high-order spatial finite-difference formulations with a staggered grid for discretization of the 3-D elastic wave equations of motion. The set of absorbing boundary conditions based on paraxial approximations of 3-D elastic wave equations are applied to the numerical boundaries. The trial re- sults for the salt model show that the numerical dispersion is decreased to a minimum extent, the accuracy high and diffracted waves abundant. It also shows that this method can be used for modeling wave propagation in complex media with the lateral variation of velocity.
文摘The overburden failure causes the changes of the constant electric field, and the resistivity is the main parameter in these changes. The experimental simulation about tbe response relation between the overburden failure and its electrical parameter changes is made by building the similar material physics model of mining. The experiment results are used to analyze and test the in-situ detection. The research indicates that the resistivity changes with the electric characteristic of the rock in cracked zone and caving zone caused by overburden failure, the response characteristics of resistivity vary with the failure degrees at different overburden failure zone and that they are corresponding. The resistivity method used in monitoring the overburden failure can determine the height and the affecting scopes of the cracked zone and caving zone. This can provide reliable techaological guarantee for mining design and safe production.
基金Supported by the National Natural Science Foundation of China for Youth (51004065) the National Basic Research Program of China (2012CB724208)
文摘This paper will present a detailed analysis of the deformation mechanism and stability assessment of the slope through field investigations, numerical modeling and measurements. Field investigation indicated that three thin coal seams encountered large mined-out area at one side and free surface of hill slope at the other side, which lead to the caving of roof strata movement, ground movement and crown crack along the preferred orientations of joints. The three-dimensional numeri- cal modeling study on the case demonstrated that the plasticity failure occurred gradually along with the extension of mined-out area in depth. When the depth of mining reached the verge defined by the seismic prospecting method, a large mount of tension failure occurred on the crown of the slope. The factor of safety was 1.36 calculated by the shear strength reduction technique, which indicated the slope was in stable state. The measurement showed that the residual deformation occurred before 1998 and became stable subsequently, which indicated that the residual deformation almost finished and the slope is in stable state.
基金the China Digital Seismic Observation Network Program and the Active Fault Detecting Sub-program of the tenth"Five-Year Plan" of China Earthquake Administration
文摘We carried out surveying on the shallow structure and faulted-stratum of the Haihe fault in Tianjin using acoustic surveying methods such as the single-channel seismic technique. The result shows that the method can obtain satisfactory results in wide and deep river courses. It also shows that in the Tanggu area of Tianjin, the upper fault point of Haihe fault is about 30m beneath the river bed and the corresponding latest active period is Qp^3- Qh^1 , which is consistent with the former borehole survey result. In the offshore area of the Bohai Sea, Haihe fault shows as a NWW-NEE strike dense fault zone and its upper fault point is less than 30m beneath the seabed. It shows that the active characteristics of Haihe fault in the Bohai Sea correspond to the Tanggu area.
基金Supported by the National High Technology Research and Development Program(863 Program)(No.2014AA06A613)by Project of Graduate Innovation Fund of Jilin University(No.2014066)
文摘Borehole gravity has been used in mineral exploration recently with the advent of slim-hole gravimeters. It is logical to perform inversion to utilize the information in the newly acquired data. The inversions were carried out by using cokriging,which is a geostatistical method of estimation that minimizes the error variance by applying cross-correlation between several variables. In this study the estimated densities are derived by using boreholes gravity and known densities along the borehole. This method does not need iterative process and computes efficiently. The selection of examples demonstrates that this method has the ability to include physical property from borehole measurements in the inversion. The synthetic examples demonstrate the density variation along a borehole can be well determined without depth constraints in the inversion. The resolution of the recovered model can be further improved by including the densities along the borehole for inversion. However,this capability decreases dramatically with the increasing of distance between the anomalous body and the borehole.
文摘By shallow seismic prospecting, the Cenozoic Group in the sea area near the Yangtze River Mouth can be divided into five seismic sequences. They correspond to the Quaternary, Pliocene, Upper Miocene, Lower Miocene and Eocene respectively. The Quaternary System covers all the detecting area. The Tertiary System overlaps and thins out from NE to SW. The sedimentary basement mainly consists of volcanic rock (J 3) and acidic rock (r 3 5). Paleogene or Late Cretaceous basins are not found there. The faults that have been detected are all normal faults. They can be divided into three groups (NE, NW, near EW) by their trend. The NE and NW trending faults are predominant, and agree with aeromagnetic anomaly. Their length and displacement are larger than that of the EW trending faults. The activity of the NE trending faults is different in different segments. The SW segment is a Quaternary fault, the middle segment is a Neogene fault, The NE is Paleogene. But the segment of the NW trending fault is not obvious. The average vertical displacement rate is about 0 015mm/a.