In the adjoint-state method, the forward-propagated source wavefield and the backward-propagated receiver wavefield must be available simultaneously either for seismic imaging in migration or for gradient calculation ...In the adjoint-state method, the forward-propagated source wavefield and the backward-propagated receiver wavefield must be available simultaneously either for seismic imaging in migration or for gradient calculation in inversion. A feasible way to avoid the excessive storage demand is to reconstruct the source wavefield backward in time by storing the entire history of the wavefield in perfectly matched layers. In this paper, we make full use of the elementwise global property of the Laplace operator of the spectral element method (SEM) and propose an efficient source wavefield reconstruction method at the cost of storing the wavefield history only at single boundary layer nodes. Numerical experiments indicate that the accuracy of the proposed method is identical to that of the conventional method and is independent of the order of the Lagrange polynomials, the element type, and the temporal discretization method. In contrast, the memory-saving ratios of the conventional method versus our method is at least N when using either quadrilateral or hexahedron elements, respectively, where N is the order of the Lagrange polynomials used in the SEM. A higher memorysaving ratio is achieved with triangular elements versus quadrilaterals. The new method is applied to reverse time migration by considering the Marmousi model as a benchmark. Numerical results demonstrate that the method is able to provide the same result as the conventional method but with about 1/25 times lower storage demand. With the proposed wavefield reconstruction method, the storage demand is dramatically reduced;therefore, in-core memory storage is feasible even for large-scale three-dimensional adjoint inversion problems.展开更多
A two-stage source reconstruction algorithm for bioluminescence tomography (BLT) is developed using hybrid finite element method (FEM). The proposed algorithm takes full advantages of linear and quadratic FEMs, which ...A two-stage source reconstruction algorithm for bioluminescence tomography (BLT) is developed using hybrid finite element method (FEM). The proposed algorithm takes full advantages of linear and quadratic FEMs, which can be used to localize and quantify bioluminescent source accurately. In the first stage, a large permissible region is roughly determined and then iteratively evolved to reduce matrix dimension using efficient linear FEM. In the final stage, high-convergence quadratic FEM is applied to improve reconstruction result. Both numerical simulation and physical experiment are performed to evaluate the proposed algorithm. The relevant results demonstrate that quantitative reconstruction can be well achieved in terms of computation efficiency, source position, power density, and total power when compared with previous studies.展开更多
The accuracy of the background optical properties has a considerable effect on the quality of reconstructed images in near-infrared functional brain imaging based on continuous wave diffuse optical tomography(CW-DOT...The accuracy of the background optical properties has a considerable effect on the quality of reconstructed images in near-infrared functional brain imaging based on continuous wave diffuse optical tomography(CW-DOT). We propose a region stepwise reconstruction method in CW-DOT scheme for reconstructing the background absorption and reduced scattering coefficients of the two-layered slab sample with the known geometric information. According to the relation between the thickness of the top layer and source– detector separation, the conventional measurement data are divided into two groups and are employed to reconstruct the top and bottom background optical properties, respectively. The numerical simulation results demonstrate that the proposed method can reconstruct the background optical properties of two-layered slab sample effectively. The region-of-interest reconstruction results are better than those of the conventional simultaneous reconstruction method.展开更多
In this work,we present a high-order discontinuous Galerkin method for the shallow water equations incorporating horizontal temperature gradients(also known as the Ripa model),which exactly maintains the lake at rest ...In this work,we present a high-order discontinuous Galerkin method for the shallow water equations incorporating horizontal temperature gradients(also known as the Ripa model),which exactly maintains the lake at rest steady state.Herein,we propose original numerical fluxes defined on the basis of the hydrostatic reconstruction idea and a simple source term approximation.This novel approach allows us to achieve the well-balancing of the discontinuous Galerkin method without complication.Moreover,the proposed method retains genuinely high-order accuracy for smooth solutions and it shows good resolution for discontinuous solutions at the same time.Rigorous numerical analysis as well as extensive numerical results all verify the good performances of the proposed method.展开更多
基金financial support for this work contributed by the National Key Research and Development Program of China (grant numbers 2016YFC0600101 and 2016YFC 0600201)the National Natural Science Foundation of China (grant numbers 41874065, 41604076, 41674102, 41674095, 41522401, 41574082, and 41774097)
摘要In the adjoint-state method, the forward-propagated source wavefield and the backward-propagated receiver wavefield must be available simultaneously either for seismic imaging in migration or for gradient calculation in inversion. A feasible way to avoid the excessive storage demand is to reconstruct the source wavefield backward in time by storing the entire history of the wavefield in perfectly matched layers. In this paper, we make full use of the elementwise global property of the Laplace operator of the spectral element method (SEM) and propose an efficient source wavefield reconstruction method at the cost of storing the wavefield history only at single boundary layer nodes. Numerical experiments indicate that the accuracy of the proposed method is identical to that of the conventional method and is independent of the order of the Lagrange polynomials, the element type, and the temporal discretization method. In contrast, the memory-saving ratios of the conventional method versus our method is at least N when using either quadrilateral or hexahedron elements, respectively, where N is the order of the Lagrange polynomials used in the SEM. A higher memorysaving ratio is achieved with triangular elements versus quadrilaterals. The new method is applied to reverse time migration by considering the Marmousi model as a benchmark. Numerical results demonstrate that the method is able to provide the same result as the conventional method but with about 1/25 times lower storage demand. With the proposed wavefield reconstruction method, the storage demand is dramatically reduced;therefore, in-core memory storage is feasible even for large-scale three-dimensional adjoint inversion problems.
基金supported by National Basic Research Program of China (973 Program) (No. 2011CB707702)National Natural Science Foundation of China (Nos. 81090272, 81000632, and 30900334)+1 种基金the Shaanxi Provincial Natural Science Foundation (No. 2009JQ8018)the Fundamental Research Funds for the Central Universities
摘要A two-stage source reconstruction algorithm for bioluminescence tomography (BLT) is developed using hybrid finite element method (FEM). The proposed algorithm takes full advantages of linear and quadratic FEMs, which can be used to localize and quantify bioluminescent source accurately. In the first stage, a large permissible region is roughly determined and then iteratively evolved to reduce matrix dimension using efficient linear FEM. In the final stage, high-convergence quadratic FEM is applied to improve reconstruction result. Both numerical simulation and physical experiment are performed to evaluate the proposed algorithm. The relevant results demonstrate that quantitative reconstruction can be well achieved in terms of computation efficiency, source position, power density, and total power when compared with previous studies.
基金supported by the National Natural Science Foundation of China(Nos.81271618 and 81371602)the Tianjin Municipal Government of China(Nos.12JCQNJC09400 and 13JCZDJC28000)the Research Fund for the Doctoral Program of Higher Education of China(No.20120032110056)
摘要The accuracy of the background optical properties has a considerable effect on the quality of reconstructed images in near-infrared functional brain imaging based on continuous wave diffuse optical tomography(CW-DOT). We propose a region stepwise reconstruction method in CW-DOT scheme for reconstructing the background absorption and reduced scattering coefficients of the two-layered slab sample with the known geometric information. According to the relation between the thickness of the top layer and source– detector separation, the conventional measurement data are divided into two groups and are employed to reconstruct the top and bottom background optical properties, respectively. The numerical simulation results demonstrate that the proposed method can reconstruct the background optical properties of two-layered slab sample effectively. The region-of-interest reconstruction results are better than those of the conventional simultaneous reconstruction method.
基金support of the Natural Science Foundation of China through Grants No.11771228The author Qiang Niu is supported by the XJTLU research enhancement fund with No.REF-18-01-04 and the Key Programme Special Fund(KSF)in XJTLU with Nos.KSF-E-32,KSF-E-21 and KSF-P-02.
摘要In this work,we present a high-order discontinuous Galerkin method for the shallow water equations incorporating horizontal temperature gradients(also known as the Ripa model),which exactly maintains the lake at rest steady state.Herein,we propose original numerical fluxes defined on the basis of the hydrostatic reconstruction idea and a simple source term approximation.This novel approach allows us to achieve the well-balancing of the discontinuous Galerkin method without complication.Moreover,the proposed method retains genuinely high-order accuracy for smooth solutions and it shows good resolution for discontinuous solutions at the same time.Rigorous numerical analysis as well as extensive numerical results all verify the good performances of the proposed method.