Seismic impedance inversion is an important technique for structure identification and reservoir prediction.Model-based and data-driven impedance inversion are the commonly used inversion methods.In practice,the geoph...Seismic impedance inversion is an important technique for structure identification and reservoir prediction.Model-based and data-driven impedance inversion are the commonly used inversion methods.In practice,the geophysical inversion problem is essentially an ill-posedness problem,which means that there are many solutions corresponding to the same seismic data.Therefore,regularization schemes,which can provide stable and unique inversion results to some extent,have been introduced into the objective function as constrain terms.Among them,given a low-frequency initial impedance model is the most commonly used regularization method,which can provide a smooth and stable solution.However,this model-based inversion method relies heavily on the initial model and the inversion result is band limited to the effective frequency bandwidth of seismic data,which cannot effectively improve the seismic vertical resolution and is difficult to be applied to complex structural regions.Therefore,we propose a data-driven approach for high-resolution impedance inversion based on the bidirectional long short-term memory recurrent neural network,which regards seismic data as time-series rather than image-like patches.Compared with the model-based inversion method,the data-driven approach provides higher resolution inversion results,which demonstrates the effectiveness of the data-driven method for recovering the high-frequency components.However,judging from the inversion results for characterization the spatial distribution of thin-layer sands,the accuracy of high-frequency components is difficult to guarantee.Therefore,we add the model constraint to the objective function to overcome the shortages of relying only on the data-driven schemes.First,constructing the supervisor1 based on the bidirectional long short-term memory recurrent neural network,which provides the predicted impedance with higher resolution.Then,convolution constraint as supervisor2 is introduced into the objective function to guarantee the reliability and accuracy of the inversion results,which makes the synthetic seismic data obtained from the inversion result consistent with the input data.Finally,we test the proposed scheme based on the synthetic and field seismic data.Compared to model-based and purely data-driven impedance inversion methods,the proposed approach provides more accurate and reliable inversion results while with higher vertical resolution and better spatial continuity.The inversion results accurately characterize the spatial distribution relationship of thin sands.The model tests demonstrate that the model-constrained and data-driven impedance inversion scheme can effectively improve the thin-layer structure characterization based on the seismic data.Moreover,tests on the oil field data indicate the practicality and adaptability of the proposed method.展开更多
The chirp sub-bottom profiler,for its high resolution,easy accessibility and cost-effectiveness,has been widely used in acoustic detection.In this paper,the acoustic impedance and grain size compositions were obtained...The chirp sub-bottom profiler,for its high resolution,easy accessibility and cost-effectiveness,has been widely used in acoustic detection.In this paper,the acoustic impedance and grain size compositions were obtained based on the chirp sub-bottom profiler data collected in the Chukchi Plateau area during the 11th Arctic Expedition of China.The time-domain adaptive search matching algorithm was used and validated on our established theoretical model.The misfit between the inversion result and the theoretical model is less than 0.067%.The grain size was calculated according to the empirical relationship between the acoustic impedance and the grain size of the sediment.The average acoustic impedance of sub-seafloor strata is 2.5026×10^(6) kg(s m^(2))^(-1)and the average grain size(θvalue)of the seafloor surface sediment is 7.1498,indicating the predominant occurrence of very fine silt sediment in the study area.Comparison of the inversion results and the laboratory measurements of nearby borehole samples shows that they are in general agreement.展开更多
The physical properties of silt in river reservoirs are important to river dynamics. Unfortunately, traditional techniques yield insufficient data. Based on porous media acoustic theory, we invert the acoustic paramet...The physical properties of silt in river reservoirs are important to river dynamics. Unfortunately, traditional techniques yield insufficient data. Based on porous media acoustic theory, we invert the acoustic parameters for the top river-bottom sediments. An explicit form of the acoustic reflection coefficient at the water-sediment interface is derived based on Biot's theory. The choice of parameters in the Blot model is discussed and the relation between acoustic and geological parameters is studied, including that between the reflection coefficient and porosity and the attenuation coefficient and permeability. The attenuation coefficient of the sound wave in the sediments is obtained by analyzing the shift of the signal frequency. The acoustic reflection coefficient at the water-sediment interface is extracted from the sonar signal. Thus, an inversion method of the physical parameters of the river- bottom surface sediments is proposed. The results of an experiment at the Sanmenxia reservoir suggest that the estimated grain size is close to the actual data. This demonstrates the ability of the proposed method to determine the physical parameters of sediments and estimate the grain size.展开更多
Measurements of bubble size distribution require the understanding of the acoustic characteristics of the medium.The bubbles show highly nonlinear properties under finite amplitude acoustic excitation,so the acoustic ...Measurements of bubble size distribution require the understanding of the acoustic characteristics of the medium.The bubbles show highly nonlinear properties under finite amplitude acoustic excitation,so the acoustic fields from bubble population are easily observed at the second harmonics as well as at the fundamental frequency,which shows that the nonlinear coefficient increases obviously.The inversion method of bubble size distribution based on nonlinear acoustic effects can peel off the influence of complex environment and obtain the size distribution coefficient information of bubbles more accurately.The previous nonlinear inversion methods of bubble size distribution are mostly based on the nonlinear scattering cross-section characteristics of bubbles.However,the stability of inversion is not high enough.In this paper,we introduce a new acoustic inversion method for bubble size distribution,which is based on the nonlinear coefficients of bubble medium.Compared with other inversion methods based on linear or nonlinear scattering cross section,the inversion method based on nonlinear coefficients of bubble medium proposed in this paper shows good robustness in both simulation and experiment.展开更多
This article compares the isotropic and anisotropic TV regularizations used in inverse acoustic scattering. It is observed that compared with the traditional Tikhonov regularization, isotropic and anisotropic TV regul...This article compares the isotropic and anisotropic TV regularizations used in inverse acoustic scattering. It is observed that compared with the traditional Tikhonov regularization, isotropic and anisotropic TV regularizations perform better in the sense of edge preserving. While anisotropic TV regularization will cause distortions along axes. To minimize the energy function with isotropic and anisotropic regularization terms, we use split Bregman scheme. We do several 2D numerical experiments to validate the above arguments.展开更多
The distributed source boundary point method (DSBPM) is used as the spatial transform algorithm for realizing nearfield acoustic holography (NAH), the sensitivity of the reconstructed solution to the measurement e...The distributed source boundary point method (DSBPM) is used as the spatial transform algorithm for realizing nearfield acoustic holography (NAH), the sensitivity of the reconstructed solution to the measurement errors is analyzed, and the regularization method is proposed to stabilize the reconstruction process, control the influence of the measurement errors and get a better approximate solution. An oscillating sphere is investigated as a numerical example, the influence of the measurement errors on the reconstruction solution is demonstrated, and the feasibility and validity of the regularization method are validated. Key words: Acoustic holography Boundary point method Inverse problem Regularization展开更多
It is proved that a sound-soft scatterer in R^N (N = 2, 3) is uniquely determined by a finite number of acoustic far-field measurements. The admissible scatterer possibly consists of finitely many solid obstacles an...It is proved that a sound-soft scatterer in R^N (N = 2, 3) is uniquely determined by a finite number of acoustic far-field measurements. The admissible scatterer possibly consists of finitely many solid obstacles and subsets of (N - 1)- dimensional hyperplanes.展开更多
This paper is concerned with inverse acoustic scattering in an inhomogeneous medium with a conductive boundary condition and the unknown buried impenetrable objects inside.Using a variational approach,we establish the...This paper is concerned with inverse acoustic scattering in an inhomogeneous medium with a conductive boundary condition and the unknown buried impenetrable objects inside.Using a variational approach,we establish the well-posedness of the direct problem.For the inverse problem,we shall numerically reconstruct the inhomogeneous medium from the far-field data for different kinds of cases.For the case when a Dirichlet boundary condition is imposed on the buried object,the classical factorization method proposed in[1]is justified as valid for reconstructing the inhomogeneous medium from the far-field data.For the case when a Neumann boundary condition is imposed on the buried object,the classical factorization method of[1]cannot be applied directly,since the middle operator of the factorization of the far-field operator is only compact.In this case,we develop a modified factorization method to locate the inhomogeneous medium with a conductive boundary condition and the unknown buried objects.Some numerical experiments are provided to demonstrate the practicability of the inversion algorithms developed.展开更多
The sand-conglomerate fans are the major depositional systems in the lower third member of Shahejie Formation in Shengtuo area, which formed in the deep lacustrine environment characterized by steep slope gradient, ne...The sand-conglomerate fans are the major depositional systems in the lower third member of Shahejie Formation in Shengtuo area, which formed in the deep lacustrine environment characterized by steep slope gradient, near sources and intensive tectonic activity. This work was focused on the sedimentary feature of the glutenite segment to conduct the seismic sedimentology research. The near-shore subaqueous fans and its relative gravity channel and slump turbidite fan depositions were identified according to observation and description of cores combining with the numerous data of seismic and logging. Then, the depositional model was built depending on the analysis of palaeogeomorphology. The seismic attributes which are related to the hydrocarbon but relative independent were chosen to conduct the analysis, the reservoir area of the glutenite segment was found performing a distribution where the amplitude value is relatively higher, and finally the RMS amplitude attribute was chosen to conduct the attribute predicting. At the same time, the horizontal distribution of the sedimentary facies was analyzed qualitatively. At last, the sparse spike inversion method was used to conduct the acoustic impedance inversion, and the inversion result can distinguish glutenite reservoir which is greater than 5 m. This method quantitatively characterizes the distribution area of the favorable reservoir sand.展开更多
It is better to use a simple configuration to enhance the matched-field inversion method based on a horizontal line applicability of ocean environment inversion in shallow water. A array (HLA) is used to retrieve th...It is better to use a simple configuration to enhance the matched-field inversion method based on a horizontal line applicability of ocean environment inversion in shallow water. A array (HLA) is used to retrieve the variation of sound speed profile. The performance of the inversion method is verified in the South China Sea in June, 2010. An HLA laid at bottom was used to receive signals from a bottom-mounted transducer. Inverted mean sound speed profiles from 9-hour long acoustic signals are in good agreement with measurements from two temperature chains at the sites of the source and receiver. The results show that an HLA can be used to monitor the variability of shallow-water sound speed profile.展开更多
This paper proposes a semi‐analytical and local meshless collocation method,the loca-lized method of fundamental solutions(LMFS),to address three‐dimensional(3D)acoustic inverse problems in complex domains.The propo...This paper proposes a semi‐analytical and local meshless collocation method,the loca-lized method of fundamental solutions(LMFS),to address three‐dimensional(3D)acoustic inverse problems in complex domains.The proposed approach is a recently developed numerical scheme with the potential of being mathematically simple,nu-merically accurate,and requiring less computational time and storage.In LMFS,an overdetermined sparse linear system is constructed by using the known data at the nodes on the accessible boundary and by making the remaining nodes satisfy the governing equation.In the numerical procedure,the pseudoinverse of a matrix is solved via the truncated singular value decomposition,and thus the regularization techniques are not needed in solving the resulting linear system with a well‐conditioned matrix.Numerical experiments,involving complicated geometry and the high noise level,confirm the ef-fectiveness and performance of the LMFS for solving 3D acoustic inverse problems.展开更多
基金funded by R&D Department of China National Petroleum Corporation(2022DQ0604-04)the Strategic Cooperation Technology Projects of CNPC and CUPB(ZLZX2020-03)the Science Research and Technology Development of PetroChina(2021DJ1206).
文摘Seismic impedance inversion is an important technique for structure identification and reservoir prediction.Model-based and data-driven impedance inversion are the commonly used inversion methods.In practice,the geophysical inversion problem is essentially an ill-posedness problem,which means that there are many solutions corresponding to the same seismic data.Therefore,regularization schemes,which can provide stable and unique inversion results to some extent,have been introduced into the objective function as constrain terms.Among them,given a low-frequency initial impedance model is the most commonly used regularization method,which can provide a smooth and stable solution.However,this model-based inversion method relies heavily on the initial model and the inversion result is band limited to the effective frequency bandwidth of seismic data,which cannot effectively improve the seismic vertical resolution and is difficult to be applied to complex structural regions.Therefore,we propose a data-driven approach for high-resolution impedance inversion based on the bidirectional long short-term memory recurrent neural network,which regards seismic data as time-series rather than image-like patches.Compared with the model-based inversion method,the data-driven approach provides higher resolution inversion results,which demonstrates the effectiveness of the data-driven method for recovering the high-frequency components.However,judging from the inversion results for characterization the spatial distribution of thin-layer sands,the accuracy of high-frequency components is difficult to guarantee.Therefore,we add the model constraint to the objective function to overcome the shortages of relying only on the data-driven schemes.First,constructing the supervisor1 based on the bidirectional long short-term memory recurrent neural network,which provides the predicted impedance with higher resolution.Then,convolution constraint as supervisor2 is introduced into the objective function to guarantee the reliability and accuracy of the inversion results,which makes the synthetic seismic data obtained from the inversion result consistent with the input data.Finally,we test the proposed scheme based on the synthetic and field seismic data.Compared to model-based and purely data-driven impedance inversion methods,the proposed approach provides more accurate and reliable inversion results while with higher vertical resolution and better spatial continuity.The inversion results accurately characterize the spatial distribution relationship of thin sands.The model tests demonstrate that the model-constrained and data-driven impedance inversion scheme can effectively improve the thin-layer structure characterization based on the seismic data.Moreover,tests on the oil field data indicate the practicality and adaptability of the proposed method.
基金supported by the National Key R&D Program of China (No.2021YFC2801202)the National Natural Science Foundation of China (No.42076224)the Fundamental Research Funds for the Central Universities (No.202262012)。
文摘The chirp sub-bottom profiler,for its high resolution,easy accessibility and cost-effectiveness,has been widely used in acoustic detection.In this paper,the acoustic impedance and grain size compositions were obtained based on the chirp sub-bottom profiler data collected in the Chukchi Plateau area during the 11th Arctic Expedition of China.The time-domain adaptive search matching algorithm was used and validated on our established theoretical model.The misfit between the inversion result and the theoretical model is less than 0.067%.The grain size was calculated according to the empirical relationship between the acoustic impedance and the grain size of the sediment.The average acoustic impedance of sub-seafloor strata is 2.5026×10^(6) kg(s m^(2))^(-1)and the average grain size(θvalue)of the seafloor surface sediment is 7.1498,indicating the predominant occurrence of very fine silt sediment in the study area.Comparison of the inversion results and the laboratory measurements of nearby borehole samples shows that they are in general agreement.
基金supported by the National Key R&D Program of China(Grant No.2016YFC0401608)the Scientific Fund of the Yellow River Institute for Hydraulic Research(Grant Nos.HKY-JBYW-2016-09 and HKY-JBYW-2016-29)
文摘The physical properties of silt in river reservoirs are important to river dynamics. Unfortunately, traditional techniques yield insufficient data. Based on porous media acoustic theory, we invert the acoustic parameters for the top river-bottom sediments. An explicit form of the acoustic reflection coefficient at the water-sediment interface is derived based on Biot's theory. The choice of parameters in the Blot model is discussed and the relation between acoustic and geological parameters is studied, including that between the reflection coefficient and porosity and the attenuation coefficient and permeability. The attenuation coefficient of the sound wave in the sediments is obtained by analyzing the shift of the signal frequency. The acoustic reflection coefficient at the water-sediment interface is extracted from the sonar signal. Thus, an inversion method of the physical parameters of the river- bottom surface sediments is proposed. The results of an experiment at the Sanmenxia reservoir suggest that the estimated grain size is close to the actual data. This demonstrates the ability of the proposed method to determine the physical parameters of sediments and estimate the grain size.
基金the National Natural Science Foundation of China(Grant Nos.11674074 and 61701133).
文摘Measurements of bubble size distribution require the understanding of the acoustic characteristics of the medium.The bubbles show highly nonlinear properties under finite amplitude acoustic excitation,so the acoustic fields from bubble population are easily observed at the second harmonics as well as at the fundamental frequency,which shows that the nonlinear coefficient increases obviously.The inversion method of bubble size distribution based on nonlinear acoustic effects can peel off the influence of complex environment and obtain the size distribution coefficient information of bubbles more accurately.The previous nonlinear inversion methods of bubble size distribution are mostly based on the nonlinear scattering cross-section characteristics of bubbles.However,the stability of inversion is not high enough.In this paper,we introduce a new acoustic inversion method for bubble size distribution,which is based on the nonlinear coefficients of bubble medium.Compared with other inversion methods based on linear or nonlinear scattering cross section,the inversion method based on nonlinear coefficients of bubble medium proposed in this paper shows good robustness in both simulation and experiment.
文摘This article compares the isotropic and anisotropic TV regularizations used in inverse acoustic scattering. It is observed that compared with the traditional Tikhonov regularization, isotropic and anisotropic TV regularizations perform better in the sense of edge preserving. While anisotropic TV regularization will cause distortions along axes. To minimize the energy function with isotropic and anisotropic regularization terms, we use split Bregman scheme. We do several 2D numerical experiments to validate the above arguments.
基金This project is supported by National Natural Science Foundation of China(No.50275044)Research Fund for Doctoral Program of Ministry of Education of China(No.20020359005).
文摘The distributed source boundary point method (DSBPM) is used as the spatial transform algorithm for realizing nearfield acoustic holography (NAH), the sensitivity of the reconstructed solution to the measurement errors is analyzed, and the regularization method is proposed to stabilize the reconstruction process, control the influence of the measurement errors and get a better approximate solution. An oscillating sphere is investigated as a numerical example, the influence of the measurement errors on the reconstruction solution is demonstrated, and the feasibility and validity of the regularization method are validated. Key words: Acoustic holography Boundary point method Inverse problem Regularization
文摘It is proved that a sound-soft scatterer in R^N (N = 2, 3) is uniquely determined by a finite number of acoustic far-field measurements. The admissible scatterer possibly consists of finitely many solid obstacles and subsets of (N - 1)- dimensional hyperplanes.
基金supported by the National Natural ScienceFoundation of China Grant(11871416,12171057)the Natural Science Foundation of Shandong Province Grant(ZR2019MA027)。
文摘This paper is concerned with inverse acoustic scattering in an inhomogeneous medium with a conductive boundary condition and the unknown buried impenetrable objects inside.Using a variational approach,we establish the well-posedness of the direct problem.For the inverse problem,we shall numerically reconstruct the inhomogeneous medium from the far-field data for different kinds of cases.For the case when a Dirichlet boundary condition is imposed on the buried object,the classical factorization method proposed in[1]is justified as valid for reconstructing the inhomogeneous medium from the far-field data.For the case when a Neumann boundary condition is imposed on the buried object,the classical factorization method of[1]cannot be applied directly,since the middle operator of the factorization of the far-field operator is only compact.In this case,we develop a modified factorization method to locate the inhomogeneous medium with a conductive boundary condition and the unknown buried objects.Some numerical experiments are provided to demonstrate the practicability of the inversion algorithms developed.
基金Project(41172109)supported by the National Natural Science Foundation of ChinaProject(20110003110014)supported by the ResearchFoundation for the Doctoral Program of Higher Education,China
文摘The sand-conglomerate fans are the major depositional systems in the lower third member of Shahejie Formation in Shengtuo area, which formed in the deep lacustrine environment characterized by steep slope gradient, near sources and intensive tectonic activity. This work was focused on the sedimentary feature of the glutenite segment to conduct the seismic sedimentology research. The near-shore subaqueous fans and its relative gravity channel and slump turbidite fan depositions were identified according to observation and description of cores combining with the numerous data of seismic and logging. Then, the depositional model was built depending on the analysis of palaeogeomorphology. The seismic attributes which are related to the hydrocarbon but relative independent were chosen to conduct the analysis, the reservoir area of the glutenite segment was found performing a distribution where the amplitude value is relatively higher, and finally the RMS amplitude attribute was chosen to conduct the attribute predicting. At the same time, the horizontal distribution of the sedimentary facies was analyzed qualitatively. At last, the sparse spike inversion method was used to conduct the acoustic impedance inversion, and the inversion result can distinguish glutenite reservoir which is greater than 5 m. This method quantitatively characterizes the distribution area of the favorable reservoir sand.
基金supported by the National Natural Science Foundation of China(Grant Nos.11434012,11404366,11125420 and 11074269)
文摘It is better to use a simple configuration to enhance the matched-field inversion method based on a horizontal line applicability of ocean environment inversion in shallow water. A array (HLA) is used to retrieve the variation of sound speed profile. The performance of the inversion method is verified in the South China Sea in June, 2010. An HLA laid at bottom was used to receive signals from a bottom-mounted transducer. Inverted mean sound speed profiles from 9-hour long acoustic signals are in good agreement with measurements from two temperature chains at the sites of the source and receiver. The results show that an HLA can be used to monitor the variability of shallow-water sound speed profile.
基金National Natural Science Foundation of China,Grant/Award Number:11802151Natural Science Foundation of Shandong Province of China,Grant/Award Number:ZR2019BA008+1 种基金supported by the National Natural Science Foundation of China(No.11802151)the Natural Science Foundation of Shandong Province of China(No.ZR2019BA008).
文摘This paper proposes a semi‐analytical and local meshless collocation method,the loca-lized method of fundamental solutions(LMFS),to address three‐dimensional(3D)acoustic inverse problems in complex domains.The proposed approach is a recently developed numerical scheme with the potential of being mathematically simple,nu-merically accurate,and requiring less computational time and storage.In LMFS,an overdetermined sparse linear system is constructed by using the known data at the nodes on the accessible boundary and by making the remaining nodes satisfy the governing equation.In the numerical procedure,the pseudoinverse of a matrix is solved via the truncated singular value decomposition,and thus the regularization techniques are not needed in solving the resulting linear system with a well‐conditioned matrix.Numerical experiments,involving complicated geometry and the high noise level,confirm the ef-fectiveness and performance of the LMFS for solving 3D acoustic inverse problems.