Mining operation, especially underground coal mining, always has the remarkable risks of ground control. Passive seismic velocity tomography based on simultaneous iterative reconstructive technique （SIRT） inversion ...Mining operation, especially underground coal mining, always has the remarkable risks of ground control. Passive seismic velocity tomography based on simultaneous iterative reconstructive technique （SIRT） inversion is used to deduce the stress redistribution around the longwall mining panel. The mining-induced microseismic events were recorded by mounting an array of receivers on the surface, above the active panel. After processing and filtering the seismic data, the three-dimensional tomography images of the p-wave velocity variations by SIRT passive seismic velocity tomography were provided. To display the velocity changes on coal seam level and subsequently to infer the stress redistribution, these three-dimensional tomograms into the coal seam level were sliced. In addition, the boundary element method （BEM） was used to simulate the stress redistribution. The results show that the inferred stresses from the passive seismic tomograms are conformed to numerical models and theoretical concept of the stress redistribution around the longwall panel. In velocity tomograms, the main zones of the stress redistribution arotmd the panel, including front and side abutment pressures, and gob stress are obvious and also the movement of stress zones along the face advancement is evident. Moreover, the effect of the advance rate of the face on the stress redistribution is demonstrated in tomography images. The research result proves that the SIRT passive seismic velocity tomography has an ultimate potential for monitoring the changes of stress redistribution around the longwall mining panel continuously and subsequently to improve safety of mining operations.展开更多
To better understand the mechanism of the Mw6.3 L＇Aquila （Central Italy） earthquake occurred in 2009, global positioning system （GPS） and interferometric synthetic aperture radar （InSAR） data were used to deriv...To better understand the mechanism of the Mw6.3 L＇Aquila （Central Italy） earthquake occurred in 2009, global positioning system （GPS） and interferometric synthetic aperture radar （InSAR） data were used to derive the coseismic slip distribution of the earthquake fault. Firstly, based on the homogeneous elastic half-space model, the fault geometric parameters were solved by the genetic algorithm. The best fitting model shows that the fault is a 13.7 km×14.1 km rectangular fault, in 139.3° strike direction and 50.2° southwest-dipping. Secondly, fixing the optimal fault geometric parameters, the fault plane was extended and discretized into 16× 16 patches, each with a size of 1 kmx 1 krn, and the non-uniform slip distribution of the fault was inverted by the steepest descent method with an appropriate smoothing ratio based on the layered crustal structure model. The preferred solution shows that the fault is mainly a normal fault with slight right-lateral strike slip, the maximum slip of 1.01 m is located in the depth of 8.28 km, the average rake is -100.9°, and the total geodetic moment is about 3.34× 1018 N.m （Mw 6.28）. The results are much closer than previous studies in comparison with the seismological estimation. These demonstrate that the coseismic fault slip distribution of the L＇Aauila earthauake inverted by the crustal model considering layered characters is reliable.展开更多
The classical elastic impedance （EI） inversion method, however, is based on the L2-norm misfit function and considerably sensitive to outliers, assuming the noise of the seismic data to be the Guassian-distribution....The classical elastic impedance （EI） inversion method, however, is based on the L2-norm misfit function and considerably sensitive to outliers, assuming the noise of the seismic data to be the Guassian-distribution. So we have developed a more robust elastic impedance inversion based on the Ll-norm misfit function, and the noise is assumed to be non-Gaussian. Meanwhile, some regularization methods including the sparse constraint regularization and elastic impedance point constraint regularization are incorporated to improve the ill-posed characteristics of the seismic inversion problem. Firstly, we create the Ll-norm misfit objective function of pre-stack inversion problem based on the Bayesian scheme within the sparse constraint regularization and elastic impedance point constraint regularization. And then, we obtain more robust elastic impedances of different angles which are less sensitive to outliers in seismic data by using the IRLS strategy. Finally, we extract the P-wave and S-wave velocity and density by using the more stable parameter extraction method. Tests on synthetic data show that the P-wave and S-wave velocity and density parameters are still estimated reasonable with moderate noise. A test on the real data set shows that compared to the results of the classical elastic impedance inversion method, the estimated results using the proposed method can get better lateral continuity and more distinct show of the gas, verifying the feasibility and stability of the method.展开更多
The aim of this work is to establish volcanic seismic reflection configuration models in the rift basins of Northeast China from a new perspective,the volcanostratigraphic structure.Accordingly,the volcanostratigraphi...The aim of this work is to establish volcanic seismic reflection configuration models in the rift basins of Northeast China from a new perspective,the volcanostratigraphic structure.Accordingly,the volcanostratigraphic structure of an outcrop near the Hailaier Rift Basin was analyzed to understand the characteristics and causal factors of physical boundaries.Further,3D seismic reflection data and analysis of deep boreholes in the Songliao Rift Basin were used to establish the relationship between volcanic seismic reflection configurations and volcanostratigraphic structures.These studies suggested that in volcanic successions,physical boundaries coincide with volcanic boundaries,and their distributions are controlled by the stacking patterns of volcanic units.Therefore,volcanic seismic reflection configurations can be interpreted in terms of the stacking patterns of volcanic units.These are also referred to as general bedding patterns in volcanostratigraphy.Furthermore,four typical seismic reflection configurations were identified,namely,the chaotic,the parallel continuous,the hummocky,the multi-mound superimposed and the composite.The corresponding interpretation models comprised single massive unit,vertical,intersectional,lateral multi-mound,and composite stacking patterns.The hummocky and composite reflection configurations with intersectional and composite stacking patterns are the most favorable for the exploration of volcanic reservoirs in rift basins.展开更多
Based on the scalar wave equation, making use of the ray approximation of the reflected seismic data (CMP or CSP gathers), the authors derive respectively the projection function of the primary waves and multiple wave...Based on the scalar wave equation, making use of the ray approximation of the reflected seismic data (CMP or CSP gathers), the authors derive respectively the projection function of the primary waves and multiple waves at the near offset (CMP or CSP gathers) in the parabolic Radon transform(PRT)domain. From the geometric point, the authors prove that the energy of the reflection still distributes along hyperbola which has higher curvature in the PRT domain and becomes some energy masses. So the primary waves and the multiple waves which interweave each other in ( x, t ) domain can be completely separated, which helps the multiple waves eliminated by filtering or muting. It is important for the analysis of velocity and the separator and elimination of multiple waves.展开更多
Replication is an approach often used to speed up the execution of queries submitted to a large dataset.A compile-time/run-time approach is presented for minimizing the response time of 2-dimensional range when a dist...Replication is an approach often used to speed up the execution of queries submitted to a large dataset.A compile-time/run-time approach is presented for minimizing the response time of 2-dimensional range when a distributed replica of a dataset exists.The aim is to partition the query payload(and its range) into subsets and distribute those to the replica nodes in a way that minimizes a client's response time.However,since query size and distribution characteristics of data(data dense/sparse regions) in varying ranges are not known a priori,performing efficient load balancing and parallel processing over the unpredictable workload is difficult.A technique based on the creation and manipulation of dynamic spatial indexes for query payload estimation in distributed queries was proposed.The effectiveness of this technique was demonstrated on queries for analysis of archived earthquake-generated seismic data records.展开更多
文摘Mining operation, especially underground coal mining, always has the remarkable risks of ground control. Passive seismic velocity tomography based on simultaneous iterative reconstructive technique （SIRT） inversion is used to deduce the stress redistribution around the longwall mining panel. The mining-induced microseismic events were recorded by mounting an array of receivers on the surface, above the active panel. After processing and filtering the seismic data, the three-dimensional tomography images of the p-wave velocity variations by SIRT passive seismic velocity tomography were provided. To display the velocity changes on coal seam level and subsequently to infer the stress redistribution, these three-dimensional tomograms into the coal seam level were sliced. In addition, the boundary element method （BEM） was used to simulate the stress redistribution. The results show that the inferred stresses from the passive seismic tomograms are conformed to numerical models and theoretical concept of the stress redistribution around the longwall panel. In velocity tomograms, the main zones of the stress redistribution arotmd the panel, including front and side abutment pressures, and gob stress are obvious and also the movement of stress zones along the face advancement is evident. Moreover, the effect of the advance rate of the face on the stress redistribution is demonstrated in tomography images. The research result proves that the SIRT passive seismic velocity tomography has an ultimate potential for monitoring the changes of stress redistribution around the longwall mining panel continuously and subsequently to improve safety of mining operations.
基金Projects(40974006,40774003) supported by the National Natural Science Foundation of ChinaProject(NCET-08-0570) supported by the Program for New Century Excellent Talents in Chinese Universities+2 种基金Projects(2011JQ001,2009QZZD004) supported by the Fundamental Research Funds for the Central Universities in ChinaProjects(09K005,09K006) supported by the Key Laboratory for Precise Engineering Surveying & Hazard Monitoring of Hunan Province,ChinaProject(1343-74334000023) supported by the Graduate DegreeThesis Innovation Foundation of Central South University,China
文摘To better understand the mechanism of the Mw6.3 L＇Aquila （Central Italy） earthquake occurred in 2009, global positioning system （GPS） and interferometric synthetic aperture radar （InSAR） data were used to derive the coseismic slip distribution of the earthquake fault. Firstly, based on the homogeneous elastic half-space model, the fault geometric parameters were solved by the genetic algorithm. The best fitting model shows that the fault is a 13.7 km×14.1 km rectangular fault, in 139.3° strike direction and 50.2° southwest-dipping. Secondly, fixing the optimal fault geometric parameters, the fault plane was extended and discretized into 16× 16 patches, each with a size of 1 kmx 1 krn, and the non-uniform slip distribution of the fault was inverted by the steepest descent method with an appropriate smoothing ratio based on the layered crustal structure model. The preferred solution shows that the fault is mainly a normal fault with slight right-lateral strike slip, the maximum slip of 1.01 m is located in the depth of 8.28 km, the average rake is -100.9°, and the total geodetic moment is about 3.34× 1018 N.m （Mw 6.28）. The results are much closer than previous studies in comparison with the seismological estimation. These demonstrate that the coseismic fault slip distribution of the L＇Aauila earthauake inverted by the crustal model considering layered characters is reliable.
基金Projects(U1562215,41674130,41404088)supported by the National Natural Science Foundation of ChinaProjects(2013CB228604,2014CB239201)supported by the National Basic Research Program of China+1 种基金Projects(2016ZX05027004-001,2016ZX05002006-009)supported by the National Oil and Gas Major Projects of ChinaProject(15CX08002A)supported by the Fundamental Research Funds for the Central Universities,China
文摘The classical elastic impedance （EI） inversion method, however, is based on the L2-norm misfit function and considerably sensitive to outliers, assuming the noise of the seismic data to be the Guassian-distribution. So we have developed a more robust elastic impedance inversion based on the Ll-norm misfit function, and the noise is assumed to be non-Gaussian. Meanwhile, some regularization methods including the sparse constraint regularization and elastic impedance point constraint regularization are incorporated to improve the ill-posed characteristics of the seismic inversion problem. Firstly, we create the Ll-norm misfit objective function of pre-stack inversion problem based on the Bayesian scheme within the sparse constraint regularization and elastic impedance point constraint regularization. And then, we obtain more robust elastic impedances of different angles which are less sensitive to outliers in seismic data by using the IRLS strategy. Finally, we extract the P-wave and S-wave velocity and density by using the more stable parameter extraction method. Tests on synthetic data show that the P-wave and S-wave velocity and density parameters are still estimated reasonable with moderate noise. A test on the real data set shows that compared to the results of the classical elastic impedance inversion method, the estimated results using the proposed method can get better lateral continuity and more distinct show of the gas, verifying the feasibility and stability of the method.
基金Projects(41472304,41430322) supported by the National Natural Science Foundation of ChinaProject(2012CB822002) supported by National Major State Basic Research Program of China
文摘The aim of this work is to establish volcanic seismic reflection configuration models in the rift basins of Northeast China from a new perspective,the volcanostratigraphic structure.Accordingly,the volcanostratigraphic structure of an outcrop near the Hailaier Rift Basin was analyzed to understand the characteristics and causal factors of physical boundaries.Further,3D seismic reflection data and analysis of deep boreholes in the Songliao Rift Basin were used to establish the relationship between volcanic seismic reflection configurations and volcanostratigraphic structures.These studies suggested that in volcanic successions,physical boundaries coincide with volcanic boundaries,and their distributions are controlled by the stacking patterns of volcanic units.Therefore,volcanic seismic reflection configurations can be interpreted in terms of the stacking patterns of volcanic units.These are also referred to as general bedding patterns in volcanostratigraphy.Furthermore,four typical seismic reflection configurations were identified,namely,the chaotic,the parallel continuous,the hummocky,the multi-mound superimposed and the composite.The corresponding interpretation models comprised single massive unit,vertical,intersectional,lateral multi-mound,and composite stacking patterns.The hummocky and composite reflection configurations with intersectional and composite stacking patterns are the most favorable for the exploration of volcanic reservoirs in rift basins.
文摘Based on the scalar wave equation, making use of the ray approximation of the reflected seismic data (CMP or CSP gathers), the authors derive respectively the projection function of the primary waves and multiple waves at the near offset (CMP or CSP gathers) in the parabolic Radon transform(PRT)domain. From the geometric point, the authors prove that the energy of the reflection still distributes along hyperbola which has higher curvature in the PRT domain and becomes some energy masses. So the primary waves and the multiple waves which interweave each other in ( x, t ) domain can be completely separated, which helps the multiple waves eliminated by filtering or muting. It is important for the analysis of velocity and the separator and elimination of multiple waves.
文摘Replication is an approach often used to speed up the execution of queries submitted to a large dataset.A compile-time/run-time approach is presented for minimizing the response time of 2-dimensional range when a distributed replica of a dataset exists.The aim is to partition the query payload(and its range) into subsets and distribute those to the replica nodes in a way that minimizes a client's response time.However,since query size and distribution characteristics of data(data dense/sparse regions) in varying ranges are not known a priori,performing efficient load balancing and parallel processing over the unpredictable workload is difficult.A technique based on the creation and manipulation of dynamic spatial indexes for query payload estimation in distributed queries was proposed.The effectiveness of this technique was demonstrated on queries for analysis of archived earthquake-generated seismic data records.
