The 2D data processing adopted by the high-density resistivity method regards the geological structures as two degrees, which makes the results of the 2D data inversion only an approximate interpretation;the accuracy ...The 2D data processing adopted by the high-density resistivity method regards the geological structures as two degrees, which makes the results of the 2D data inversion only an approximate interpretation;the accuracy and effect can not meet the precise requirement of the inversion. Two typical models of the geological bodies were designed, and forward calculation was carried out using finite element method. The forward-modeled profiles were obtained. 1% Gaussian random error was added in the forward models and then 2D and 3D inversions using a high-density resistivity method were undertaken to realistically simulate field data and analyze the sensitivity of the 2D and 3D inversion algorithms to noise. Contrast between the 2D and 3D inversion results of least squares inversion shows that two inversion results of high-density resistivity method all can basically reflect the spatial position of an anomalous body. However, the 3D inversion can more effectively eliminate the influence of interference from Gaussian random error and better reflect the distribution of resistivity in the anomalous bodies. Overall, the 3D inversion was better than 2D inversion in terms of embodying anomalous body positions, morphology and resistivity properties.展开更多
The synthesis of new 4-imino-4H-chromeno[2,3-d]pyrimidin-3(5H)-amine in four steps including one step under microwave dielectric heating is reported. The structural identity of the synthesized compounds was establishe...The synthesis of new 4-imino-4H-chromeno[2,3-d]pyrimidin-3(5H)-amine in four steps including one step under microwave dielectric heating is reported. The structural identity of the synthesized compounds was established according to their spectroscopic analysis, such as FT-IR, NMR and mass spectroscopy. These new compounds were tested for their antiproliferative activities on seven representative human tumoral cell lines (Huh7 D12, Caco2, MDA-MB231, MDA-MB468, HCT116, PC3 and MCF7) and also on fibroblasts. Among them, only the compounds 6c showed micromolar cytotoxic activity on tumor cell lines (1.8 50 50 > 25 μM). Finally, in silico ADMET studies ware performed to investigate the possibility of using of the identified compound 6c as potential anti-tumor compound.展开更多
An engineering system approach of 2-D cylindrical model of transient mass balance calculations of ozone and other concerned chemicals along with fourteen photolysis, ozone-generating and ozone-depleting chemical react...An engineering system approach of 2-D cylindrical model of transient mass balance calculations of ozone and other concerned chemicals along with fourteen photolysis, ozone-generating and ozone-depleting chemical reaction equations was developed, validated, and used for studying the ozone concentrations, distribution and peak of the layer, ozone depletion and total ozone abundance in the stratosphere. The calculated ozone concentrations and profile at both the Equator and a 60˚N location were found to follow closely with the measured data. The calculated average ozone concentration was within 1% of the measured average, and the deviation of ozone profiles was within 14%. The monthly evolution of stratospheric ozone concentrations and distribution above the Equator was studied with results discussed in details. The influences of slow air movement in both altitudinal and radial directions on ozone concentrations and profile in the stratosphere were explored and discussed. Parametric studies of the influences of gas diffusivities of ozone D<sub>O3</sub> and active atomic oxygen D<sub>O</sub> on ozone concentrations and distributions were also studied and delineated. Having both influences through physical diffusion and chemical reactions, the diffusivity (and diffusion) of atomic oxygen D<sub>O</sub> was found to be more sensitive and important than that of ozone D<sub>O3</sub> on ozone concentrations and distribution. The 2-D ozone model present in this paper for stratospheric ozone and its layer and depletion is shown to be robust, convenient, efficient, and executable for analyzing the complex ozone phenomena in the stratosphere. .展开更多
2-D velocity structure and tectonics of the crust and upper mantle is revealed by inversion of seismic refraction and wide-angle reflection traveltimes acquired along the profile L1 in the Changbaishan-Tianchi volcani...2-D velocity structure and tectonics of the crust and upper mantle is revealed by inversion of seismic refraction and wide-angle reflection traveltimes acquired along the profile L1 in the Changbaishan-Tianchi volcanic region. It is used in this study that seismic traveltime inversion for simultaneous determination of 2-D velocity and interface structure of the crust and upper mantle. The result shows that, under Changbaishan-Tianchi crater, there exists a low-velocity body in the shape of an inverted triangle, and the crustal reflecting boundaries and Moho all become lower by a varying margin of 2-6 km, forming a crustal root which is assumed to be the Changbaishan-Tianchi volcanic system. Finally, we make a comparison between our 2-D velocity model and the result from the studies by using trial-and-error forward modeling with SEIS83.展开更多
D-T_(2)two-dimensional nuclear magnetic resonance(2D NMR)logging technology can distinguish pore fluid types intuitively,and it is widely used in oil and gas exploration.Many 2D NMR inversion methods(e.g.,truncated si...D-T_(2)two-dimensional nuclear magnetic resonance(2D NMR)logging technology can distinguish pore fluid types intuitively,and it is widely used in oil and gas exploration.Many 2D NMR inversion methods(e.g.,truncated singular value decomposition(TSVD),Butler-Reds-Dawson(BRD),LM-norm smoothing,and TIST-L1 regularization methods)have been proposed successively,but most are limited to numerical simulations.This study focused on the applicability of different inversion methods for NMR logging data of various acquisition sequences,from which the optimal inversion method was selected based on the comparative analysis.First,the two-dimensional NMR logging principle was studied.Then,these inversion methods were studied in detail,and the precision and computational efficiency of CPMG and diffusion editing(DE)sequences obtained from oil-water and gas-water models were compared,respectively.The inversion results and calculation time of truncated singular value decomposition(TSVD),Butler-Reds-Dawson(BRD),LM-norm smoothing,and TIST-L1 regularization were compared and analyzed through numerical simulations.The inversion method was optimized to process SP mode logging data from the MR Scanner instrument.The results showed that the TIST-regularization and LM-norm smoothing methods were more accurate for the CPMG and DE sequence echo trains of the oil-water and gas-water models.However,the LM-norm smoothing method was less time-consuming,making it more suitable for logging data processing.A case study in well A25 showed that the processing results by the LM-norm smoothing method were consistent with GEOLOG software.This demonstrates that the LM-norm smoothing method is applicable in practical NMR logging processing.展开更多
The ratio between vertical and radial amplitudes of Rayleigh waves(hereafter,the Rayleigh wave ZH ratio)is an important parameter used to constrain structures beneath seismic stations.Some previous studies have explor...The ratio between vertical and radial amplitudes of Rayleigh waves(hereafter,the Rayleigh wave ZH ratio)is an important parameter used to constrain structures beneath seismic stations.Some previous studies have explored crust and upper mantle structures by joint inversion of the Rayleigh wave ZH ratio and surface wave dispersion.However,all these studies have used a 1-D depth sensitivity kernel,and this kernel may lack precision when the structure varies a great deal laterally.Here,we present a systematic investigation of the two-dimensional(2-D)Rayleigh wave ZH ratio kernel based on the adjoint-wavefield method and perform two synthetic tests using the new kernel.The 2-D ZH ratio kernel is consistent with the traditional 1-D sensitivity kernel but has an asymmetric pattern with a preferred orientation toward the source.The predominant effect caused by heterogeneity can clearly be seen from kernels calculated from models with 2-D heterogeneities,which confirms the necessity of using the new 2-D kernel in some complex regions.Inversion tests using synthetic data show that the 2-D ZH ratio kernel has the potential to resolve small anomalies as well as complex lateral structures.展开更多
Traditional two-dimensional(2D) complex resistivity forward modeling is based on Poisson's equation but spectral induced polarization(SIP) data are the coproducts of the induced polarization(IP) and the electro...Traditional two-dimensional(2D) complex resistivity forward modeling is based on Poisson's equation but spectral induced polarization(SIP) data are the coproducts of the induced polarization(IP) and the electromagnetic induction(EMI) effects.This is especially true under high frequencies,where the EMI effect can exceed the IP effect.2D inversion that only considers the IP effect reduces the reliability of the inversion data.In this paper,we derive differential equations using Maxwell's equations.With the introduction of the Cole-Cole model,we use the finite-element method to conduct2 D SIP forward modeling that considers the EMI and IP effects simultaneously.The data-space Occam method,in which different constraints to the model smoothness and parametric boundaries are introduced,is then used to simultaneously obtain the four parameters of the Cole-Cole model using multi-array electric field data.This approach not only improves the stability of the inversion but also significantly reduces the solution ambiguity.To improve the computational efficiency,message passing interface programming was used to accelerate the 2D SIP forward modeling and inversion.Synthetic datasets were tested using both serial and parallel algorithms,and the tests suggest that the proposed parallel algorithm is robust and efficient.展开更多
Low-dimensional halide perovskites have become the most promising candidates for X-ray imaging,yet the issues of the poor chemical stability of hybrid halide perovskite,the high poisonousness of lead halides and the r...Low-dimensional halide perovskites have become the most promising candidates for X-ray imaging,yet the issues of the poor chemical stability of hybrid halide perovskite,the high poisonousness of lead halides and the relatively low detectivity of the lead-free halide perovskites which seriously restrain its commercialization.Here,we developed a solution inverse temperature crystal growth(ITCG)method to bring-up high quality Cs_(3)Cu_(2)I_(5)crystals with large size of centimeter order,in which the oleic acid(OA)is introduced as an antioxidative ligand to inhibit the oxidation of cuprous ions effieiently,as well as to decelerate the crystallization rate remarkalby.Based on these fine crystals,the vapor deposition technique is empolyed to prepare high quality Cs_(3)Cu_(2)I_(5)films for efficient X-ray imaging.Smooth surface morphology,high light yields and short decay time endow the Cs_(3)Cu_(2)I_(5)films with strong radioluminescence,high resolution(12 lp/mm),low detection limits(53 nGyair/s)and desirable stability.Subsequently,the Cs_(3)Cu_(2)I_(5)films have been applied to the practical radiography which exhibit superior X-ray imaging performance.Our work provides a paradigm to fabricate nonpoisonous and chemically stable inorganic halide perovskite for X-ray imaging.展开更多
In this work,a sponge-like polysulfonamide(PSA)/SiO_2 composite membrane is unprecedentedly prepared by the phase inversion method,and successfully demonstrated as a novel separator of lithium-ion batteries(LIBs).Comp...In this work,a sponge-like polysulfonamide(PSA)/SiO_2 composite membrane is unprecedentedly prepared by the phase inversion method,and successfully demonstrated as a novel separator of lithium-ion batteries(LIBs).Compared to the commercial polypropylene(PP) separator,the sponge-like PSA/SiO_2 composite possesses better physical and electrochemical properties,such as higher porosity,ionic conductivity,thermal stability and flame retarding ability.The LiCoO_2/Li half-cells using the sponge-like composite separator demonstrate superior rate capability and cyclability over those using the commercial PP separator.Moreover,the sponge-like composite separator can ensure the normal operation of LiCoO_2/Li half-cell at an extremely high temperature of 90 °C,while the commercial PP separator cannot.All these encouraging results suggest that this phase inversion based sponge-like PSA/SiO_2 composite separator is really a promising separator for high performance LIBs.展开更多
The diagnosis of the recurrent optic neuritis is commonly established clinically,and sometimes it could be challenging because the involved optic nerve does not always show significant enhancement on conventional cont...The diagnosis of the recurrent optic neuritis is commonly established clinically,and sometimes it could be challenging because the involved optic nerve does not always show significant enhancement on conventional contrast enhanced-T1 weighted imaging(CE-T1W1).In this paper,we reported a middle-aged female with early diagnosis of recurrent optic neuritis using contrast-enhanced T2 fluid-attenuated inversion recovery imaging(CET2FLAIR).The involved optic nerve presented evident enhancement on CE-T2FLAIR and no enhancement on CE-T1W1.This case suggested that the CE-T2FLAIR may be a useful diagnostic tool specifically for the recurrent optic neuritis in clinical practice.展开更多
Nuclear Magnetic inversion is the basis of NMR Resonance (NMR) T2 logging interpretation. The regularization parameter selection of the penalty term directly influences the NMR T2 inversion result. We implemented b...Nuclear Magnetic inversion is the basis of NMR Resonance (NMR) T2 logging interpretation. The regularization parameter selection of the penalty term directly influences the NMR T2 inversion result. We implemented both norm smoothing and curvature smoothing methods for NMR T2 inversion, and compared the inversion results with respect to the optimal regular- ization parameters ((Xopt) which were selected by the dis- crepancy principle (DP), generalized cross-validation (GCV), S-curve, L-curve, and the slope of L-curve methods, respectively. The numerical results indicate that the DP method can lead to an oscillating or oversmoothed solution which is caused by an inaccurately estimated noise level. The (Xopt selected by the L-curve method is occa- sionally small or large which causes an undersmoothed or oversmoothed T2 distribution. The inversion results from GCV, S-curve and the slope of L-curve methods show satisfying inversion results. The slope of the L-curve method with less computation is more suitable for NMR T2 inversion. The inverted T2 distribution from norm smoothing is better than that from curvature smoothing when the noise level is high.展开更多
基金Projects(41074085,41374118)supported by the National Natural Science Foundation of ChinaProject(20120162110015)supported by Doctoral Fund of Ministry of Education of ChinaProject(NCET-12-0551)supported by Program for New Century Excellent Talents in University,China
文摘The 2D data processing adopted by the high-density resistivity method regards the geological structures as two degrees, which makes the results of the 2D data inversion only an approximate interpretation;the accuracy and effect can not meet the precise requirement of the inversion. Two typical models of the geological bodies were designed, and forward calculation was carried out using finite element method. The forward-modeled profiles were obtained. 1% Gaussian random error was added in the forward models and then 2D and 3D inversions using a high-density resistivity method were undertaken to realistically simulate field data and analyze the sensitivity of the 2D and 3D inversion algorithms to noise. Contrast between the 2D and 3D inversion results of least squares inversion shows that two inversion results of high-density resistivity method all can basically reflect the spatial position of an anomalous body. However, the 3D inversion can more effectively eliminate the influence of interference from Gaussian random error and better reflect the distribution of resistivity in the anomalous bodies. Overall, the 3D inversion was better than 2D inversion in terms of embodying anomalous body positions, morphology and resistivity properties.
文摘The synthesis of new 4-imino-4H-chromeno[2,3-d]pyrimidin-3(5H)-amine in four steps including one step under microwave dielectric heating is reported. The structural identity of the synthesized compounds was established according to their spectroscopic analysis, such as FT-IR, NMR and mass spectroscopy. These new compounds were tested for their antiproliferative activities on seven representative human tumoral cell lines (Huh7 D12, Caco2, MDA-MB231, MDA-MB468, HCT116, PC3 and MCF7) and also on fibroblasts. Among them, only the compounds 6c showed micromolar cytotoxic activity on tumor cell lines (1.8 50 50 > 25 μM). Finally, in silico ADMET studies ware performed to investigate the possibility of using of the identified compound 6c as potential anti-tumor compound.
文摘An engineering system approach of 2-D cylindrical model of transient mass balance calculations of ozone and other concerned chemicals along with fourteen photolysis, ozone-generating and ozone-depleting chemical reaction equations was developed, validated, and used for studying the ozone concentrations, distribution and peak of the layer, ozone depletion and total ozone abundance in the stratosphere. The calculated ozone concentrations and profile at both the Equator and a 60˚N location were found to follow closely with the measured data. The calculated average ozone concentration was within 1% of the measured average, and the deviation of ozone profiles was within 14%. The monthly evolution of stratospheric ozone concentrations and distribution above the Equator was studied with results discussed in details. The influences of slow air movement in both altitudinal and radial directions on ozone concentrations and profile in the stratosphere were explored and discussed. Parametric studies of the influences of gas diffusivities of ozone D<sub>O3</sub> and active atomic oxygen D<sub>O</sub> on ozone concentrations and distributions were also studied and delineated. Having both influences through physical diffusion and chemical reactions, the diffusivity (and diffusion) of atomic oxygen D<sub>O</sub> was found to be more sensitive and important than that of ozone D<sub>O3</sub> on ozone concentrations and distribution. The 2-D ozone model present in this paper for stratospheric ozone and its layer and depletion is shown to be robust, convenient, efficient, and executable for analyzing the complex ozone phenomena in the stratosphere. .
基金Key Project (95-11-02-01) from China Seismological Bureau.Contribution No. RCEG200129, Research Center of Exploration Geophysi
文摘2-D velocity structure and tectonics of the crust and upper mantle is revealed by inversion of seismic refraction and wide-angle reflection traveltimes acquired along the profile L1 in the Changbaishan-Tianchi volcanic region. It is used in this study that seismic traveltime inversion for simultaneous determination of 2-D velocity and interface structure of the crust and upper mantle. The result shows that, under Changbaishan-Tianchi crater, there exists a low-velocity body in the shape of an inverted triangle, and the crustal reflecting boundaries and Moho all become lower by a varying margin of 2-6 km, forming a crustal root which is assumed to be the Changbaishan-Tianchi volcanic system. Finally, we make a comparison between our 2-D velocity model and the result from the studies by using trial-and-error forward modeling with SEIS83.
基金sponsored by the National Natural Science Foundation of China(Nos.42174149,41774144)the National Major Projects(No.2016ZX05014-001).
文摘D-T_(2)two-dimensional nuclear magnetic resonance(2D NMR)logging technology can distinguish pore fluid types intuitively,and it is widely used in oil and gas exploration.Many 2D NMR inversion methods(e.g.,truncated singular value decomposition(TSVD),Butler-Reds-Dawson(BRD),LM-norm smoothing,and TIST-L1 regularization methods)have been proposed successively,but most are limited to numerical simulations.This study focused on the applicability of different inversion methods for NMR logging data of various acquisition sequences,from which the optimal inversion method was selected based on the comparative analysis.First,the two-dimensional NMR logging principle was studied.Then,these inversion methods were studied in detail,and the precision and computational efficiency of CPMG and diffusion editing(DE)sequences obtained from oil-water and gas-water models were compared,respectively.The inversion results and calculation time of truncated singular value decomposition(TSVD),Butler-Reds-Dawson(BRD),LM-norm smoothing,and TIST-L1 regularization were compared and analyzed through numerical simulations.The inversion method was optimized to process SP mode logging data from the MR Scanner instrument.The results showed that the TIST-regularization and LM-norm smoothing methods were more accurate for the CPMG and DE sequence echo trains of the oil-water and gas-water models.However,the LM-norm smoothing method was less time-consuming,making it more suitable for logging data processing.A case study in well A25 showed that the processing results by the LM-norm smoothing method were consistent with GEOLOG software.This demonstrates that the LM-norm smoothing method is applicable in practical NMR logging processing.
基金This study was funded by the National Key R&D Program of China(2016YFC0600301,2018YFC1503400)the National Natural Science Foundation of China(41790464)+1 种基金Natural Science Foundation of Jiangsu Province of China(BK20190499)the Fundamental Research Funds for the Central Universities(2019B0071428).
文摘The ratio between vertical and radial amplitudes of Rayleigh waves(hereafter,the Rayleigh wave ZH ratio)is an important parameter used to constrain structures beneath seismic stations.Some previous studies have explored crust and upper mantle structures by joint inversion of the Rayleigh wave ZH ratio and surface wave dispersion.However,all these studies have used a 1-D depth sensitivity kernel,and this kernel may lack precision when the structure varies a great deal laterally.Here,we present a systematic investigation of the two-dimensional(2-D)Rayleigh wave ZH ratio kernel based on the adjoint-wavefield method and perform two synthetic tests using the new kernel.The 2-D ZH ratio kernel is consistent with the traditional 1-D sensitivity kernel but has an asymmetric pattern with a preferred orientation toward the source.The predominant effect caused by heterogeneity can clearly be seen from kernels calculated from models with 2-D heterogeneities,which confirms the necessity of using the new 2-D kernel in some complex regions.Inversion tests using synthetic data show that the 2-D ZH ratio kernel has the potential to resolve small anomalies as well as complex lateral structures.
基金jointly sponsored by the National Natural Science Foundation of China(Grant No.41374078)the Geological Survey Projects of the Ministry of Land and Resources of China(Grant Nos.12120113086100 and 12120113101300)Beijing Higher Education Young Elite Teacher Project
文摘Traditional two-dimensional(2D) complex resistivity forward modeling is based on Poisson's equation but spectral induced polarization(SIP) data are the coproducts of the induced polarization(IP) and the electromagnetic induction(EMI) effects.This is especially true under high frequencies,where the EMI effect can exceed the IP effect.2D inversion that only considers the IP effect reduces the reliability of the inversion data.In this paper,we derive differential equations using Maxwell's equations.With the introduction of the Cole-Cole model,we use the finite-element method to conduct2 D SIP forward modeling that considers the EMI and IP effects simultaneously.The data-space Occam method,in which different constraints to the model smoothness and parametric boundaries are introduced,is then used to simultaneously obtain the four parameters of the Cole-Cole model using multi-array electric field data.This approach not only improves the stability of the inversion but also significantly reduces the solution ambiguity.To improve the computational efficiency,message passing interface programming was used to accelerate the 2D SIP forward modeling and inversion.Synthetic datasets were tested using both serial and parallel algorithms,and the tests suggest that the proposed parallel algorithm is robust and efficient.
基金the financially support of the National Natural Science Foundation of China(12164051)the Joint Foundation of Provincial Science and Technology Department-Double First-class Construction of Yunnan University(2019FY003016)+4 种基金the Young Top Talent Project of Yunnan Province(YNWR-QNBJ-2018-229)the financially support by Yunnan Major Scientific and Technological Projects(202202AG050016)Advanced Analysis and Measurement Center of Yunnan University for the sample characterization service and the Postgraduate Research and Innovation Foundation of Yunnan University(2021Y036)the financially support of the National Natural Science Foundation of China(62064013)the Application Basic Research Project of Yunnan Province[2019FB130]。
文摘Low-dimensional halide perovskites have become the most promising candidates for X-ray imaging,yet the issues of the poor chemical stability of hybrid halide perovskite,the high poisonousness of lead halides and the relatively low detectivity of the lead-free halide perovskites which seriously restrain its commercialization.Here,we developed a solution inverse temperature crystal growth(ITCG)method to bring-up high quality Cs_(3)Cu_(2)I_(5)crystals with large size of centimeter order,in which the oleic acid(OA)is introduced as an antioxidative ligand to inhibit the oxidation of cuprous ions effieiently,as well as to decelerate the crystallization rate remarkalby.Based on these fine crystals,the vapor deposition technique is empolyed to prepare high quality Cs_(3)Cu_(2)I_(5)films for efficient X-ray imaging.Smooth surface morphology,high light yields and short decay time endow the Cs_(3)Cu_(2)I_(5)films with strong radioluminescence,high resolution(12 lp/mm),low detection limits(53 nGyair/s)and desirable stability.Subsequently,the Cs_(3)Cu_(2)I_(5)films have been applied to the practical radiography which exhibit superior X-ray imaging performance.Our work provides a paradigm to fabricate nonpoisonous and chemically stable inorganic halide perovskite for X-ray imaging.
基金Supported by the funding from "135" Projects Fund of CAS-QIBEBT Director Innovation FoundationThink-Tank Mutual Fund of Qingdao Energy Storage Industry Scientific Research+3 种基金Qingdao Key Lab of Solar Energy Utilization and Energy Storage Technologythe Strategic Priority Research Program of the Chinese Academy of Sciences(XDA09010105)National Natural Science Foundation of China(51502319)Shandong Provincial Natural Science Foundation(ZR2016BQ18)
文摘In this work,a sponge-like polysulfonamide(PSA)/SiO_2 composite membrane is unprecedentedly prepared by the phase inversion method,and successfully demonstrated as a novel separator of lithium-ion batteries(LIBs).Compared to the commercial polypropylene(PP) separator,the sponge-like PSA/SiO_2 composite possesses better physical and electrochemical properties,such as higher porosity,ionic conductivity,thermal stability and flame retarding ability.The LiCoO_2/Li half-cells using the sponge-like composite separator demonstrate superior rate capability and cyclability over those using the commercial PP separator.Moreover,the sponge-like composite separator can ensure the normal operation of LiCoO_2/Li half-cell at an extremely high temperature of 90 °C,while the commercial PP separator cannot.All these encouraging results suggest that this phase inversion based sponge-like PSA/SiO_2 composite separator is really a promising separator for high performance LIBs.
文摘The diagnosis of the recurrent optic neuritis is commonly established clinically,and sometimes it could be challenging because the involved optic nerve does not always show significant enhancement on conventional contrast enhanced-T1 weighted imaging(CE-T1W1).In this paper,we reported a middle-aged female with early diagnosis of recurrent optic neuritis using contrast-enhanced T2 fluid-attenuated inversion recovery imaging(CET2FLAIR).The involved optic nerve presented evident enhancement on CE-T2FLAIR and no enhancement on CE-T1W1.This case suggested that the CE-T2FLAIR may be a useful diagnostic tool specifically for the recurrent optic neuritis in clinical practice.
基金funded by Shell International Exploration and Production Inc.(PT45371)the National Natural Science Foundation of China-China National Petroleum Corporation Petrochemical Engineering United Fund(U1262114)the National Natural Science Foundation of China(41272163)
文摘Nuclear Magnetic inversion is the basis of NMR Resonance (NMR) T2 logging interpretation. The regularization parameter selection of the penalty term directly influences the NMR T2 inversion result. We implemented both norm smoothing and curvature smoothing methods for NMR T2 inversion, and compared the inversion results with respect to the optimal regular- ization parameters ((Xopt) which were selected by the dis- crepancy principle (DP), generalized cross-validation (GCV), S-curve, L-curve, and the slope of L-curve methods, respectively. The numerical results indicate that the DP method can lead to an oscillating or oversmoothed solution which is caused by an inaccurately estimated noise level. The (Xopt selected by the L-curve method is occa- sionally small or large which causes an undersmoothed or oversmoothed T2 distribution. The inversion results from GCV, S-curve and the slope of L-curve methods show satisfying inversion results. The slope of the L-curve method with less computation is more suitable for NMR T2 inversion. The inverted T2 distribution from norm smoothing is better than that from curvature smoothing when the noise level is high.