This paper focuses on estimating a new high-resolution Earth’s gravity field model named SGG-UGM-2 from satellite gravimetry,satellite altimetry,and Earth Gravitational Model 2008(EGM2008)-derived gravity data based ...This paper focuses on estimating a new high-resolution Earth’s gravity field model named SGG-UGM-2 from satellite gravimetry,satellite altimetry,and Earth Gravitational Model 2008(EGM2008)-derived gravity data based on the theory of the ellipsoidal harmonic analysis and coefficient transformation(EHA-CT).We first derive the related formulas of the EHA-CT method,which is used for computing the spherical harmonic coefficients from grid area-mean and point gravity anomalies on the ellipsoid.The derived formulas are successfully evaluated based on numerical experiments.Then,based on the derived least-squares formulas of the EHA-CT method,we develop the new model SGG-UGM-2 up to degree 2190 and order 2159 by combining the observations of the Gravity Field and Steady-State Ocean Circulation Explorer(GOCE),the normal equation of the Gravity Recovery and Climate Experiment(GRACE),marine gravity data derived from satellite altimetry data,and EGM2008-derived continental gravity data.The coefficients of degrees 251–2159 are estimated by solving the block-diagonal form normal equations of surface gravity anomalies(including the marine gravity data).The coefficients of degrees 2–250 are determined by combining the normal equations of satellite observations and surface gravity anomalies.The variance component estimation technique is used to estimate the relative weights of different observations.Finally,global positioning system(GPS)/leveling data in the mainland of China and the United States are used to validate SGG-UGM-2 together with other models,such as European improved gravity model of the earth by new techniques(EIGEN)-6C4,GECO,EGM2008,and SGG-UGM-1(the predecessor of SGG-UGM-2).Compared to other models,the model SGG-UGM-2 shows a promising performance in the GPS/leveling validation.All GOCE-related models have similar performances both in the mainland of China and the United States,and better performances than that of EGM2008 in the mainland of China.Due to the contribution of GRACE data and the new marine gravity anomalies,SGG-UGM-2 is slightly better than SGG-UGM-1 both in the mainland of China and the United States.展开更多
To create a new prediction model, the unbiased GM (1,1) model is optimized by the five-point slide method in this paper. Then, based on the occurrence areas of dce blast in Enshi District during 1995 -2004, the new ...To create a new prediction model, the unbiased GM (1,1) model is optimized by the five-point slide method in this paper. Then, based on the occurrence areas of dce blast in Enshi District during 1995 -2004, the new model and unbiased GM (1, 1 ) model are applied to predict the occurrence areas of rice blast during 2005 -2010. Predicting outcomes show that the prediction accuracy of five-point unbiased sliding optimized GM (1, 1 ) model is higher than the unbiased GM (1,1) model. Finally, combined with the prediction results, the author provides some suggestion for Enshi District in the prevention and control of rice blast in 2010.展开更多
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.展开更多
Using a modified C D function and stochastic frontier model, the paper analyzed China's cotton yield capacity and found that the yield and technical efficiency of China's cotton planting system can be increas...Using a modified C D function and stochastic frontier model, the paper analyzed China's cotton yield capacity and found that the yield and technical efficiency of China's cotton planting system can be increased by the use of genetically modified (GM) varieties.展开更多
Firstly, the research progress of grey model GM (1,1) is summarized, which is divided into three development stages: assimilation, alienation and melting stages. Then, the matrix analysis theory is used to study th...Firstly, the research progress of grey model GM (1,1) is summarized, which is divided into three development stages: assimilation, alienation and melting stages. Then, the matrix analysis theory is used to study the modeling mechanism of GM (1,1), which decomposes the modeling data matrix into raw data transformation matrix, accumulated generating operation matrix and background value selection matrix. The changes of these three matrices are the essential reasons affecting the modeling and the accuracy of GM (1,1). Finally, the paper proposes a generalization grey model GGM (1,1), which is a extended form of GM (1,1) and also a unified form of model GM (1,1), model GM (1,1,α), stage grey model, hopping grey model, generalized accumulated model, strengthening operator model, weakening operator model and unequal interval model. And the theory and practical significance of the extended model is analyzed.展开更多
In gravity-anomaly-based prospecting, the computational and memory requirements for practical numerical modeling are potentially enormous. Achieving an efficient and precise inversion for gravity anomaly imaging over ...In gravity-anomaly-based prospecting, the computational and memory requirements for practical numerical modeling are potentially enormous. Achieving an efficient and precise inversion for gravity anomaly imaging over large-scale and complex terrain requires additional methods. To this end, we have proposed a new topography-capable By performing a two-dimensional Fourier transform in the horizontal directions, threedimensional partial differential equations in the spatial domain were transformed into a group of independent, one-dimensional differential equations engaged with different wave numbers. These independent differential equations are highly parallel across different wave numbers. differential equations with different wave numbers, and the efficiency of solving fixedbandwidth linear equations was further improved by a chasing method. In a synthetic test, a prism model was used to verify the accuracy and reliability of the proposed algorithm by comparing the numerical solution with the analytical solution. We studied the computational precision and efficiency with and without topography using different Fourier transform methods. The results showed that the Guass-FFT method has higher numerical precision, while the standard FFT method is superior, in terms of computation time, for inversion and quantitative interpretation under complicated terrain.展开更多
Based on the optimization method, a new modified GM (1,1) model is presented, which is characterized by more accuracy prediction for the grey modeling.
Theoretical analysis and practical observations show that fault dislocations can change the gravity field around the fault. Gravity changes which were caused by the repeated dislocations over a long period of time wer...Theoretical analysis and practical observations show that fault dislocations can change the gravity field around the fault. Gravity changes which were caused by the repeated dislocations over a long period of time were superimposed on the Bougeur gravity anomalies. These anomalies became the evidence of historical movement of fault as well as provide a way for the study of paleo earthquakes. This paper investigates inversion methods for the geological dislocation modeling of faults using the local Bouguer's gravity anomalies. To remove the effects of the irrelevant part of gravity anomalies to fault movements, we propose the robust nonlinear inversion method and set up the corresponding algorithm. Modeling examples indicate that the Marquardt's and Baye's least squares solutions depart from the true solution due to the attraction of gross errors in the data. The more seriously the data is contaminated, the more seriously the solutions are biased. In contrast, the proposed robust Marquardt's and Baye's inversion solutions can still maintain consistency with the solution without gross errors, even though 50 percent of the data is contaminated. This indicates that the proposed robust methods are effective. Using the proposed methods, we invert the geological dislocation models of the faults around the Erhai Lake in West Yunnan. The results show that the Northern Cangdong fault and the Erhai fault are normal dip slip faults with about 4 to 5 km dislocations; and that the Southern Cangdong fault has a less dip slip compared with the former two. A satisfactory fitting between the theoretical values of the inversion solution and the actual local gravity field is achievable.展开更多
The effects of rotation and gravity on an electro-magneto-thermoelastic medium with diffusion and voids in a generalized thermoplastic half-space are studied by using the Lord-Shulman (L-S) model and the dual-phase-la...The effects of rotation and gravity on an electro-magneto-thermoelastic medium with diffusion and voids in a generalized thermoplastic half-space are studied by using the Lord-Shulman (L-S) model and the dual-phase-lag (DPL) model. The analytical solutions for the displacements, stresses, temperature, diffusion concentration, and volume fraction field with different values of the magnetic field, the rotation, the gravity, and the initial stress are obtained and portrayed graphically. The results indicate that the effects of gravity, rotation, voids, diffusion, initial stress, and electromagnetic field are very pronounced on the physical properties of the material.展开更多
As a kind of mathematical model, grey systems predi ct ion model has been widely applied to economy, management and engineering technol ogy. In 1982, Professor Deng Ju-long presented GM prediction model. Then some o t...As a kind of mathematical model, grey systems predi ct ion model has been widely applied to economy, management and engineering technol ogy. In 1982, Professor Deng Ju-long presented GM prediction model. Then some o ther scholars made improvements on GM model. Of course, much still should be don e to develop it. What the scholars have done is to take the first component of X (1) as the starting conditions of the grey differential model. It occ urs that the new information can not be used enough. This paper is addressed to choose the nth component of X (1) as the starting conditions to improv e the models. The main results of the paper is given in Theorem 2: The time response function of the grey differential equation x (0)(k)+az (1)(k)=b is given by x (1)(k)=x (1)(n)-ba e -a(k-n )+ba. and Theorem4: The time response of the grey Verhulst model is given by (1)(k) =ax (1)(n)bx (1)(n)+(a-bx (1)(n))ae a(k-n). As the new information is fully used, the accuracy of prediction is improved gre atly. Therefore, the new model with a certain theoretical and practical value.展开更多
Regional surface gravity data and global satellite magnetic data have been utilized to generate a preliminary model of the crustal structure along a southwest-northeast profile (Gadra-Fatehpur) through western Rajas...Regional surface gravity data and global satellite magnetic data have been utilized to generate a preliminary model of the crustal structure along a southwest-northeast profile (Gadra-Fatehpur) through western Rajasthan.The study area represents the western part of the Indian continental landmass which has undergone several major episodes of repeated subduction/collision,plume traces and rifting from Archaean to recent times.The temporal and spatial relationship between the various geotectonic provinces is quite complex,thereby limiting the emergence of a suitable crustal structure model for this region.Exposures of the Malani Igneous Suite (MIS),a product of bimodal volcanism (~780 Ma),and considered to be the third largest felsic magmatic province of the world,is evident along the profile and also to the southwest of the study area.The easternmost part of the profile is close to the DAFB (Delhi Aravalli Fold Belt),a Proterozoic orogenic belt.This study probes the geometry of the different crustal units in terms of density and susceptibility variations in order to decipher the imprints of the major tectonic processes the region has undergone.In order to decipher the crustal geometry of the Gadra-Fatehpur profile,two NW-SE gravity and magnetic profile vertical sections (A-A' in the south and B-B' in the north) are modelled on the basis of the constraints provided from previous seismic models.The crustal model of the Gadra-Fatehpur profile is composed of alluvium,Tertiary sediments,MIS,Marwar Supergroup,low-density layers (LDLs) and the middle-lower crustal layers,with a distinct change in configuration from the southwest to northeast.The Moho dips from SW to NE,the MIS in the SW gives way to the thick pile of the Marwar Supergroup to the NE.The evolution of MIS has been suggested to have occurred as a consequence of delamination of the upper mantle.LDLs are incorporated in Gadra-Fatehpur model.In the SW,LDL (2550 kg/m3) lies below the MIS in the NE,another LDL (2604 kg/m3) is depicted below the mid-crustal layer.展开更多
The spatial interaction model is an effective way to explore the geographical disparities inherent in the Belt and Road Initiative(BRI) by simulating spatial flows. The traditional gravity model implies the hypothesis...The spatial interaction model is an effective way to explore the geographical disparities inherent in the Belt and Road Initiative(BRI) by simulating spatial flows. The traditional gravity model implies the hypothesis of equilibrium points without any reference to when or how to achieve it. In this paper, a dynamic gravity model was established based on the Maximum Entropy(MaxEnt) theory to estimate and monitor the interconnection intensity and dynamic characters of bilateral relations. In order to detect the determinants of interconnection intensity, a Geodetector method was applied to identify and evaluate the determinants of spatial networks in five dimensions. The empirical study clearly demonstrates a heterogeneous and non-circular spatial structure. The main driving forces of spatial-temporal evolution are foreign direct investment, tourism and railway infrastructure construction, while determinants in different sub-regions show obvious spatial differentiation. Southeast Asian countries are typically multi-island area where aviation infrastructure plays a more important role. North and Central Asian countries regard oil as a pillar industry where power and port facilities have a greater impact on the interconnection. While Western Asian countries are mostly influenced by the railway infrastructure, Eastern European countries already have relatively robust infrastructure where tariff policies provide a greater impetus.展开更多
For the classical GM(1,1)model,the prediction accuracy is not high,and the optimization of the initial and background values is one-sided.In this paper,the Lagrange mean value theorem is used to construct the backgrou...For the classical GM(1,1)model,the prediction accuracy is not high,and the optimization of the initial and background values is one-sided.In this paper,the Lagrange mean value theorem is used to construct the background value as a variable related to k.At the same time,the initial value is set as a variable,and the corresponding optimal parameter and the time response formula are determined according to the minimum value of mean relative error(MRE).Combined with the domestic natural gas annual consumption data,the classical model and the improved GM(1,1)model are applied to the calculation and error comparison respectively.It proves that the improved model is better than any other models.展开更多
基金We appreciate the help from Torsten Mayer-Gürr and Andreas Kvas for providing us the NEQ system of the ITSG-Grace2018 model.This research was financially supported by the National Natural Science Foundation of China(41574019 and 41774020)the German Academic Exchange Service(DAAD)Thematic Network Project(57421148)+2 种基金the Major Project of High-Resolution Earth Observation System,and Science Fund for Creative Research Groups of the National Natural Science Foundation of China(41721003)the Fundamental Research Funds for the Central Universities(N170103009)We also thank the editor and the anonymous reviewers for their constructive remarks that helped us to improve the quality of the manuscript.
文摘This paper focuses on estimating a new high-resolution Earth’s gravity field model named SGG-UGM-2 from satellite gravimetry,satellite altimetry,and Earth Gravitational Model 2008(EGM2008)-derived gravity data based on the theory of the ellipsoidal harmonic analysis and coefficient transformation(EHA-CT).We first derive the related formulas of the EHA-CT method,which is used for computing the spherical harmonic coefficients from grid area-mean and point gravity anomalies on the ellipsoid.The derived formulas are successfully evaluated based on numerical experiments.Then,based on the derived least-squares formulas of the EHA-CT method,we develop the new model SGG-UGM-2 up to degree 2190 and order 2159 by combining the observations of the Gravity Field and Steady-State Ocean Circulation Explorer(GOCE),the normal equation of the Gravity Recovery and Climate Experiment(GRACE),marine gravity data derived from satellite altimetry data,and EGM2008-derived continental gravity data.The coefficients of degrees 251–2159 are estimated by solving the block-diagonal form normal equations of surface gravity anomalies(including the marine gravity data).The coefficients of degrees 2–250 are determined by combining the normal equations of satellite observations and surface gravity anomalies.The variance component estimation technique is used to estimate the relative weights of different observations.Finally,global positioning system(GPS)/leveling data in the mainland of China and the United States are used to validate SGG-UGM-2 together with other models,such as European improved gravity model of the earth by new techniques(EIGEN)-6C4,GECO,EGM2008,and SGG-UGM-1(the predecessor of SGG-UGM-2).Compared to other models,the model SGG-UGM-2 shows a promising performance in the GPS/leveling validation.All GOCE-related models have similar performances both in the mainland of China and the United States,and better performances than that of EGM2008 in the mainland of China.Due to the contribution of GRACE data and the new marine gravity anomalies,SGG-UGM-2 is slightly better than SGG-UGM-1 both in the mainland of China and the United States.
基金Supported by Science Research Project of Department of Education of Hubei Province (B20092901)~~
文摘To create a new prediction model, the unbiased GM (1,1) model is optimized by the five-point slide method in this paper. Then, based on the occurrence areas of dce blast in Enshi District during 1995 -2004, the new model and unbiased GM (1, 1 ) model are applied to predict the occurrence areas of rice blast during 2005 -2010. Predicting outcomes show that the prediction accuracy of five-point unbiased sliding optimized GM (1, 1 ) model is higher than the unbiased GM (1,1) model. Finally, combined with the prediction results, the author provides some suggestion for Enshi District in the prevention and control of rice blast in 2010.
基金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.
文摘Using a modified C D function and stochastic frontier model, the paper analyzed China's cotton yield capacity and found that the yield and technical efficiency of China's cotton planting system can be increased by the use of genetically modified (GM) varieties.
基金supported by the National Natural Science Foundation of China(70971103)the Specialized Research Fund for the Doctora Program of Higher Education(20120143110001)
文摘Firstly, the research progress of grey model GM (1,1) is summarized, which is divided into three development stages: assimilation, alienation and melting stages. Then, the matrix analysis theory is used to study the modeling mechanism of GM (1,1), which decomposes the modeling data matrix into raw data transformation matrix, accumulated generating operation matrix and background value selection matrix. The changes of these three matrices are the essential reasons affecting the modeling and the accuracy of GM (1,1). Finally, the paper proposes a generalization grey model GGM (1,1), which is a extended form of GM (1,1) and also a unified form of model GM (1,1), model GM (1,1,α), stage grey model, hopping grey model, generalized accumulated model, strengthening operator model, weakening operator model and unequal interval model. And the theory and practical significance of the extended model is analyzed.
基金supported by the Natural Science Foundation of China(No.41574127)the China Postdoctoral Science Foundation(No.2017M622608)the project for the independent exploration of graduate students at Central South University(No.2017zzts008)
文摘In gravity-anomaly-based prospecting, the computational and memory requirements for practical numerical modeling are potentially enormous. Achieving an efficient and precise inversion for gravity anomaly imaging over large-scale and complex terrain requires additional methods. To this end, we have proposed a new topography-capable By performing a two-dimensional Fourier transform in the horizontal directions, threedimensional partial differential equations in the spatial domain were transformed into a group of independent, one-dimensional differential equations engaged with different wave numbers. These independent differential equations are highly parallel across different wave numbers. differential equations with different wave numbers, and the efficiency of solving fixedbandwidth linear equations was further improved by a chasing method. In a synthetic test, a prism model was used to verify the accuracy and reliability of the proposed algorithm by comparing the numerical solution with the analytical solution. We studied the computational precision and efficiency with and without topography using different Fourier transform methods. The results showed that the Guass-FFT method has higher numerical precision, while the standard FFT method is superior, in terms of computation time, for inversion and quantitative interpretation under complicated terrain.
文摘Based on the optimization method, a new modified GM (1,1) model is presented, which is characterized by more accuracy prediction for the grey modeling.
文摘Theoretical analysis and practical observations show that fault dislocations can change the gravity field around the fault. Gravity changes which were caused by the repeated dislocations over a long period of time were superimposed on the Bougeur gravity anomalies. These anomalies became the evidence of historical movement of fault as well as provide a way for the study of paleo earthquakes. This paper investigates inversion methods for the geological dislocation modeling of faults using the local Bouguer's gravity anomalies. To remove the effects of the irrelevant part of gravity anomalies to fault movements, we propose the robust nonlinear inversion method and set up the corresponding algorithm. Modeling examples indicate that the Marquardt's and Baye's least squares solutions depart from the true solution due to the attraction of gross errors in the data. The more seriously the data is contaminated, the more seriously the solutions are biased. In contrast, the proposed robust Marquardt's and Baye's inversion solutions can still maintain consistency with the solution without gross errors, even though 50 percent of the data is contaminated. This indicates that the proposed robust methods are effective. Using the proposed methods, we invert the geological dislocation models of the faults around the Erhai Lake in West Yunnan. The results show that the Northern Cangdong fault and the Erhai fault are normal dip slip faults with about 4 to 5 km dislocations; and that the Southern Cangdong fault has a less dip slip compared with the former two. A satisfactory fitting between the theoretical values of the inversion solution and the actual local gravity field is achievable.
文摘The effects of rotation and gravity on an electro-magneto-thermoelastic medium with diffusion and voids in a generalized thermoplastic half-space are studied by using the Lord-Shulman (L-S) model and the dual-phase-lag (DPL) model. The analytical solutions for the displacements, stresses, temperature, diffusion concentration, and volume fraction field with different values of the magnetic field, the rotation, the gravity, and the initial stress are obtained and portrayed graphically. The results indicate that the effects of gravity, rotation, voids, diffusion, initial stress, and electromagnetic field are very pronounced on the physical properties of the material.
文摘As a kind of mathematical model, grey systems predi ct ion model has been widely applied to economy, management and engineering technol ogy. In 1982, Professor Deng Ju-long presented GM prediction model. Then some o ther scholars made improvements on GM model. Of course, much still should be don e to develop it. What the scholars have done is to take the first component of X (1) as the starting conditions of the grey differential model. It occ urs that the new information can not be used enough. This paper is addressed to choose the nth component of X (1) as the starting conditions to improv e the models. The main results of the paper is given in Theorem 2: The time response function of the grey differential equation x (0)(k)+az (1)(k)=b is given by x (1)(k)=x (1)(n)-ba e -a(k-n )+ba. and Theorem4: The time response of the grey Verhulst model is given by (1)(k) =ax (1)(n)bx (1)(n)+(a-bx (1)(n))ae a(k-n). As the new information is fully used, the accuracy of prediction is improved gre atly. Therefore, the new model with a certain theoretical and practical value.
文摘Regional surface gravity data and global satellite magnetic data have been utilized to generate a preliminary model of the crustal structure along a southwest-northeast profile (Gadra-Fatehpur) through western Rajasthan.The study area represents the western part of the Indian continental landmass which has undergone several major episodes of repeated subduction/collision,plume traces and rifting from Archaean to recent times.The temporal and spatial relationship between the various geotectonic provinces is quite complex,thereby limiting the emergence of a suitable crustal structure model for this region.Exposures of the Malani Igneous Suite (MIS),a product of bimodal volcanism (~780 Ma),and considered to be the third largest felsic magmatic province of the world,is evident along the profile and also to the southwest of the study area.The easternmost part of the profile is close to the DAFB (Delhi Aravalli Fold Belt),a Proterozoic orogenic belt.This study probes the geometry of the different crustal units in terms of density and susceptibility variations in order to decipher the imprints of the major tectonic processes the region has undergone.In order to decipher the crustal geometry of the Gadra-Fatehpur profile,two NW-SE gravity and magnetic profile vertical sections (A-A' in the south and B-B' in the north) are modelled on the basis of the constraints provided from previous seismic models.The crustal model of the Gadra-Fatehpur profile is composed of alluvium,Tertiary sediments,MIS,Marwar Supergroup,low-density layers (LDLs) and the middle-lower crustal layers,with a distinct change in configuration from the southwest to northeast.The Moho dips from SW to NE,the MIS in the SW gives way to the thick pile of the Marwar Supergroup to the NE.The evolution of MIS has been suggested to have occurred as a consequence of delamination of the upper mantle.LDLs are incorporated in Gadra-Fatehpur model.In the SW,LDL (2550 kg/m3) lies below the MIS in the NE,another LDL (2604 kg/m3) is depicted below the mid-crustal layer.
基金the auspices of A Category of Strategic Priority Research Program of Chinese Academy of Sciences(No.XDA20010101)。
文摘The spatial interaction model is an effective way to explore the geographical disparities inherent in the Belt and Road Initiative(BRI) by simulating spatial flows. The traditional gravity model implies the hypothesis of equilibrium points without any reference to when or how to achieve it. In this paper, a dynamic gravity model was established based on the Maximum Entropy(MaxEnt) theory to estimate and monitor the interconnection intensity and dynamic characters of bilateral relations. In order to detect the determinants of interconnection intensity, a Geodetector method was applied to identify and evaluate the determinants of spatial networks in five dimensions. The empirical study clearly demonstrates a heterogeneous and non-circular spatial structure. The main driving forces of spatial-temporal evolution are foreign direct investment, tourism and railway infrastructure construction, while determinants in different sub-regions show obvious spatial differentiation. Southeast Asian countries are typically multi-island area where aviation infrastructure plays a more important role. North and Central Asian countries regard oil as a pillar industry where power and port facilities have a greater impact on the interconnection. While Western Asian countries are mostly influenced by the railway infrastructure, Eastern European countries already have relatively robust infrastructure where tariff policies provide a greater impetus.
基金supported by the National Natural Science Foundation of China (71871106)the Blue and Green Project in Jiangsu Provincethe Six Talent Peaks Project in Jiangsu Province (2016-JY-011)
文摘For the classical GM(1,1)model,the prediction accuracy is not high,and the optimization of the initial and background values is one-sided.In this paper,the Lagrange mean value theorem is used to construct the background value as a variable related to k.At the same time,the initial value is set as a variable,and the corresponding optimal parameter and the time response formula are determined according to the minimum value of mean relative error(MRE).Combined with the domestic natural gas annual consumption data,the classical model and the improved GM(1,1)model are applied to the calculation and error comparison respectively.It proves that the improved model is better than any other models.