A new method to accelerate the convergent rate of the space-alternatinggeneralized expectation-maximization (SAGE) algorithm is proposed. The new rescaled block-iterativeSAGE (RBI-SAGE) algorithm combines the RBI algo...A new method to accelerate the convergent rate of the space-alternatinggeneralized expectation-maximization (SAGE) algorithm is proposed. The new rescaled block-iterativeSAGE (RBI-SAGE) algorithm combines the RBI algorithm with the SAGE algorithm for PET imagereconstruction. In the new approach, the projection data is partitioned into disjoint blocks; eachiteration step involves only one of these blocks. SAGE updates the parameters sequentially in eachblock. In experiments, the RBI-SAGE algorithm and classical SAGE algorithm are compared in theapplication on positron emission tomography (PET) image reconstruction. Simulation results show thatRBI-SAGE has better performance than SAGE in both convergence and image quality.展开更多
The one-block version of ordered subsets (OS) techniques was used to accelerate the convergent rate of the space-alternating generalized expectation-maximization (SAGE) algorithm. The new row-action SAGE (RA-SAGE) alg...The one-block version of ordered subsets (OS) techniques was used to accelerate the convergent rate of the space-alternating generalized expectation-maximization (SAGE) algorithm. The new row-action SAGE (RA-SAGE) algorithm processed projections in sequentially orthogonal order which reduced the dependency among the projections and speeds up the convergences. Additionally, the over-relaxation parameter in the direction defined by the RA-SAGE algorithm was also applied to obtain fast convergence to a globally maximum likelihood (ML) solution. In experiments, the RA-SAGE algorithm and the classical SAGE algorithm were compared in the application to positron emission tomography (PET) image reconstruction. Simulation results showed that RA-SAGE had better performance than SAGE in both convergence and image quality.展开更多
A new method that uses a modified ordered subsets （MOS） algorithm to improve the convergence rate of space-alternating generalized expectation-maximization （SAGE） algorithm for positron emission tomography （PET）...A new method that uses a modified ordered subsets （MOS） algorithm to improve the convergence rate of space-alternating generalized expectation-maximization （SAGE） algorithm for positron emission tomography （PET） image reconstruction is proposed.In the MOS-SAGE algorithm,the number of projections and the access order of the subsets are modified in order to improve the quality of the reconstructed images and accelerate the convergence speed.The number of projections in a subset increases as follows：2,4,8,16,32 and 64.This sequence means that the high frequency component is recovered first and the low frequency component is recovered in the succeeding iteration steps.In addition,the neighboring subsets are separated as much as possible so that the correlation of projections can be decreased and the convergences can be speeded up.The application of the proposed method to simulated and real images shows that the MOS-SAGE algorithm has better performance than the SAGE algorithm and the OSEM algorithm in convergence and image quality.展开更多
The 5.12 Wenchuan Earthquake in 2008 induced hundreds of large-scale landslides. This paper systematically analyzes 112 large-scale landslides (surface area > 50000 m2), which were identified by interpretation of r...The 5.12 Wenchuan Earthquake in 2008 induced hundreds of large-scale landslides. This paper systematically analyzes 112 large-scale landslides (surface area > 50000 m2), which were identified by interpretation of remote sensing imagery and field investigations. The analysis suggests that the distribution of large-scale landslides is affected by the following four factors: (a) distance effect: 80% of studied large-scale landslides are located within a distance of 5 km from the seismic faults. The farther the distance to the faults, the lower the number of large-scale landslides; (b) locked segment effect: the large-scale landslides are mainly located in five concentration zones closely related with the crossing, staggering and transfer sections between one seismic fault section and the next one, as well as the end of the NE fault section. The zone with the highest concentration was the Hongbai-Chaping segment, where a great number of large-scale landslides including the two largest landslides were located. The second highest concentration of large-scale landslides was observed in the Nanba-Donghekou segment at the end of NE fault, where the Donghekou landslide and the Woqian landslide occurred; (c) Hanging wall effect: about 70% of the large-scale landslides occurred on the hanging wall of the seismic faults; and (d) direction effect: in valleys perpendicular to the seismic faults, the density of large-scale landslides on the slopes facing the seismic wave is obviously higher than that on the slopes dipping in the same direction as the direction of propagation of the seismic wave. Meanwhile, it is found that the sliding and moving directions of large-scale landslides are related to the staggering direction of the faults in each section. In Qingchuan County where the main fault activity was horizontal twisting and staggering, a considerable number of landslides showed the feature of sliding and moving in NE direction which coincides with the staggering direction of the seismic faults.展开更多
Reconstructed synsedimentary paleogeomorpholgy is an effective method in predicting the distribution of sandbodies and can provide valid evidence in the search for reservoirs. Based on a synsedimentary paleogeomorphol...Reconstructed synsedimentary paleogeomorpholgy is an effective method in predicting the distribution of sandbodies and can provide valid evidence in the search for reservoirs. Based on a synsedimentary paleogeomorpholgy reconstruction of the third member of the Dongying formation (Ed3) in the Qikou sag,the basic paleogeomorphic characteristics of Ed3 are described and the spatial distributions of denudation and subsidence areas are discussed. Key boundary faults controlling the deposition of important sediment bodies are proposed and the interrelations among faulted slope-break belts or flexure belts,intra-depression uplifts and sags,sediment sources and sediments input points have been investigated.展开更多
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.展开更多
In order to explore the influence of coal mining disturbance on the rockburst occurring in fault zone, this research constructed a mechanical model for the evolution of fault stress, and analyzed the influence of the ...In order to explore the influence of coal mining disturbance on the rockburst occurring in fault zone, this research constructed a mechanical model for the evolution of fault stress, and analyzed the influence of the ratio of horizontal stress to vertical stress on the stability of fault, and the spatial distribution of the stress in fault zone as well as its evolution rule. Furthermore, the rockburst danger at different spatial areas of fault zone was predicated. Results show that: when both sides of the working face are mined out, the fault zone in the working face presents greater horizontal and vertical stresses at its boundaries but exhibits smaller stresses in its middle section; however, the ratio of horizontal stress to vertical stress is found to be greater at middle section and smaller at boundaries. As the working face advances towards the fault, the horizontal and vertical stresses of the fault firstly increases and then decreases; conversely,the ratio of horizontal stress to vertical stress keeps decreasing all the time. Therefore, if the fault zones are closer to the goaf and the coal wall, the stress ratio will be smaller, and the fault slip will be more prone to occur, therefore rockburst danger will be greater. This research results provide guidance for the rockburst prevention and hazard control of the coal mining in fault zone.展开更多
Simulation and optimization were applied to a capacitive sensor system based on electrical tomography technology. Sensors, consisting of Morgantown Energy Technology Center (METC) axial synchro driving guard electrode...Simulation and optimization were applied to a capacitive sensor system based on electrical tomography technology. Sensors, consisting of Morgantown Energy Technology Center (METC) axial synchro driving guard electrodes and two sets of detecting electrodes, make it possible to obtain simultaneously two groups of signals of the void fraction in oil-gas two-phase flow. The computational and experimental results showed that available sensors, charactered by high resolution and fast real-time response can be used for real-time liquid-gas two-phase flow pattern determination.展开更多
For the diagnosis of glaucoma,optical coherence tomography(OCT)is a noninvasive imaging technique for the assessment of retinal layers.To accurately segment intraretinal layers in an optic nerve head(ONH)region,we pro...For the diagnosis of glaucoma,optical coherence tomography(OCT)is a noninvasive imaging technique for the assessment of retinal layers.To accurately segment intraretinal layers in an optic nerve head(ONH)region,we proposed an automatic method for the segmentation of three intraretinal layers in eye OCT scans centered on ONH.The internal limiting membrane,inner segment and outer segment,Bruch’s membrane surfaces under vascular shadows,and interaction of multiple high-reflectivity regions in the OCT image can be accurately segmented through this method.Then,we constructed a novel spatial-gradient continuity constraint,termed spatial-gradient continuity constraint,for the correction of discontinuity between adjacent image segmentation results.In our experiment,we randomly selected 20 B-scans,each annotated three retinal layers by experts.Signed distance errors of?0.80μm obtained through this method are lower than those obtained through the state-of-art method(?1.43μm).Meanwhile,the segmentation results can be used as bases for the diagnosis of glaucoma.展开更多
In the proximity of an active fault, spatial variation of peak ground motion is significantly affected by the faulting mechanism. It has been observed that near fault ground motions consists of different characteristi...In the proximity of an active fault, spatial variation of peak ground motion is significantly affected by the faulting mechanism. It has been observed that near fault ground motions consists of different characteristics compared to the far fault ground motions. Near fault records, in the distance range of less than 100 m from the faults are not available except for few cases. Therefore numerical simulation of ground motions for such near-fault situations is necessary. In addition to the understanding of the phenomenon of near fault ground motion there is a need to enhance our understanding of the possible potential hazard that can be caused due to the future rupture activity by understanding the phenomenon of surface faulting. In this paper we propose numerical simulation based on discrete modeling to investigate the fault rupture propagation and its effect on the surface peak ground acceleration. In the present two dimensional study rupture propagation due to bedrock motion has been observed for different shear wave velocity. A model of size 1000× 150 m is selected for this purpose. It has been observed that as the stiffness of the media is decreasing, the affected surface is decreasing and also width of the shear crack zone is decreasing. Secondly, we attempted to study the ground motion on the surface due to the bedrock motion in presence of boulders in the soil media by keeping the boulder at different positions. We find that there is an increase in the shear zone as well as the PGA on the surface when the boulder is present on the foot wall and in the vicinity of the rupture zone. Finally, we performed the analysis using layered media and studied the affect of crack propagation and also the variation of peak accelerations. Findings from the study can be utilized to assess the damage potential of the near fault areas.展开更多
The rupture dimensions of earthquake faults are important parameters for characterizing earthquake ruptures and ground motions. Two key parameters to be determined are the rupture depth and dip angle of earthquake fau...The rupture dimensions of earthquake faults are important parameters for characterizing earthquake ruptures and ground motions. Two key parameters to be determined are the rupture depth and dip angle of earthquake faults. Dislocation theory in an elastic half space indicates that if a seismic rupture directly runs up to the ground surface, there exist zero points of horizontal strain in the surface deformation, which correspond to the rupture depths, except for pure strike-slip faults. In this study, we use numerical simulations to investigate the possibility of inferring rupture depths from zero-strain points for cases of buried faults and heterogeneous media. The results show that the correspondence of zero-strain points to the rupture depths can be influenced by the heterogeneity of the underground media and the stress field. For buried faults, the correspondence relationship is approximately valid when the fault depth is <1 km. In addition, the range of earthquake fault dip angles can be estimated by horizontal displacements on the ground. We also study how to determine the rupture depths of faults from InSAR data after large earthquakes, and successfully apply the method to the 2008 Wenchuan earthquake. The method proposed here, which determines the parameters of fault geometry according to surface deformation, is simple and easy to perform. With independent of aftershocks, it can provide valuable constraints to kinematic inversions.展开更多
文摘A new method to accelerate the convergent rate of the space-alternatinggeneralized expectation-maximization (SAGE) algorithm is proposed. The new rescaled block-iterativeSAGE (RBI-SAGE) algorithm combines the RBI algorithm with the SAGE algorithm for PET imagereconstruction. In the new approach, the projection data is partitioned into disjoint blocks; eachiteration step involves only one of these blocks. SAGE updates the parameters sequentially in eachblock. In experiments, the RBI-SAGE algorithm and classical SAGE algorithm are compared in theapplication on positron emission tomography (PET) image reconstruction. Simulation results show thatRBI-SAGE has better performance than SAGE in both convergence and image quality.
文摘The one-block version of ordered subsets (OS) techniques was used to accelerate the convergent rate of the space-alternating generalized expectation-maximization (SAGE) algorithm. The new row-action SAGE (RA-SAGE) algorithm processed projections in sequentially orthogonal order which reduced the dependency among the projections and speeds up the convergences. Additionally, the over-relaxation parameter in the direction defined by the RA-SAGE algorithm was also applied to obtain fast convergence to a globally maximum likelihood (ML) solution. In experiments, the RA-SAGE algorithm and the classical SAGE algorithm were compared in the application to positron emission tomography (PET) image reconstruction. Simulation results showed that RA-SAGE had better performance than SAGE in both convergence and image quality.
基金The National Basic Research Program of China (973Program) (No.2003CB716102).
文摘A new method that uses a modified ordered subsets （MOS） algorithm to improve the convergence rate of space-alternating generalized expectation-maximization （SAGE） algorithm for positron emission tomography （PET） image reconstruction is proposed.In the MOS-SAGE algorithm,the number of projections and the access order of the subsets are modified in order to improve the quality of the reconstructed images and accelerate the convergence speed.The number of projections in a subset increases as follows：2,4,8,16,32 and 64.This sequence means that the high frequency component is recovered first and the low frequency component is recovered in the succeeding iteration steps.In addition,the neighboring subsets are separated as much as possible so that the correlation of projections can be decreased and the convergences can be speeded up.The application of the proposed method to simulated and real images shows that the MOS-SAGE algorithm has better performance than the SAGE algorithm and the OSEM algorithm in convergence and image quality.
基金sponsored by the project of the Chinese National Key Basic Research Program on "The failure mechanism and distribution rule of slopes under strong earthquakes" (Grant No. 2008CB425801)the Education Department Innovation Research Team Program (Grant No. IRT0812)
文摘The 5.12 Wenchuan Earthquake in 2008 induced hundreds of large-scale landslides. This paper systematically analyzes 112 large-scale landslides (surface area > 50000 m2), which were identified by interpretation of remote sensing imagery and field investigations. The analysis suggests that the distribution of large-scale landslides is affected by the following four factors: (a) distance effect: 80% of studied large-scale landslides are located within a distance of 5 km from the seismic faults. The farther the distance to the faults, the lower the number of large-scale landslides; (b) locked segment effect: the large-scale landslides are mainly located in five concentration zones closely related with the crossing, staggering and transfer sections between one seismic fault section and the next one, as well as the end of the NE fault section. The zone with the highest concentration was the Hongbai-Chaping segment, where a great number of large-scale landslides including the two largest landslides were located. The second highest concentration of large-scale landslides was observed in the Nanba-Donghekou segment at the end of NE fault, where the Donghekou landslide and the Woqian landslide occurred; (c) Hanging wall effect: about 70% of the large-scale landslides occurred on the hanging wall of the seismic faults; and (d) direction effect: in valleys perpendicular to the seismic faults, the density of large-scale landslides on the slopes facing the seismic wave is obviously higher than that on the slopes dipping in the same direction as the direction of propagation of the seismic wave. Meanwhile, it is found that the sliding and moving directions of large-scale landslides are related to the staggering direction of the faults in each section. In Qingchuan County where the main fault activity was horizontal twisting and staggering, a considerable number of landslides showed the feature of sliding and moving in NE direction which coincides with the staggering direction of the seismic faults.
基金Projects 40872077 supported by the National Natural Science Foundation of China2008CDA098 by the Key Natural Science Foundation of Hubei Province
文摘Reconstructed synsedimentary paleogeomorpholgy is an effective method in predicting the distribution of sandbodies and can provide valid evidence in the search for reservoirs. Based on a synsedimentary paleogeomorpholgy reconstruction of the third member of the Dongying formation (Ed3) in the Qikou sag,the basic paleogeomorphic characteristics of Ed3 are described and the spatial distributions of denudation and subsidence areas are discussed. Key boundary faults controlling the deposition of important sediment bodies are proposed and the interrelations among faulted slope-break belts or flexure belts,intra-depression uplifts and sags,sediment sources and sediments input points have been investigated.
基金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.
基金Financial support for this work,provided by the National Natural Science Foundation of China(Nos.51104150,51204165,51404243)the Twelfth Five-Year National Science and Technology Support Program of China(No.2012BAK09B01)the Priority Academic Program Development of Jiangsu Higher Education Institutions of China,are gratefully acknowledged
文摘In order to explore the influence of coal mining disturbance on the rockburst occurring in fault zone, this research constructed a mechanical model for the evolution of fault stress, and analyzed the influence of the ratio of horizontal stress to vertical stress on the stability of fault, and the spatial distribution of the stress in fault zone as well as its evolution rule. Furthermore, the rockburst danger at different spatial areas of fault zone was predicated. Results show that: when both sides of the working face are mined out, the fault zone in the working face presents greater horizontal and vertical stresses at its boundaries but exhibits smaller stresses in its middle section; however, the ratio of horizontal stress to vertical stress is found to be greater at middle section and smaller at boundaries. As the working face advances towards the fault, the horizontal and vertical stresses of the fault firstly increases and then decreases; conversely,the ratio of horizontal stress to vertical stress keeps decreasing all the time. Therefore, if the fault zones are closer to the goaf and the coal wall, the stress ratio will be smaller, and the fault slip will be more prone to occur, therefore rockburst danger will be greater. This research results provide guidance for the rockburst prevention and hazard control of the coal mining in fault zone.
基金Project (No. 2002AA616050) supported by the Hi-Tech Researchand Development Program (863) of China
文摘Simulation and optimization were applied to a capacitive sensor system based on electrical tomography technology. Sensors, consisting of Morgantown Energy Technology Center (METC) axial synchro driving guard electrodes and two sets of detecting electrodes, make it possible to obtain simultaneously two groups of signals of the void fraction in oil-gas two-phase flow. The computational and experimental results showed that available sensors, charactered by high resolution and fast real-time response can be used for real-time liquid-gas two-phase flow pattern determination.
基金Projects(61672542,61573380)supported by the National Natural Science Foundation of China
文摘For the diagnosis of glaucoma,optical coherence tomography(OCT)is a noninvasive imaging technique for the assessment of retinal layers.To accurately segment intraretinal layers in an optic nerve head(ONH)region,we proposed an automatic method for the segmentation of three intraretinal layers in eye OCT scans centered on ONH.The internal limiting membrane,inner segment and outer segment,Bruch’s membrane surfaces under vascular shadows,and interaction of multiple high-reflectivity regions in the OCT image can be accurately segmented through this method.Then,we constructed a novel spatial-gradient continuity constraint,termed spatial-gradient continuity constraint,for the correction of discontinuity between adjacent image segmentation results.In our experiment,we randomly selected 20 B-scans,each annotated three retinal layers by experts.Signed distance errors of?0.80μm obtained through this method are lower than those obtained through the state-of-art method(?1.43μm).Meanwhile,the segmentation results can be used as bases for the diagnosis of glaucoma.
文摘In the proximity of an active fault, spatial variation of peak ground motion is significantly affected by the faulting mechanism. It has been observed that near fault ground motions consists of different characteristics compared to the far fault ground motions. Near fault records, in the distance range of less than 100 m from the faults are not available except for few cases. Therefore numerical simulation of ground motions for such near-fault situations is necessary. In addition to the understanding of the phenomenon of near fault ground motion there is a need to enhance our understanding of the possible potential hazard that can be caused due to the future rupture activity by understanding the phenomenon of surface faulting. In this paper we propose numerical simulation based on discrete modeling to investigate the fault rupture propagation and its effect on the surface peak ground acceleration. In the present two dimensional study rupture propagation due to bedrock motion has been observed for different shear wave velocity. A model of size 1000× 150 m is selected for this purpose. It has been observed that as the stiffness of the media is decreasing, the affected surface is decreasing and also width of the shear crack zone is decreasing. Secondly, we attempted to study the ground motion on the surface due to the bedrock motion in presence of boulders in the soil media by keeping the boulder at different positions. We find that there is an increase in the shear zone as well as the PGA on the surface when the boulder is present on the foot wall and in the vicinity of the rupture zone. Finally, we performed the analysis using layered media and studied the affect of crack propagation and also the variation of peak accelerations. Findings from the study can be utilized to assess the damage potential of the near fault areas.
基金supported by the National Natural Science Foundation of China (Grant Nos. 41074070, 41174035)the SinoProbe Program (Grant No. SinoProbe-08-01)
文摘The rupture dimensions of earthquake faults are important parameters for characterizing earthquake ruptures and ground motions. Two key parameters to be determined are the rupture depth and dip angle of earthquake faults. Dislocation theory in an elastic half space indicates that if a seismic rupture directly runs up to the ground surface, there exist zero points of horizontal strain in the surface deformation, which correspond to the rupture depths, except for pure strike-slip faults. In this study, we use numerical simulations to investigate the possibility of inferring rupture depths from zero-strain points for cases of buried faults and heterogeneous media. The results show that the correspondence of zero-strain points to the rupture depths can be influenced by the heterogeneity of the underground media and the stress field. For buried faults, the correspondence relationship is approximately valid when the fault depth is <1 km. In addition, the range of earthquake fault dip angles can be estimated by horizontal displacements on the ground. We also study how to determine the rupture depths of faults from InSAR data after large earthquakes, and successfully apply the method to the 2008 Wenchuan earthquake. The method proposed here, which determines the parameters of fault geometry according to surface deformation, is simple and easy to perform. With independent of aftershocks, it can provide valuable constraints to kinematic inversions.
