Local and global optimization methods are widely used in geophysical inversion but each has its own advantages and disadvantages. The combination of the two methods will make it possible to overcome their weaknesses. ...Local and global optimization methods are widely used in geophysical inversion but each has its own advantages and disadvantages. The combination of the two methods will make it possible to overcome their weaknesses. Based on the simulated annealing genetic algorithm （SAGA） and the simplex algorithm, an efficient and robust 2-D nonlinear method for seismic travel-time inversion is presented in this paper. First we do a global search over a large range by SAGA and then do a rapid local search using the simplex method. A multi-scale tomography method is adopted in order to reduce non-uniqueness. The velocity field is divided into different spatial scales and velocities at the grid nodes are taken as unknown parameters. The model is parameterized by a bi-cubic spline function. The finite-difference method is used to solve the forward problem while the hybrid method combining multi-scale SAGA and simplex algorithms is applied to the inverse problem. The algorithm has been applied to a numerical test and a travel-time perturbation test using an anomalous low-velocity body. For a practical example, it is used in the study of upper crustal velocity structure of the A＇nyemaqen suture zone at the north-east edge of the Qinghai-Tibet Plateau. The model test and practical application both prove that the method is effective and robust.展开更多
The crustal structure of Xiachayu-Gonghe geophysical profile in eastern Tibetan plateau is simulated with Bouguer anomaly corrected for sediments and lithosphere. The forward simulation shows that the thickness of upp...The crustal structure of Xiachayu-Gonghe geophysical profile in eastern Tibetan plateau is simulated with Bouguer anomaly corrected for sediments and lithosphere. The forward simulation shows that the thickness of upper crust in eastem Tibetan plateau is about 20 km, and the density is 2.78 × 10^3 kg/m^3. The bottom interface of middle crust changes from 30 km to 40 km, the density of middle crust is 2.89 × 10^3 kg/m^3. The materials with low density of 2.78 × 10^3 kg/m^3 exist in middle crust, and those with high density of 3.33 × 10^3 kg/m^3 exist at the bottom of middle crust between Wenquan and Tanggemu. The density is 3.10× 10^3 kg/m^3 in lower crust. The shallowest depth of Moho interface is about 56 km, and the deepest one is about 74 km, the undulation of interface is large, the deep Moho is located in Xiachayu, Chayu, Nujiang, and Wenquan. The crustal density of eastern Tibetan plateau is larger than that of central section; the low velocity layers are located in middle crust and bottom in eastern Tibetan plateau and at the bottom of the upper crust in the central plateau.展开更多
基金supported by the National Natural Science Foundation of China (Grant Nos.40334040 and 40974033)the Promoting Foundation for Advanced Persons of Talent of NCWU
文摘Local and global optimization methods are widely used in geophysical inversion but each has its own advantages and disadvantages. The combination of the two methods will make it possible to overcome their weaknesses. Based on the simulated annealing genetic algorithm （SAGA） and the simplex algorithm, an efficient and robust 2-D nonlinear method for seismic travel-time inversion is presented in this paper. First we do a global search over a large range by SAGA and then do a rapid local search using the simplex method. A multi-scale tomography method is adopted in order to reduce non-uniqueness. The velocity field is divided into different spatial scales and velocities at the grid nodes are taken as unknown parameters. The model is parameterized by a bi-cubic spline function. The finite-difference method is used to solve the forward problem while the hybrid method combining multi-scale SAGA and simplex algorithms is applied to the inverse problem. The algorithm has been applied to a numerical test and a travel-time perturbation test using an anomalous low-velocity body. For a practical example, it is used in the study of upper crustal velocity structure of the A＇nyemaqen suture zone at the north-east edge of the Qinghai-Tibet Plateau. The model test and practical application both prove that the method is effective and robust.
基金Supported by the National Natural Science Foundation of China (No. 40874036,No. 90814009,No. 40774011)the National Science & Technology Major Project of China (No. 2008ZX05008-006-062)
文摘The crustal structure of Xiachayu-Gonghe geophysical profile in eastern Tibetan plateau is simulated with Bouguer anomaly corrected for sediments and lithosphere. The forward simulation shows that the thickness of upper crust in eastem Tibetan plateau is about 20 km, and the density is 2.78 × 10^3 kg/m^3. The bottom interface of middle crust changes from 30 km to 40 km, the density of middle crust is 2.89 × 10^3 kg/m^3. The materials with low density of 2.78 × 10^3 kg/m^3 exist in middle crust, and those with high density of 3.33 × 10^3 kg/m^3 exist at the bottom of middle crust between Wenquan and Tanggemu. The density is 3.10× 10^3 kg/m^3 in lower crust. The shallowest depth of Moho interface is about 56 km, and the deepest one is about 74 km, the undulation of interface is large, the deep Moho is located in Xiachayu, Chayu, Nujiang, and Wenquan. The crustal density of eastern Tibetan plateau is larger than that of central section; the low velocity layers are located in middle crust and bottom in eastern Tibetan plateau and at the bottom of the upper crust in the central plateau.
