Since its inception in the 1970s,multi-dimensional magnetic resonance(MR)has emerged as a powerful tool for non-invasive investigations of structures and molecular interactions.MR spectroscopy beyond one dimension all...Since its inception in the 1970s,multi-dimensional magnetic resonance(MR)has emerged as a powerful tool for non-invasive investigations of structures and molecular interactions.MR spectroscopy beyond one dimension allows the study of the correlation,exchange processes,and separation of overlapping spectral information.The multi-dimensional concept has been re-implemented over the last two decades to explore molecular motion and spin dynamics in porous media.Apart from Fourier transform,methods have been developed for processing the multi-dimensional time-domain data,identifying the fluid components,and estimating pore surface permeability via joint relaxation and diffusion spectra.Through the resolution of spectroscopic signals with spatial encoding gradients,multi-dimensional MR imaging has been widely used to investigate the microscopic environment of living tissues and distinguish diseases.Signals in each voxel are usually expressed as multi-exponential decay,representing microstructures or environments along multiple pore scales.The separation of contributions from different environments is a common ill-posed problem,which can be resolved numerically.Moreover,the inversion methods and experimental parameters determine the resolution of multi-dimensional spectra.This paper reviews the algorithms that have been proposed to process multidimensional MR datasets in different scenarios.Detailed information at the microscopic level,such as tissue components,fluid types and food structures in multi-disciplinary sciences,could be revealed through multi-dimensional MR.展开更多
Tikhonov regularization（TR） method has played a very important role in the gravity data and magnetic data process. In this paper, the Tikhonov regularization method with respect to the inversion of gravity data is d...Tikhonov regularization（TR） method has played a very important role in the gravity data and magnetic data process. In this paper, the Tikhonov regularization method with respect to the inversion of gravity data is discussed. and the extrapolated TR method（EXTR） is introduced to improve the fitting error. Furthermore, the effect of the parameters in the EXTR method on the fitting error, number of iterations, and inversion results are discussed in details. The computation results using a synthetic model with the same and different densities indicated that. compared with the TR method, the EXTR method not only achieves the a priori fitting error level set by the interpreter but also increases the fitting precision, although it increases the computation time and number of iterations. And the EXTR inversion results are more compact than the TR inversion results, which are more divergent. The range of the inversion data is closer to the default range of the model parameters, and the model features and default model density distribution agree well.展开更多
With the continuous development of full tensor gradiometer （FTG） measurement techniques, three-dimensional （3D） inversion of FTG data is becoming increasingly used in oil and gas exploration. In the fast processin...With the continuous development of full tensor gradiometer （FTG） measurement techniques, three-dimensional （3D） inversion of FTG data is becoming increasingly used in oil and gas exploration. In the fast processing and interpretation of large-scale high-precision data, the use of the graphics processing unit process unit （GPU） and preconditioning methods are very important in the data inversion. In this paper, an improved preconditioned conjugate gradient algorithm is proposed by combining the symmetric successive over-relaxation （SSOR） technique and the incomplete Choleksy decomposition conjugate gradient algorithm （ICCG）. Since preparing the preconditioner requires extra time, a parallel implement based on GPU is proposed. The improved method is then applied in the inversion of noise- contaminated synthetic data to prove its adaptability in the inversion of 3D FTG data. Results show that the parallel SSOR-ICCG algorithm based on NVIDIA Tesla C2050 GPU achieves a speedup of approximately 25 times that of a serial program using a 2.0 GHz Central Processing Unit （CPU）. Real airbome gravity-gradiometry data from Vinton salt dome （south- west Louisiana, USA） are also considered. Good results are obtained, which verifies the efficiency and feasibility of the proposed parallel method in fast inversion of 3D FTG data.展开更多
Conventional time-space domain and frequency-space domain prediction filtering methods assume that seismic data consists of two parts, signal and random noise. That is, the so-called additive noise model. However, whe...Conventional time-space domain and frequency-space domain prediction filtering methods assume that seismic data consists of two parts, signal and random noise. That is, the so-called additive noise model. However, when estimating random noise, it is assumed that random noise can be predicted from the seismic data by convolving with a prediction error filter. That is, the source-noise model. Model inconsistencies, before and after denoising, compromise the noise attenuation and signal-preservation performances of prediction filtering methods. Therefore, this study presents an inversion-based time-space domain random noise attenuation method to overcome the model inconsistencies. In this method, a prediction error filter （PEF）, is first estimated from seismic data; the filter characterizes the predictability of the seismic data and adaptively describes the seismic data＇s space structure. After calculating PEF, it can be applied as a regularized constraint in the inversion process for seismic signal from noisy data. Unlike conventional random noise attenuation methods, the proposed method solves a seismic data inversion problem using regularization constraint; this overcomes the model inconsistency of the prediction filtering method. The proposed method was tested on both synthetic and real seismic data, and results from the prediction filtering method and the proposed method are compared. The testing demonstrated that the proposed method suppresses noise effectively and provides better signal-preservation performance.展开更多
A data-space inversion(DSI)method has been recently proposed and successfully applied to the history matching and production prediction of reservoirs.Based on Bayesian theory,DSI can directly and effectively obtain go...A data-space inversion(DSI)method has been recently proposed and successfully applied to the history matching and production prediction of reservoirs.Based on Bayesian theory,DSI can directly and effectively obtain good posterior flow predictions without inversion of geological parameters of reservoir model.This paper presents an improved DSI method to fast predict reservoir state fields(e.g.saturation and pressure profiles)via observed production data.Firstly,a large number of production curves and state data are generated by reservoir model simulation to expand the data space of original DSI.Then,efficient history matching only on the observed production data is carried out via the original DSI to obtain related parameters which reflects the weight of the real reservoir model relative to prior reservoir models.Finally,those parameters are used to predict the oil saturation and pressure profiles of the real reservoir model by combining large amounts of state data of prior reservoir models.Two examples including conventional heterogeneous and unconventional fractured reservoir are implemented to test the performances of predicting saturation and pressure profiles of this improved DSI method.Besides,this method is also tested in a real field and the obtained results show the high computational efficiency and high accuracy of the practical application of this method.展开更多
Seismic inversion and basic theory are briefly presented and the main idea of this method is introduced. Both non-linear wave equation inversion technique and Complete Utilization of Samples Information (CUSI) neural ...Seismic inversion and basic theory are briefly presented and the main idea of this method is introduced. Both non-linear wave equation inversion technique and Complete Utilization of Samples Information (CUSI) neural network analysis are used in lithological interpretation in Jibei coal field. The prediction results indicate that this method can provide reliable data for thin coal exploitation and promising area evaluation.展开更多
Based on the synchronous joint gravity and magnetic inversion of single interface by Pilkington and the need of revealing Cenozoic and crystalline basement thickness in the new round of oil-gas exploration, we propose...Based on the synchronous joint gravity and magnetic inversion of single interface by Pilkington and the need of revealing Cenozoic and crystalline basement thickness in the new round of oil-gas exploration, we propose a joint gravity and magnetic inversion methodfor two-layer models by concentrating on the relationship between the change of thicknessI and position of the middle layer and anomaly and discuss the effects of the key parameters. Model tests and application to field data show the validity of this method.展开更多
A two-dimensional forward and backward algorithm for the controlled-source audio-frequency magnetotelluric （CSAMT） method is developed to invert data in the entire region （near, transition, and far） and deal with ...A two-dimensional forward and backward algorithm for the controlled-source audio-frequency magnetotelluric （CSAMT） method is developed to invert data in the entire region （near, transition, and far） and deal with the effects of artificial sources. First, a regularization factor is introduced in the 2D magnetic inversion, and the magnetic susceptibility is updated in logarithmic form so that the inversion magnetic susceptibility is always positive. Second, the joint inversion of the CSAMT and magnetic methods is completed with the introduction of the cross gradient. By searching for the weight of the cross-gradient term in the objective function, the mutual influence between two different physical properties at different locations are avoided. Model tests show that the joint inversion based on cross-gradient theory offers better results than the single-method inversion. The 2D forward and inverse algorithm for CSAMT with source can effectively deal with artificial sources and ensures the reliability of the final joint inversion algorithm.展开更多
Geodetic functional models,stochastic models,and model parameter estimation theory are fundamental for geodetic data processing.In the past five years,through the unremitting efforts of Chinese scholars in the field o...Geodetic functional models,stochastic models,and model parameter estimation theory are fundamental for geodetic data processing.In the past five years,through the unremitting efforts of Chinese scholars in the field of geodetic data processing,according to the application and practice of geodesy,they have made significant contributions in the fields of hypothesis testing theory,un-modeled error,outlier detection,and robust estimation,variance component estimation,complex least squares,and ill-posed problems treatment.Many functional models such as the nonlinear adjustment model,EIV model,and mixed additive and multiplicative random error model are also constructed and improved.Geodetic data inversion is an important part of geodetic data processing,and Chinese scholars have done a lot of work in geodetic data inversion in the past five years,such as seismic slide distribution inversion,intelligent inversion algorithm,multi-source data joint inversion,water reserve change and satellite gravity inversion.This paper introduces the achievements of Chinese scholars in the field of geodetic data processing in the past five years,analyzes the methods used by scholars and the problems solved,and looks forward to the unsolved problems in geodetic data processing and the direction that needs further research in the future.展开更多
In direct sequence spread spectrum communication both for satelliteto-ground and inter-satellite links, the system constrains due to radio frequency spectral occupation, channel data throughput and link performances i...In direct sequence spread spectrum communication both for satelliteto-ground and inter-satellite links, the system constrains due to radio frequency spectral occupation, channel data throughput and link performances in terms of data channel coding which might result in a signal structure where the symbol duration is shorter than the pseudo code period. This can generate some difficulties in the DSSS signal acquisition due to the polarity inversion caused by the data modulation. To eliminate the influence due to polarity inversion, this paper proposes a novel acquisition algorithm based on the simultaneous search of the code phase, data phase and Doppler frequency. In the proposed algorithm the data phase is predicted and the correlation period for the coherent integration can be set equal to the symbol duration. Then non-coherent accumulation over different symbol is implemented in order to enhance the acquisition algorithm sensitivity; the interval of non-coherent accumulation is the least common multiple between the symbol duration and the pseudo code period. The algorithm proposed can largely minimize the SNR loss caused by data polarity inversion and enhance acquisition performance without a noticeable increase in hardware complexity. Theoretical analysis, simulation and measured results verify the validity of the algorithm.展开更多
Based on surfaced-related multiple elimination （SRME） , this research has derived the methods on multiples elimination in the inverse data space. Inverse data processing means moving seismic data from forwar...Based on surfaced-related multiple elimination （SRME） , this research has derived the methods on multiples elimination in the inverse data space. Inverse data processing means moving seismic data from forward data space （FDS） to inverse data space （ IDS） . The surface-related multiples and primaries can then be sepa-rated in the IDS, since surface-related multiples wi l l form a focus region in the IDS. Muting the multiples ener-gy can achieve the purpose of multiples elimination and avoid the damage to primaries energy during the process of adaptive subtraction. Randomized singular value decomposition （ RSYD） is used to enhance calculation speed and improve the accuracy in the conversion of FDS to IDS. The synthetic shot record of the salt dome model shows that the relationship between primaries and multiples is simple and clear, and RSVD can easily eliminate multiples and save primaries energy. Compared with conventional multiples elimination methods and ordinary methods of multiples elimination in the inverse data space, this technique has an advantage of high cal-culation speed and reliable outcomes.展开更多
In order to accurately evaluate the fracturing stimulation effect of horizontal wells in shale gas reservoirs,it is suggested to establish a model to predict the temperature distribution in the wellbore of horizontal ...In order to accurately evaluate the fracturing stimulation effect of horizontal wells in shale gas reservoirs,it is suggested to establish a model to predict the temperature distribution in the wellbore of horizontal well and an inversion model of DTS data by using MCMC algorithm,and optimize the production profile interpretation process.In this paper,the characteristics of temperature profile of fractured horizontal wells in shale gas reservoirs were analyzed and the main factors affecting the temperature profile were figured out.Finally,the newly established inversion model was applied to the production profile interpretation of one case well in a certain shale gas reservoir.And the following research results were obtained.First,the temperature profile of fractured horizontal wells is in the shape of irregular“saw tooth”,and each“saw tooth”corresponds to an effective hydraulic fracture with fluid inflow.Second,the longer the fracture is,the greater the wellbore temperature drop at the corresponding fracture location is,and the gas flow rate in the fracture is positively correlated with the temperature drop.Third,from the perspective of influence degree,the factors influencing the temperature profile of fractured horizontal wells in shale gas reservoirs are ranked from the strong to the weak as follows:fracture half-length,gas flow rate,permeability of stimulated area,wellbore diameter,fracture conductivity,horizontal inclination angle,and comprehensive thermal conductivity,among which,the first three are main controlling factors.Fourth,the MCMC inversion method is applied to invert the DST temperature data of the case well.The temperature profile predicted in the model is better accordant with the measured DTS profile,and the absolute error of the predicted temperature at different levels of effective hydraulic fractures is less than 0.02°C.The interpreted gas flow rate of each fracturing stage is closer to the field measurement,and the deviation of the maximum gas flow rate of a single fracturing stage is only 180.35 m^(3)/d.The absolute error between single-well gas production rate and gas production rate measured at the wellhead is less than 3 m^(3)/d,which proves the reliability of this newly developed inversion model.展开更多
In order to explore the nonlinear structure hidden in high-dimensional data, some dimen-sion reduction techniques have been developed, such as the Projection Pursuit technique (PP).However, PP will involve enormous co...In order to explore the nonlinear structure hidden in high-dimensional data, some dimen-sion reduction techniques have been developed, such as the Projection Pursuit technique (PP).However, PP will involve enormous computational load. To overcome this, an inverse regressionmethod is proposed. In this paper, we discuss and develop this method. To seek the interestingprojective direction, the minimization of the residual sum of squares is used as a criterion, andspline functions are applied to approximate the general nonlinear transform function. The algo-rithm is simple, and saves the computational load. Under certain proper conditions, consistencyof the estimators of the interesting direction is shown.展开更多
A simple data assimilation method for improving estimation of moderate resolution imaging spectroradiometer （MODIS） leaf area index （LAI） time-series data products based on the gradient inverse weighted filter and...A simple data assimilation method for improving estimation of moderate resolution imaging spectroradiometer （MODIS） leaf area index （LAI） time-series data products based on the gradient inverse weighted filter and object analysis is proposed. The properties and quality control （QC） of MODIS LAI data products are introduced. Also, the gradient inverse weighted filter and object analysis are analyzed. An experiment based on the simple data assimilation method is performed using MODIS LAI data sets from 2000 to 2005 of Guizhou Province in China.展开更多
基金supported by the National Natural Science Foundation of China(No.61901465,82222032,82172050).
文摘Since its inception in the 1970s,multi-dimensional magnetic resonance(MR)has emerged as a powerful tool for non-invasive investigations of structures and molecular interactions.MR spectroscopy beyond one dimension allows the study of the correlation,exchange processes,and separation of overlapping spectral information.The multi-dimensional concept has been re-implemented over the last two decades to explore molecular motion and spin dynamics in porous media.Apart from Fourier transform,methods have been developed for processing the multi-dimensional time-domain data,identifying the fluid components,and estimating pore surface permeability via joint relaxation and diffusion spectra.Through the resolution of spectroscopic signals with spatial encoding gradients,multi-dimensional MR imaging has been widely used to investigate the microscopic environment of living tissues and distinguish diseases.Signals in each voxel are usually expressed as multi-exponential decay,representing microstructures or environments along multiple pore scales.The separation of contributions from different environments is a common ill-posed problem,which can be resolved numerically.Moreover,the inversion methods and experimental parameters determine the resolution of multi-dimensional spectra.This paper reviews the algorithms that have been proposed to process multidimensional MR datasets in different scenarios.Detailed information at the microscopic level,such as tissue components,fluid types and food structures in multi-disciplinary sciences,could be revealed through multi-dimensional MR.
基金supported by the National Scientific and Technological Plan(Nos.2009BAB43B00 and 2009BAB43B01)
文摘Tikhonov regularization（TR） method has played a very important role in the gravity data and magnetic data process. In this paper, the Tikhonov regularization method with respect to the inversion of gravity data is discussed. and the extrapolated TR method（EXTR） is introduced to improve the fitting error. Furthermore, the effect of the parameters in the EXTR method on the fitting error, number of iterations, and inversion results are discussed in details. The computation results using a synthetic model with the same and different densities indicated that. compared with the TR method, the EXTR method not only achieves the a priori fitting error level set by the interpreter but also increases the fitting precision, although it increases the computation time and number of iterations. And the EXTR inversion results are more compact than the TR inversion results, which are more divergent. The range of the inversion data is closer to the default range of the model parameters, and the model features and default model density distribution agree well.
基金the Sub-project of National Science and Technology Major Project of China(No.2016ZX05027-002-003)the National Natural Science Foundation of China(No.41404089)+1 种基金the State Key Program of National Natural Science of China(No.41430322)the National Basic Research Program of China(973 Program)(No.2015CB45300)
文摘With the continuous development of full tensor gradiometer （FTG） measurement techniques, three-dimensional （3D） inversion of FTG data is becoming increasingly used in oil and gas exploration. In the fast processing and interpretation of large-scale high-precision data, the use of the graphics processing unit process unit （GPU） and preconditioning methods are very important in the data inversion. In this paper, an improved preconditioned conjugate gradient algorithm is proposed by combining the symmetric successive over-relaxation （SSOR） technique and the incomplete Choleksy decomposition conjugate gradient algorithm （ICCG）. Since preparing the preconditioner requires extra time, a parallel implement based on GPU is proposed. The improved method is then applied in the inversion of noise- contaminated synthetic data to prove its adaptability in the inversion of 3D FTG data. Results show that the parallel SSOR-ICCG algorithm based on NVIDIA Tesla C2050 GPU achieves a speedup of approximately 25 times that of a serial program using a 2.0 GHz Central Processing Unit （CPU）. Real airbome gravity-gradiometry data from Vinton salt dome （south- west Louisiana, USA） are also considered. Good results are obtained, which verifies the efficiency and feasibility of the proposed parallel method in fast inversion of 3D FTG data.
基金supported by the National Natural Science Foundation of China(No.41474109)the China National Petroleum Corporation under grant number 2016A-33
文摘Conventional time-space domain and frequency-space domain prediction filtering methods assume that seismic data consists of two parts, signal and random noise. That is, the so-called additive noise model. However, when estimating random noise, it is assumed that random noise can be predicted from the seismic data by convolving with a prediction error filter. That is, the source-noise model. Model inconsistencies, before and after denoising, compromise the noise attenuation and signal-preservation performances of prediction filtering methods. Therefore, this study presents an inversion-based time-space domain random noise attenuation method to overcome the model inconsistencies. In this method, a prediction error filter （PEF）, is first estimated from seismic data; the filter characterizes the predictability of the seismic data and adaptively describes the seismic data＇s space structure. After calculating PEF, it can be applied as a regularized constraint in the inversion process for seismic signal from noisy data. Unlike conventional random noise attenuation methods, the proposed method solves a seismic data inversion problem using regularization constraint; this overcomes the model inconsistency of the prediction filtering method. The proposed method was tested on both synthetic and real seismic data, and results from the prediction filtering method and the proposed method are compared. The testing demonstrated that the proposed method suppresses noise effectively and provides better signal-preservation performance.
基金supported by Southern Marine Science and Engineering Guangdong Laboratory(Zhanjiang)(No.ZJW-2019-04)Cooperative Innovation Center of Unconventional Oil and Gas(Ministry of Education&Hubei Province),Yangtze University(No.UOG2020-17)the National Natural Science Foundation of China(No.51874044,51922007)。
文摘A data-space inversion(DSI)method has been recently proposed and successfully applied to the history matching and production prediction of reservoirs.Based on Bayesian theory,DSI can directly and effectively obtain good posterior flow predictions without inversion of geological parameters of reservoir model.This paper presents an improved DSI method to fast predict reservoir state fields(e.g.saturation and pressure profiles)via observed production data.Firstly,a large number of production curves and state data are generated by reservoir model simulation to expand the data space of original DSI.Then,efficient history matching only on the observed production data is carried out via the original DSI to obtain related parameters which reflects the weight of the real reservoir model relative to prior reservoir models.Finally,those parameters are used to predict the oil saturation and pressure profiles of the real reservoir model by combining large amounts of state data of prior reservoir models.Two examples including conventional heterogeneous and unconventional fractured reservoir are implemented to test the performances of predicting saturation and pressure profiles of this improved DSI method.Besides,this method is also tested in a real field and the obtained results show the high computational efficiency and high accuracy of the practical application of this method.
文摘Seismic inversion and basic theory are briefly presented and the main idea of this method is introduced. Both non-linear wave equation inversion technique and Complete Utilization of Samples Information (CUSI) neural network analysis are used in lithological interpretation in Jibei coal field. The prediction results indicate that this method can provide reliable data for thin coal exploitation and promising area evaluation.
基金Supported by the National Natural Science Foundation of China(Grant No.40674063)National Hi-tech Research and Development Program of China(863Program)(Grant No.2006AA09Z311)
文摘Based on the synchronous joint gravity and magnetic inversion of single interface by Pilkington and the need of revealing Cenozoic and crystalline basement thickness in the new round of oil-gas exploration, we propose a joint gravity and magnetic inversion methodfor two-layer models by concentrating on the relationship between the change of thicknessI and position of the middle layer and anomaly and discuss the effects of the key parameters. Model tests and application to field data show the validity of this method.
基金jointly sponsored by the Fundamental Research Funds for the Central Universitiesthe National Natural Science Foundation of China(No.41374078)
文摘A two-dimensional forward and backward algorithm for the controlled-source audio-frequency magnetotelluric （CSAMT） method is developed to invert data in the entire region （near, transition, and far） and deal with the effects of artificial sources. First, a regularization factor is introduced in the 2D magnetic inversion, and the magnetic susceptibility is updated in logarithmic form so that the inversion magnetic susceptibility is always positive. Second, the joint inversion of the CSAMT and magnetic methods is completed with the introduction of the cross gradient. By searching for the weight of the cross-gradient term in the objective function, the mutual influence between two different physical properties at different locations are avoided. Model tests show that the joint inversion based on cross-gradient theory offers better results than the single-method inversion. The 2D forward and inverse algorithm for CSAMT with source can effectively deal with artificial sources and ensures the reliability of the final joint inversion algorithm.
基金National Natural Science Foundation of China(No.42174011)。
文摘Geodetic functional models,stochastic models,and model parameter estimation theory are fundamental for geodetic data processing.In the past five years,through the unremitting efforts of Chinese scholars in the field of geodetic data processing,according to the application and practice of geodesy,they have made significant contributions in the fields of hypothesis testing theory,un-modeled error,outlier detection,and robust estimation,variance component estimation,complex least squares,and ill-posed problems treatment.Many functional models such as the nonlinear adjustment model,EIV model,and mixed additive and multiplicative random error model are also constructed and improved.Geodetic data inversion is an important part of geodetic data processing,and Chinese scholars have done a lot of work in geodetic data inversion in the past five years,such as seismic slide distribution inversion,intelligent inversion algorithm,multi-source data joint inversion,water reserve change and satellite gravity inversion.This paper introduces the achievements of Chinese scholars in the field of geodetic data processing in the past five years,analyzes the methods used by scholars and the problems solved,and looks forward to the unsolved problems in geodetic data processing and the direction that needs further research in the future.
基金the support of the National High Technology Research and Development Program of China (863) (Grant No. 2012AA1406)
文摘In direct sequence spread spectrum communication both for satelliteto-ground and inter-satellite links, the system constrains due to radio frequency spectral occupation, channel data throughput and link performances in terms of data channel coding which might result in a signal structure where the symbol duration is shorter than the pseudo code period. This can generate some difficulties in the DSSS signal acquisition due to the polarity inversion caused by the data modulation. To eliminate the influence due to polarity inversion, this paper proposes a novel acquisition algorithm based on the simultaneous search of the code phase, data phase and Doppler frequency. In the proposed algorithm the data phase is predicted and the correlation period for the coherent integration can be set equal to the symbol duration. Then non-coherent accumulation over different symbol is implemented in order to enhance the acquisition algorithm sensitivity; the interval of non-coherent accumulation is the least common multiple between the symbol duration and the pseudo code period. The algorithm proposed can largely minimize the SNR loss caused by data polarity inversion and enhance acquisition performance without a noticeable increase in hardware complexity. Theoretical analysis, simulation and measured results verify the validity of the algorithm.
文摘Based on surfaced-related multiple elimination （SRME） , this research has derived the methods on multiples elimination in the inverse data space. Inverse data processing means moving seismic data from forward data space （FDS） to inverse data space （ IDS） . The surface-related multiples and primaries can then be sepa-rated in the IDS, since surface-related multiples wi l l form a focus region in the IDS. Muting the multiples ener-gy can achieve the purpose of multiples elimination and avoid the damage to primaries energy during the process of adaptive subtraction. Randomized singular value decomposition （ RSYD） is used to enhance calculation speed and improve the accuracy in the conversion of FDS to IDS. The synthetic shot record of the salt dome model shows that the relationship between primaries and multiples is simple and clear, and RSVD can easily eliminate multiples and save primaries energy. Compared with conventional multiples elimination methods and ordinary methods of multiples elimination in the inverse data space, this technique has an advantage of high cal-culation speed and reliable outcomes.
基金supported by the National Major Science and Technology Project“Comprehensive evaluation of horizontal well completion and key technique with high-temperature-resistant and anti-migration slurry”(No.:2016ZX 05021-005-009HZ).
文摘In order to accurately evaluate the fracturing stimulation effect of horizontal wells in shale gas reservoirs,it is suggested to establish a model to predict the temperature distribution in the wellbore of horizontal well and an inversion model of DTS data by using MCMC algorithm,and optimize the production profile interpretation process.In this paper,the characteristics of temperature profile of fractured horizontal wells in shale gas reservoirs were analyzed and the main factors affecting the temperature profile were figured out.Finally,the newly established inversion model was applied to the production profile interpretation of one case well in a certain shale gas reservoir.And the following research results were obtained.First,the temperature profile of fractured horizontal wells is in the shape of irregular“saw tooth”,and each“saw tooth”corresponds to an effective hydraulic fracture with fluid inflow.Second,the longer the fracture is,the greater the wellbore temperature drop at the corresponding fracture location is,and the gas flow rate in the fracture is positively correlated with the temperature drop.Third,from the perspective of influence degree,the factors influencing the temperature profile of fractured horizontal wells in shale gas reservoirs are ranked from the strong to the weak as follows:fracture half-length,gas flow rate,permeability of stimulated area,wellbore diameter,fracture conductivity,horizontal inclination angle,and comprehensive thermal conductivity,among which,the first three are main controlling factors.Fourth,the MCMC inversion method is applied to invert the DST temperature data of the case well.The temperature profile predicted in the model is better accordant with the measured DTS profile,and the absolute error of the predicted temperature at different levels of effective hydraulic fractures is less than 0.02°C.The interpreted gas flow rate of each fracturing stage is closer to the field measurement,and the deviation of the maximum gas flow rate of a single fracturing stage is only 180.35 m^(3)/d.The absolute error between single-well gas production rate and gas production rate measured at the wellhead is less than 3 m^(3)/d,which proves the reliability of this newly developed inversion model.
基金This project is supported by the National Natural Science Foundation of China
文摘In order to explore the nonlinear structure hidden in high-dimensional data, some dimen-sion reduction techniques have been developed, such as the Projection Pursuit technique (PP).However, PP will involve enormous computational load. To overcome this, an inverse regressionmethod is proposed. In this paper, we discuss and develop this method. To seek the interestingprojective direction, the minimization of the residual sum of squares is used as a criterion, andspline functions are applied to approximate the general nonlinear transform function. The algo-rithm is simple, and saves the computational load. Under certain proper conditions, consistencyof the estimators of the interesting direction is shown.
基金This work was supported by the China Postdoctoral Science Foundation(No.20060390326)the key international S&T cooperation project of China(No.2004DFA06300).
文摘A simple data assimilation method for improving estimation of moderate resolution imaging spectroradiometer （MODIS） leaf area index （LAI） time-series data products based on the gradient inverse weighted filter and object analysis is proposed. The properties and quality control （QC） of MODIS LAI data products are introduced. Also, the gradient inverse weighted filter and object analysis are analyzed. An experiment based on the simple data assimilation method is performed using MODIS LAI data sets from 2000 to 2005 of Guizhou Province in China.
