This paper calculates the parameters of image position and orientation,proposes a mathematical model and adopts a new method with three steps of transformations based on parallel ray projection.Every step of the model...This paper calculates the parameters of image position and orientation,proposes a mathematical model and adopts a new method with three steps of transformations based on parallel ray projection.Every step of the model is strict,and the map function of each transformation is the first order polynomials and other simple function.The final calculation of the parameters is for the linear equations with good status.As a result,the problem of the relativity of image parameter calculation is solved completely.Some experiments are carried out.展开更多
Full waveform inversion size of full waveform inversion will and the limitation of full waveform is mainly used to obtain high resolution velocity models of subsurface. The lead to a gigantic computation cost. Under t...Full waveform inversion size of full waveform inversion will and the limitation of full waveform is mainly used to obtain high resolution velocity models of subsurface. The lead to a gigantic computation cost. Under the available computer resource inversion, the authors propose L-BFGS algorithm as the optimization method to solve this problem. In order to demonstrate the flexibility of the method, three different numerical experi- ments have been done to analyze the properties of full waveform inversion based on L-BFGS.展开更多
Full waveform inversion( FWI) is a high resolution inversion method,which can reveal detailed information of the structure and lithology under complex geological background. It is limited by many kinds of noises when ...Full waveform inversion( FWI) is a high resolution inversion method,which can reveal detailed information of the structure and lithology under complex geological background. It is limited by many kinds of noises when the method applied to the real seismic data. Based on Huber function criterion,the objective function combinates the anti-noise of L1 norm and the stability of L2 norm in theory,the authors derive the gradient formula of the Huber function by using L-BFGS algorithm for FWI. The new method is proved by synthetic seismic data with the Gaussian noise and the impulse noise. Numerical test results show that L-BFGS algorithm is applied to the frequency domain FWI with the convergence speed and high calculation accuracy,and can effectively reduce computer memory usage; and the Huber function is more robust and stable than L2 norm even with the noises.展开更多
The new geoid model of Tanzania is a major breakthrough in the determination of vertical heights for the past 60 years. A new geoid model has been computed using precise gravity data supplemented with marine gravity d...The new geoid model of Tanzania is a major breakthrough in the determination of vertical heights for the past 60 years. A new geoid model has been computed using precise gravity data supplemented with marine gravity data from Gravity Recovery and Climate Experiment （GRACE） satellite and new high-resolution satellite gravity data from Gravity field and steady-state Ocean Circulation Explorer （GOCE）. A digital terrain model was also used based on 15＂ Shuttle and Rader Terrain Model （SRTM） data. The use of gravity data makes an orthometric height easily computed using Global Navigation Satellite System （GNSS）. The GNSS is being used in many African countries including Tanzania and soon it will replace conventional leveling technique and avoid frequent maintenance of benchmarks. This paper presents preliminary result of a new geoid model which shows that by using precious gravity data and Remove-Compute-Restore technique, an accuracy of 10 cm can be reached.展开更多
The decoding algorithms of two-dimensional Vernier anodes are deduced theoretically.The precision of decoding and uniqueness of encoding are proved.The influencing factors of detection sensitivity and spatial resoluti...The decoding algorithms of two-dimensional Vernier anodes are deduced theoretically.The precision of decoding and uniqueness of encoding are proved.The influencing factors of detection sensitivity and spatial resolution are discussed.The single photon imaging system is constructed,and the two-dimensional Vernier collector is fabricated.The image of the ultra-weak emission source is reconstructed.The spatial resolution of the system is about 100μm.展开更多
Recently, Sandia Laboratories developed a neutron scatter camera to detect special nuclear materials. This camera exhibits the following advantages: high efficiency, direction discrimination, neutron-gamma discriminat...Recently, Sandia Laboratories developed a neutron scatter camera to detect special nuclear materials. This camera exhibits the following advantages: high efficiency, direction discrimination, neutron-gamma discrimination ability, and wide field of view. However, using the direct projection method, the angular resolution of this camera is limited by uncertainties in the energies estimated from pulse height and time of flight measurements. In this study, we established an eight-element neutron scatter camera and conducted the experiment with a ^(252)Cf neutron source. The results show that it has an angular resolution better than 8°(1s) and a detection efficiency of approximately 2.6′10-4. Using maximum likelihood expectation maximization method, the image artifact was eliminated, and the angular resolution was improved. We proposed an average scattering angle method to estimate the scattering energy of neutrons and Compton gamma rays. As such, we can obtain a recognizable image and energy spectrum of the source with some degradation of energy and image resolutions. Finally, a newly measured light response function based on the MPD^(-4) device was used for image reconstruction. Although we did not obtain a better result than that of the standard light response function, we have observed the effects of light response function on image reconstruction.展开更多
In this work we compare the DNS results(Fabregat et al.2021,Fabregat et al.2021)for a mild cough already reported in the literarure with those obtained with a compressible URANS equations with a k-ϵturbulence model.In...In this work we compare the DNS results(Fabregat et al.2021,Fabregat et al.2021)for a mild cough already reported in the literarure with those obtained with a compressible URANS equations with a k-ϵturbulence model.In both cases,the dispersed phase has been modelled as spherical Lagrangian particles using the one-way coupling assumption.Overall,the URANS model is capable of reproducing the observed tendency of light particles under 64µm in diameter to rise due to the action of the drag exerted by the buoyant puff generated by the cough.Both DNS and URANS found that particles above 64µm will tend to describe parabolic trajectories under the action of gravitational forces.Grid independence analysis allows to qualify the impact of increasing mesh resolution on the particle cloud statistics as flow evolves.Results suggest that the k-ϵmodel overpredicts the horizontal displacement of the particles smaller than 64µm while the opposite occurs for the particles larger than 64µm.展开更多
文摘This paper calculates the parameters of image position and orientation,proposes a mathematical model and adopts a new method with three steps of transformations based on parallel ray projection.Every step of the model is strict,and the map function of each transformation is the first order polynomials and other simple function.The final calculation of the parameters is for the linear equations with good status.As a result,the problem of the relativity of image parameter calculation is solved completely.Some experiments are carried out.
文摘Full waveform inversion size of full waveform inversion will and the limitation of full waveform is mainly used to obtain high resolution velocity models of subsurface. The lead to a gigantic computation cost. Under the available computer resource inversion, the authors propose L-BFGS algorithm as the optimization method to solve this problem. In order to demonstrate the flexibility of the method, three different numerical experi- ments have been done to analyze the properties of full waveform inversion based on L-BFGS.
基金Supported by the National "863" Project(No.2014AA06A605)
文摘Full waveform inversion( FWI) is a high resolution inversion method,which can reveal detailed information of the structure and lithology under complex geological background. It is limited by many kinds of noises when the method applied to the real seismic data. Based on Huber function criterion,the objective function combinates the anti-noise of L1 norm and the stability of L2 norm in theory,the authors derive the gradient formula of the Huber function by using L-BFGS algorithm for FWI. The new method is proved by synthetic seismic data with the Gaussian noise and the impulse noise. Numerical test results show that L-BFGS algorithm is applied to the frequency domain FWI with the convergence speed and high calculation accuracy,and can effectively reduce computer memory usage; and the Huber function is more robust and stable than L2 norm even with the noises.
文摘The new geoid model of Tanzania is a major breakthrough in the determination of vertical heights for the past 60 years. A new geoid model has been computed using precise gravity data supplemented with marine gravity data from Gravity Recovery and Climate Experiment （GRACE） satellite and new high-resolution satellite gravity data from Gravity field and steady-state Ocean Circulation Explorer （GOCE）. A digital terrain model was also used based on 15＂ Shuttle and Rader Terrain Model （SRTM） data. The use of gravity data makes an orthometric height easily computed using Global Navigation Satellite System （GNSS）. The GNSS is being used in many African countries including Tanzania and soon it will replace conventional leveling technique and avoid frequent maintenance of benchmarks. This paper presents preliminary result of a new geoid model which shows that by using precious gravity data and Remove-Compute-Restore technique, an accuracy of 10 cm can be reached.
基金supported by the National Natural Science Foundation of China (Grant No. 10878005/A03)
文摘The decoding algorithms of two-dimensional Vernier anodes are deduced theoretically.The precision of decoding and uniqueness of encoding are proved.The influencing factors of detection sensitivity and spatial resolution are discussed.The single photon imaging system is constructed,and the two-dimensional Vernier collector is fabricated.The image of the ultra-weak emission source is reconstructed.The spatial resolution of the system is about 100μm.
基金supported by the National Natural Science Fundation of China(Grant Nos.1110510611375144&11275153)
文摘Recently, Sandia Laboratories developed a neutron scatter camera to detect special nuclear materials. This camera exhibits the following advantages: high efficiency, direction discrimination, neutron-gamma discrimination ability, and wide field of view. However, using the direct projection method, the angular resolution of this camera is limited by uncertainties in the energies estimated from pulse height and time of flight measurements. In this study, we established an eight-element neutron scatter camera and conducted the experiment with a ^(252)Cf neutron source. The results show that it has an angular resolution better than 8°(1s) and a detection efficiency of approximately 2.6′10-4. Using maximum likelihood expectation maximization method, the image artifact was eliminated, and the angular resolution was improved. We proposed an average scattering angle method to estimate the scattering energy of neutrons and Compton gamma rays. As such, we can obtain a recognizable image and energy spectrum of the source with some degradation of energy and image resolutions. Finally, a newly measured light response function based on the MPD^(-4) device was used for image reconstruction. Although we did not obtain a better result than that of the standard light response function, we have observed the effects of light response function on image reconstruction.
基金This work was supported by Spanish Ministerio de Ciencia,Innovación y Universidades(Grants Nos.RTI2018-100907-A-I00 and PID2020-113303GB-C21)the Generalitat de Catalunya(Grant No.2017-SGR-1234).
文摘In this work we compare the DNS results(Fabregat et al.2021,Fabregat et al.2021)for a mild cough already reported in the literarure with those obtained with a compressible URANS equations with a k-ϵturbulence model.In both cases,the dispersed phase has been modelled as spherical Lagrangian particles using the one-way coupling assumption.Overall,the URANS model is capable of reproducing the observed tendency of light particles under 64µm in diameter to rise due to the action of the drag exerted by the buoyant puff generated by the cough.Both DNS and URANS found that particles above 64µm will tend to describe parabolic trajectories under the action of gravitational forces.Grid independence analysis allows to qualify the impact of increasing mesh resolution on the particle cloud statistics as flow evolves.Results suggest that the k-ϵmodel overpredicts the horizontal displacement of the particles smaller than 64µm while the opposite occurs for the particles larger than 64µm.
