The nonsplitting perfectly matched layer (NPML) absorbing boundary condition (ABC) was first provided by Wang and Tang (2003) for the finite-difference simulation of elastic wave propagation in solids. In this p...The nonsplitting perfectly matched layer (NPML) absorbing boundary condition (ABC) was first provided by Wang and Tang (2003) for the finite-difference simulation of elastic wave propagation in solids. In this paper, the method is developed to extend the NPML to simulating elastic wave propagation in poroelastic media. Biot's equations are discretized and approximated to a staggered-grid by applying a fourth-order accurate central difference in space and a second-order accurate central difference in time. A cylindrical twolayer seismic model and a borehole model are chosen to validate the effectiveness of the NPML. The results show that the numerical solutions agree well with the solutions of the discrete wavenumber (DW) method.展开更多
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.展开更多
Electromagnetic holographic data are characterized by two modes, suggesting that image reconstruction requires a dual-mode sensitivity field as well. We analyze an electromagnetic holographic field based on tomography...Electromagnetic holographic data are characterized by two modes, suggesting that image reconstruction requires a dual-mode sensitivity field as well. We analyze an electromagnetic holographic field based on tomography theory and Radon inverse transform to derive the expression of the electromagnetic holographic sensitivity field (EMHSF). Then, we apply the EMHSF calculated by using finite-element methods to flow simulations and holographic imaging. The results suggest that the EMHSF based on the partial derivative of radius of the complex electric potential φ is closely linked to the Radon inverse transform and encompasses the sensitivities of the amplitude and phase data. The flow images obtained with inversion using EMHSF better agree with the actual flow patterns. The EMHSF overcomes the limitations of traditional single-mode sensitivity fields.展开更多
Ground-penetrating radar (GPR) has been used predominantly for environments with low electrical conductivity like freshwater aquifers, glaciers, or dry sandy soils. The objective of the present study was to explore it...Ground-penetrating radar (GPR) has been used predominantly for environments with low electrical conductivity like freshwater aquifers, glaciers, or dry sandy soils. The objective of the present study was to explore its application for mapping in subsurface agricultural soils to a depth of several meters. For a loamy sand and a clayey site on the North China Plain, clay inclusions in the sand were detected; the thickness, inclination, and continuity of the confining clay and silt layers was assessed; and a local water table was mapped. Direct sampling (soil coring and profiling) in the top meter and independent measurement of the water table were utilized to confirm the findings. Also, effective estimates of the dielectric number for the site with the dielectric number of moist clayey soils depending strongly on frequency were obtained. Thus, important properties of soils, like the arrangement and type of layers and in particular their continuity and inclination, could be explored with moderate efforts for rather large areas to help find optimal locations for the time-consuming and expensive measurements which would be necessary to detail a model of the subsurface.展开更多
On the standpoint of the disaster prevention from water inrush,discussed the genesis and geologic condition of karstic collapse column in one coal field,analyzed the geophysical characteristics of karstic collapse col...On the standpoint of the disaster prevention from water inrush,discussed the genesis and geologic condition of karstic collapse column in one coal field,analyzed the geophysical characteristics of karstic collapse column by using high resolution 3D seismic data.It shows the effective result of the technology of high resolution 3D seismic pros- pecting in the exploration of the karstic collapse column,and presents some prediction methods and prevention measures.展开更多
Field SZ36-1 is a water-flooded heavy oil reservoir with high porosity and unconsolidated sand. The recovery rate is low so that it becomes a challenge for production. Time-lapse seismic data is studied to improve the...Field SZ36-1 is a water-flooded heavy oil reservoir with high porosity and unconsolidated sand. The recovery rate is low so that it becomes a challenge for production. Time-lapse seismic data is studied to improve the oil recovery for this field, This feasibility study analyzes the possible time-lapse seismic attribute spatial distribution using dynamic data and the reservoir model to determine the optimum time to acquire a new seismic survey. Based on the study, it is found that the time-lapse seismic response for this unconsolidated sand has a strong signature due to solution gas when the reservoir pressure is below the bubble point. This indicates that acquiring a new survey after 10 years of production is appropriate for a time-lapse seismic application.展开更多
In order to detect fault exactly and quickly, cusp catastrophe theory is used to interpret 3D coal seismic data in this paper. By establishing a cusp model, seismic signal is transformed into standard form of cusp cat...In order to detect fault exactly and quickly, cusp catastrophe theory is used to interpret 3D coal seismic data in this paper. By establishing a cusp model, seismic signal is transformed into standard form of cusp catastrophe and catastrophe parameters, including time-domain catastrophe potential, time-domain catastrophe time, frequency-domain catastrophe potential and frequency- domain degree, are calculated. Catastrophe theory is used in 3D seismic structural interpretation in coal mine. The results show that the position of abnormality of the catastrophe parameter profile or curve is related to the location of fault, and the cusp catastrophe theory is effective to automatically pick up geology information and improve the interpretation precision in 3D seismic data.展开更多
The Connolly (1999) elastic impedance (EI) equation is a function of P-wave velocity, S-wave velocity, density, and incidence angle. Conventional inversion methods based on this equation can only extract P-velocit...The Connolly (1999) elastic impedance (EI) equation is a function of P-wave velocity, S-wave velocity, density, and incidence angle. Conventional inversion methods based on this equation can only extract P-velocity, S-velocity, and density data directly and the elastic impedance at different incidence angles are not at the same scale, which makes comparison difficult. We propose a new elastic impedance equation based on the Gray et al. (1999) Zoeppritz approximation using Lamé parameters to address the conventional inversion method's deficiencies. This equation has been normalized to unify the elastic impedance dimensions at different angles and used for inversion. Lamé parameters can be extracted directly from the elastic impedance data obtained from inversion using the linear relation between Lamé parameters and elastic impedance. The application example shows that the elastic parameters extracted using this new method are more stable and correct and can recover the reservoir information very well. The new method is an improvement on the conventional method based on Connolly's equation.展开更多
The Formation Evaluation Tool (FET) introduced in the paper represents a new generation of formation evaluation systems developed and manufactured by China Oilfield Services Limited (COSL), CNOOC, using a FET tech...The Formation Evaluation Tool (FET) introduced in the paper represents a new generation of formation evaluation systems developed and manufactured by China Oilfield Services Limited (COSL), CNOOC, using a FET technology transfer from Crocker Research, Australia. The system has been applied successfully in the Bohai Sea and South China Sea. For instance, a multilayered oil and water system has been confirmed with the aid of accurate formation pressure tests, even in very thin beds and edge water reservoirs, overcoming the difficulty of determining this kind of oil-water and gas-water contacts. Moreover, the FET pumping and real-time fluid monitoring function allows acquiring a true sample of formation fluid unpolluted by drilling mud which plays an important role in determining the fluid properties of the target stratum and analyzing the fluid component. The principles and purpose of the Formation Evaluation Tool (FET) will be briefly introduced and successful examples of the application of the technology will be described in detail in this paper.展开更多
In order to increase the exploration depth of Rayleigh wave, new idea that dif-ferent from the former principles in data acquisition was applied. Suitable data acquisition parameter was given out on the basis of large...In order to increase the exploration depth of Rayleigh wave, new idea that dif-ferent from the former principles in data acquisition was applied. Suitable data acquisition parameter was given out on the basis of large amount of experiments. By reducing the group interval, the low frequency signal are enhanced instead of been attenuated. Fur-thermore, to solve the problem that the precision of Rayleigh wave exploration method count much to the signal-to-noise ratio, some preprocessing methods were put forward. By using zero shift rectifying, digital F-K filtering and cutting, noises can be effectively eliminated.展开更多
The depth of upper fault point is the key data for ascertaining the active age of a buried fault on a plain. The difference of depth obtained from same fault may be dozens to several hundred meters when using differen...The depth of upper fault point is the key data for ascertaining the active age of a buried fault on a plain. The difference of depth obtained from same fault may be dozens to several hundred meters when using different geophysical methods. It can result in the absolutely opposite conclusions when judging fault activity. Because of a lack of an artificial earthquake source with wide band and high central frequency, many kinds of methods have to be used together. The higher the frequency of the artificial earthquake wave, electromagnetic wave and sonic wave, the higher the resolution. However the attenuation is also very fast and the exploration depth is very shallow. The reverse is also true. The frequency of artificial seismic waves is in the tens of Hz. Its exploration depth is big and the resolution is poor. The frequency of radar electromagnetic waves is about a million Hz, indicating that the resolving power is better, but the exploration depth is very shallow. However, the acoustic frequency is thousands of Hz, its resolving power is better than that of the artificial earthquake method and the exploration depth is larger than that of the radar method. So it is suitable for extra shallow exploration in the thick deposit strata of the Quaternary. The preliminary results detected using the high frequency acoustic method in extra shallow layers indicates that previous inferences about some fault activity in the eastern part of the North China plain may need to be greatly corrected.展开更多
Branching river channels and the coexistence of valleys, ridges, hiils, and slopes'as the result of long-term weathering and erosion form the unique loess topography. The Changqing Geophysical Company, working in the...Branching river channels and the coexistence of valleys, ridges, hiils, and slopes'as the result of long-term weathering and erosion form the unique loess topography. The Changqing Geophysical Company, working in these complex conditions, has established a suite of technologies for high-fidelity processing and fine interpretation of seismic data. This article introduces the processes involved in the data processing and interpretation and illustrates the results.展开更多
All coal mine disasters are dynamic geological phenomenon and affected by many factors. However, locating the enriched areas of CSM (coal seam methane) may be the precondition for the successful prediction of such dis...All coal mine disasters are dynamic geological phenomenon and affected by many factors. However, locating the enriched areas of CSM (coal seam methane) may be the precondition for the successful prediction of such disasters. Traditional methods of investigating CSM enriched areas use limited data and only consider a few important factors. Their success rate is low and cannot meet practical needs. In this paper, an alternative method is proposed. The proce- dure is given as follows: 1) fracture attributes derived from azimuth variations of P-wave data in coal seams and wall rocks can be extracted; 2) AVO attributes, such as the intercept P and gradient G parameters can be extracted from different azimuths from 3D seismic data; 3) seismic cubes can be inverted and the relative attributes of imped- ance cubes can be extracted; 4) using a GIS platform, multi-source information can be obtained and analyzed; these include fracture attributes of coal seams and wall rocks, the thickness of coal seams, the distribution of faults and structures, the depth of coal seams, the inclination and exposure of coal seams and the coal rank. Through this processing procedure, methane enriched areas can be systematically detected.展开更多
Based on the biological resources surveys in China’s exclusive economic zone, the paper discusses the application of the GIS technology in evaluating the marine environmental quality, establishes the basic concept an...Based on the biological resources surveys in China’s exclusive economic zone, the paper discusses the application of the GIS technology in evaluating the marine environmental quality, establishes the basic concept and technological process for evaluating the environmental quality of the habitats of marine living resources, proposes the key points for the establishment of the relevant evaluation attributes database by applying the GIS technology and analyses several typical evaluation models and the data interpolation method for the model in GIS. Results show that the combination of the GIS technology with the quantitative evaluation model can be applied to the fast and accurate evaluation of the marine environmental quality.展开更多
A convenient numerical calculation method (inverse spline interpolation) for all-time apparent resistivity intransient electromagnetic method (TEM) is proposed in this paper. Characteristic of early and late normalize...A convenient numerical calculation method (inverse spline interpolation) for all-time apparent resistivity intransient electromagnetic method (TEM) is proposed in this paper. Characteristic of early and late normalized inductiveelectromotive force was investigated. According to the turning point, the transient process is divided into the earlyphase, the turning point, and the late phase. Afterwards, apparent resistivity is obtained through inverse spline interpo-lation in the early and the late phases, respectively. Finally, the resistivities of the early-time and the late-time wereconnected together by the turning point. The result shows that the inverse spline method is feasible and the method alsolays a foundation for initial model construction in the TEM automatic inversion.展开更多
The present paper investigates the application of high resolution magnetic survey to detecting igneous bodies. The slight difference in magnetism between igneous bodies and their surrounding rocks is measured first an...The present paper investigates the application of high resolution magnetic survey to detecting igneous bodies. The slight difference in magnetism between igneous bodies and their surrounding rocks is measured first and then the magnetic survey data are processed to determine whether there exist igneous bodies by analog among several measuring lines, and finally the modified Marquart inversion was used to determine the occurrence and distribution of the igneous bodies.展开更多
基金This research was supported by Natural Science Foundation of China (No. 403740043).
文摘The nonsplitting perfectly matched layer (NPML) absorbing boundary condition (ABC) was first provided by Wang and Tang (2003) for the finite-difference simulation of elastic wave propagation in solids. In this paper, the method is developed to extend the NPML to simulating elastic wave propagation in poroelastic media. Biot's equations are discretized and approximated to a staggered-grid by applying a fourth-order accurate central difference in space and a second-order accurate central difference in time. A cylindrical twolayer seismic model and a borehole model are chosen to validate the effectiveness of the NPML. The results show that the numerical solutions agree well with the solutions of the discrete wavenumber (DW) method.
基金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 the National Science and Technology Major Project(No.2011ZX05020-006)
文摘Electromagnetic holographic data are characterized by two modes, suggesting that image reconstruction requires a dual-mode sensitivity field as well. We analyze an electromagnetic holographic field based on tomography theory and Radon inverse transform to derive the expression of the electromagnetic holographic sensitivity field (EMHSF). Then, we apply the EMHSF calculated by using finite-element methods to flow simulations and holographic imaging. The results suggest that the EMHSF based on the partial derivative of radius of the complex electric potential φ is closely linked to the Radon inverse transform and encompasses the sensitivities of the amplitude and phase data. The flow images obtained with inversion using EMHSF better agree with the actual flow patterns. The EMHSF overcomes the limitations of traditional single-mode sensitivity fields.
基金Project supported in part by the Deutsche Forschungsgemeinschaft (DFG), Germany (No. RO 1080/8-1) jointly by Max-Planck Gesellschaft and the Chinese Academy of Sciences through a travel grant to the first author.
文摘Ground-penetrating radar (GPR) has been used predominantly for environments with low electrical conductivity like freshwater aquifers, glaciers, or dry sandy soils. The objective of the present study was to explore its application for mapping in subsurface agricultural soils to a depth of several meters. For a loamy sand and a clayey site on the North China Plain, clay inclusions in the sand were detected; the thickness, inclination, and continuity of the confining clay and silt layers was assessed; and a local water table was mapped. Direct sampling (soil coring and profiling) in the top meter and independent measurement of the water table were utilized to confirm the findings. Also, effective estimates of the dielectric number for the site with the dielectric number of moist clayey soils depending strongly on frequency were obtained. Thus, important properties of soils, like the arrangement and type of layers and in particular their continuity and inclination, could be explored with moderate efforts for rather large areas to help find optimal locations for the time-consuming and expensive measurements which would be necessary to detail a model of the subsurface.
基金the National Natural Science Foundation of China(2007CB209600)
文摘On the standpoint of the disaster prevention from water inrush,discussed the genesis and geologic condition of karstic collapse column in one coal field,analyzed the geophysical characteristics of karstic collapse column by using high resolution 3D seismic data.It shows the effective result of the technology of high resolution 3D seismic pros- pecting in the exploration of the karstic collapse column,and presents some prediction methods and prevention measures.
文摘Field SZ36-1 is a water-flooded heavy oil reservoir with high porosity and unconsolidated sand. The recovery rate is low so that it becomes a challenge for production. Time-lapse seismic data is studied to improve the oil recovery for this field, This feasibility study analyzes the possible time-lapse seismic attribute spatial distribution using dynamic data and the reservoir model to determine the optimum time to acquire a new seismic survey. Based on the study, it is found that the time-lapse seismic response for this unconsolidated sand has a strong signature due to solution gas when the reservoir pressure is below the bubble point. This indicates that acquiring a new survey after 10 years of production is appropriate for a time-lapse seismic application.
文摘In order to detect fault exactly and quickly, cusp catastrophe theory is used to interpret 3D coal seismic data in this paper. By establishing a cusp model, seismic signal is transformed into standard form of cusp catastrophe and catastrophe parameters, including time-domain catastrophe potential, time-domain catastrophe time, frequency-domain catastrophe potential and frequency- domain degree, are calculated. Catastrophe theory is used in 3D seismic structural interpretation in coal mine. The results show that the position of abnormality of the catastrophe parameter profile or curve is related to the location of fault, and the cusp catastrophe theory is effective to automatically pick up geology information and improve the interpretation precision in 3D seismic data.
文摘The Connolly (1999) elastic impedance (EI) equation is a function of P-wave velocity, S-wave velocity, density, and incidence angle. Conventional inversion methods based on this equation can only extract P-velocity, S-velocity, and density data directly and the elastic impedance at different incidence angles are not at the same scale, which makes comparison difficult. We propose a new elastic impedance equation based on the Gray et al. (1999) Zoeppritz approximation using Lamé parameters to address the conventional inversion method's deficiencies. This equation has been normalized to unify the elastic impedance dimensions at different angles and used for inversion. Lamé parameters can be extracted directly from the elastic impedance data obtained from inversion using the linear relation between Lamé parameters and elastic impedance. The application example shows that the elastic parameters extracted using this new method are more stable and correct and can recover the reservoir information very well. The new method is an improvement on the conventional method based on Connolly's equation.
文摘The Formation Evaluation Tool (FET) introduced in the paper represents a new generation of formation evaluation systems developed and manufactured by China Oilfield Services Limited (COSL), CNOOC, using a FET technology transfer from Crocker Research, Australia. The system has been applied successfully in the Bohai Sea and South China Sea. For instance, a multilayered oil and water system has been confirmed with the aid of accurate formation pressure tests, even in very thin beds and edge water reservoirs, overcoming the difficulty of determining this kind of oil-water and gas-water contacts. Moreover, the FET pumping and real-time fluid monitoring function allows acquiring a true sample of formation fluid unpolluted by drilling mud which plays an important role in determining the fluid properties of the target stratum and analyzing the fluid component. The principles and purpose of the Formation Evaluation Tool (FET) will be briefly introduced and successful examples of the application of the technology will be described in detail in this paper.
文摘In order to increase the exploration depth of Rayleigh wave, new idea that dif-ferent from the former principles in data acquisition was applied. Suitable data acquisition parameter was given out on the basis of large amount of experiments. By reducing the group interval, the low frequency signal are enhanced instead of been attenuated. Fur-thermore, to solve the problem that the precision of Rayleigh wave exploration method count much to the signal-to-noise ratio, some preprocessing methods were put forward. By using zero shift rectifying, digital F-K filtering and cutting, noises can be effectively eliminated.
文摘The depth of upper fault point is the key data for ascertaining the active age of a buried fault on a plain. The difference of depth obtained from same fault may be dozens to several hundred meters when using different geophysical methods. It can result in the absolutely opposite conclusions when judging fault activity. Because of a lack of an artificial earthquake source with wide band and high central frequency, many kinds of methods have to be used together. The higher the frequency of the artificial earthquake wave, electromagnetic wave and sonic wave, the higher the resolution. However the attenuation is also very fast and the exploration depth is very shallow. The reverse is also true. The frequency of artificial seismic waves is in the tens of Hz. Its exploration depth is big and the resolution is poor. The frequency of radar electromagnetic waves is about a million Hz, indicating that the resolving power is better, but the exploration depth is very shallow. However, the acoustic frequency is thousands of Hz, its resolving power is better than that of the artificial earthquake method and the exploration depth is larger than that of the radar method. So it is suitable for extra shallow exploration in the thick deposit strata of the Quaternary. The preliminary results detected using the high frequency acoustic method in extra shallow layers indicates that previous inferences about some fault activity in the eastern part of the North China plain may need to be greatly corrected.
文摘Branching river channels and the coexistence of valleys, ridges, hiils, and slopes'as the result of long-term weathering and erosion form the unique loess topography. The Changqing Geophysical Company, working in these complex conditions, has established a suite of technologies for high-fidelity processing and fine interpretation of seismic data. This article introduces the processes involved in the data processing and interpretation and illustrates the results.
基金Project 40574057 supported by the National Natural Science Foundation of China and CUMT Youth Foundation
文摘All coal mine disasters are dynamic geological phenomenon and affected by many factors. However, locating the enriched areas of CSM (coal seam methane) may be the precondition for the successful prediction of such disasters. Traditional methods of investigating CSM enriched areas use limited data and only consider a few important factors. Their success rate is low and cannot meet practical needs. In this paper, an alternative method is proposed. The proce- dure is given as follows: 1) fracture attributes derived from azimuth variations of P-wave data in coal seams and wall rocks can be extracted; 2) AVO attributes, such as the intercept P and gradient G parameters can be extracted from different azimuths from 3D seismic data; 3) seismic cubes can be inverted and the relative attributes of imped- ance cubes can be extracted; 4) using a GIS platform, multi-source information can be obtained and analyzed; these include fracture attributes of coal seams and wall rocks, the thickness of coal seams, the distribution of faults and structures, the depth of coal seams, the inclination and exposure of coal seams and the coal rank. Through this processing procedure, methane enriched areas can be systematically detected.
文摘Based on the biological resources surveys in China’s exclusive economic zone, the paper discusses the application of the GIS technology in evaluating the marine environmental quality, establishes the basic concept and technological process for evaluating the environmental quality of the habitats of marine living resources, proposes the key points for the establishment of the relevant evaluation attributes database by applying the GIS technology and analyses several typical evaluation models and the data interpolation method for the model in GIS. Results show that the combination of the GIS technology with the quantitative evaluation model can be applied to the fast and accurate evaluation of the marine environmental quality.
基金Project 40344022 supported by National Natural Science Foundation of China
文摘A convenient numerical calculation method (inverse spline interpolation) for all-time apparent resistivity intransient electromagnetic method (TEM) is proposed in this paper. Characteristic of early and late normalized inductiveelectromotive force was investigated. According to the turning point, the transient process is divided into the earlyphase, the turning point, and the late phase. Afterwards, apparent resistivity is obtained through inverse spline interpo-lation in the early and the late phases, respectively. Finally, the resistivities of the early-time and the late-time wereconnected together by the turning point. The result shows that the inverse spline method is feasible and the method alsolays a foundation for initial model construction in the TEM automatic inversion.
文摘The present paper investigates the application of high resolution magnetic survey to detecting igneous bodies. The slight difference in magnetism between igneous bodies and their surrounding rocks is measured first and then the magnetic survey data are processed to determine whether there exist igneous bodies by analog among several measuring lines, and finally the modified Marquart inversion was used to determine the occurrence and distribution of the igneous bodies.