A method of three dimensional (3-D) model parameterization is presented that makes forward and inverse problems become easy. The velocity and interface structure of crust and upper mantle are described by a set of hig...A method of three dimensional (3-D) model parameterization is presented that makes forward and inverse problems become easy. The velocity and interface structure of crust and upper mantle are described by a set of highly smoothed functions. Shooting ray tracing method is chosen to calculate the ray paths for both forward and inverse problems. The partial derivatives of traveltime with respect to parameters of the model grids are calculated analytically while rays are being traced. Because velocity and interface functions have second-order continuous partial derivatives, the geometrical shadow zones at the surface caused by scattering and focusing of ray paths can be prevented. After ray tracing, an equation consisting of matrix and vectors for inverse problem is obtained. We use singular value decomposition method with damped factor to solve the equation. A synthetic data set which consists of several in-line profiles is used to test the methods. The results show that the methods are robust. Compared with the two dimensional method, the 3-D inversion method can give the right position of interfaces and the velocity structure when the crustal model is complicated.展开更多
3-D velocity structure of P wave in the upper mantle beneath southwestern China and its adjacent areas (10°N [similar to] 36°N, 70°E [similar to] 110°E) down to the depth of 400 km has been studied...3-D velocity structure of P wave in the upper mantle beneath southwestern China and its adjacent areas (10°N [similar to] 36°N, 70°E [similar to] 110°E) down to the depth of 400 km has been studied by using 80 974 P-wave first arrival times recorded at 165 stations from 7 053 events both within the studying areas, selected from the ISC bulletin and the Bulletin of China and NEIC fundamental seismic network. With a resolution of grid spacing of 2°×2°, the velocity heterogeneity on the horizontal profile is obvious though it attenuates with the depth increasing. On the vertical profiles of velocity along the latitude of 16°N and 24°N, the collision and extrusion of India plate to Eurasia plate is displayed, and a remarkable velocity difference between India plate and Eurasia plate is shown. In the vertical profile along the longitude of 90°E, the subducting of India plate northward beneath Eurasia plate (Tibet plateau) is also obvious. On the horizontal profile at the depth of 90 km, a slow velocity stripe from Myitkyina, Myanmar to Donghai, Vietnam seems to be related to Honghe fault belt. An illustration method of describing the resolution more directly and exactly has been proposed and utilized in this paper.展开更多
A simplified numerical model of groundwater and solute transport is developed. At large scale area there exists a big spatial scale difference between horizontal and vertical length scales. In the resultant model, the...A simplified numerical model of groundwater and solute transport is developed. At large scale area there exists a big spatial scale difference between horizontal and vertical length scales. In the resultant model, the seepage region is particularly divided into several virtual layers along the z direction and vertical 1-D columns covering x-y 2-D area according to stratum properties. The numerical algorithm is replacing the full 3-D water and mass balance analysis as the 2-D Galerkin finite element method in x- and y-directions and 1-D finite differential approach in the z direction. The reasonable method of giving minimum thickness is successfully used to handle transient change of water table, drying cells and problem of rewetting. The solution of the simplified model is preconditioned conjugate gradient and ORTHOMIN method. The validity of the developed 3-D groundwater model is tested with the typical pumping and backwater scenarios. Results of water balance of the computed example reveal the model computation reliability. Based on a representative 3-D pollution case, the solute transport module is tested against computing results using the MT3DMS. The capability and high efficiency to predict non-stationary situations of free groundwater surface and solute plume in regional scale problem is quantitatively investigated. It is shown that the proposed model is computationally effective.展开更多
By processing the CSND Rayleigh wave data with the matched filter FTAN technique, Rayleigh wave dispersion for southeast China is obtained. The 4°×4°S wave dispersion of the pure path is calculated usin...By processing the CSND Rayleigh wave data with the matched filter FTAN technique, Rayleigh wave dispersion for southeast China is obtained. The 4°×4°S wave dispersion of the pure path is calculated using random inversion scheme, and 3-D S wave velocity structure is set up. Incorporating the above-mentioned results with wide angle seismic sounding data, we studied structure framework and the extending of faults in this area, which demonstrates that the depth of Moho in South China varies from 30 to 40 km, shallower from west to east. The depth of Moho varies from 25 to 28 km for the offshore. The depth of the asthenosphere in upper mantle varies from 60 to 100 km. The depth difference of layers at the two sides of Tanlu fault is more than 10 km at the south part of the Yangtze River, and the fault extends downward more than 170 km. The fault exceeds the main land at Hainan Island and slips into the southern China Sea. Both Tanlu fault and the huge bend of gravity gradient anomaly are influenced by展开更多
On the basis of synthetic studies of geology, geophysics and geochemistry, the present Qinling Orogenic Belt can be described as a 3-D "flyover-type" geometric model with rheological layering structures. Fur...On the basis of synthetic studies of geology, geophysics and geochemistry, the present Qinling Orogenic Belt can be described as a 3-D "flyover-type" geometric model with rheological layering structures. Furthermore, the tectonic dynamic analyses have been done based on the structural geometry and kinematic features. Thus it can be concluded that its present structure has resulted from an adjustment of deep-seated mantle dynamics and lithosphere coupling relationship since the Mesozoic-Cenozoic time.展开更多
文摘A method of three dimensional (3-D) model parameterization is presented that makes forward and inverse problems become easy. The velocity and interface structure of crust and upper mantle are described by a set of highly smoothed functions. Shooting ray tracing method is chosen to calculate the ray paths for both forward and inverse problems. The partial derivatives of traveltime with respect to parameters of the model grids are calculated analytically while rays are being traced. Because velocity and interface functions have second-order continuous partial derivatives, the geometrical shadow zones at the surface caused by scattering and focusing of ray paths can be prevented. After ray tracing, an equation consisting of matrix and vectors for inverse problem is obtained. We use singular value decomposition method with damped factor to solve the equation. A synthetic data set which consists of several in-line profiles is used to test the methods. The results show that the methods are robust. Compared with the two dimensional method, the 3-D inversion method can give the right position of interfaces and the velocity structure when the crustal model is complicated.
基金State Natural Science Foundation of China (49734150).
文摘3-D velocity structure of P wave in the upper mantle beneath southwestern China and its adjacent areas (10°N [similar to] 36°N, 70°E [similar to] 110°E) down to the depth of 400 km has been studied by using 80 974 P-wave first arrival times recorded at 165 stations from 7 053 events both within the studying areas, selected from the ISC bulletin and the Bulletin of China and NEIC fundamental seismic network. With a resolution of grid spacing of 2°×2°, the velocity heterogeneity on the horizontal profile is obvious though it attenuates with the depth increasing. On the vertical profiles of velocity along the latitude of 16°N and 24°N, the collision and extrusion of India plate to Eurasia plate is displayed, and a remarkable velocity difference between India plate and Eurasia plate is shown. In the vertical profile along the longitude of 90°E, the subducting of India plate northward beneath Eurasia plate (Tibet plateau) is also obvious. On the horizontal profile at the depth of 90 km, a slow velocity stripe from Myitkyina, Myanmar to Donghai, Vietnam seems to be related to Honghe fault belt. An illustration method of describing the resolution more directly and exactly has been proposed and utilized in this paper.
基金supported by the National Science Fund for Distinguished Young Scholars(Grant No.40701071)
文摘A simplified numerical model of groundwater and solute transport is developed. At large scale area there exists a big spatial scale difference between horizontal and vertical length scales. In the resultant model, the seepage region is particularly divided into several virtual layers along the z direction and vertical 1-D columns covering x-y 2-D area according to stratum properties. The numerical algorithm is replacing the full 3-D water and mass balance analysis as the 2-D Galerkin finite element method in x- and y-directions and 1-D finite differential approach in the z direction. The reasonable method of giving minimum thickness is successfully used to handle transient change of water table, drying cells and problem of rewetting. The solution of the simplified model is preconditioned conjugate gradient and ORTHOMIN method. The validity of the developed 3-D groundwater model is tested with the typical pumping and backwater scenarios. Results of water balance of the computed example reveal the model computation reliability. Based on a representative 3-D pollution case, the solute transport module is tested against computing results using the MT3DMS. The capability and high efficiency to predict non-stationary situations of free groundwater surface and solute plume in regional scale problem is quantitatively investigated. It is shown that the proposed model is computationally effective.
基金the "95" Key Project (Grant No. KZ2952-51-410) from the Chinese Academy of Sciences and a project from the National Natural Science Foundation of China.
文摘By processing the CSND Rayleigh wave data with the matched filter FTAN technique, Rayleigh wave dispersion for southeast China is obtained. The 4°×4°S wave dispersion of the pure path is calculated using random inversion scheme, and 3-D S wave velocity structure is set up. Incorporating the above-mentioned results with wide angle seismic sounding data, we studied structure framework and the extending of faults in this area, which demonstrates that the depth of Moho in South China varies from 30 to 40 km, shallower from west to east. The depth of Moho varies from 25 to 28 km for the offshore. The depth of the asthenosphere in upper mantle varies from 60 to 100 km. The depth difference of layers at the two sides of Tanlu fault is more than 10 km at the south part of the Yangtze River, and the fault extends downward more than 170 km. The fault exceeds the main land at Hainan Island and slips into the southern China Sea. Both Tanlu fault and the huge bend of gravity gradient anomaly are influenced by
基金Project supported by the National Natural Science Foundation of China.
文摘On the basis of synthetic studies of geology, geophysics and geochemistry, the present Qinling Orogenic Belt can be described as a 3-D "flyover-type" geometric model with rheological layering structures. Furthermore, the tectonic dynamic analyses have been done based on the structural geometry and kinematic features. Thus it can be concluded that its present structure has resulted from an adjustment of deep-seated mantle dynamics and lithosphere coupling relationship since the Mesozoic-Cenozoic time.
