We develop a new approach to estimating bottom parameters based on the Bayesian theory in deep ocean. The solution in a Bayesian inversion is characterized by its posterior probability density (PPD), which combines ...We develop a new approach to estimating bottom parameters based on the Bayesian theory in deep ocean. The solution in a Bayesian inversion is characterized by its posterior probability density (PPD), which combines prior information about the model with information from an observed data set. Bottom parameters are sensitive to the transmission loss (TL) data in shadow zones of deep ocean. In this study, TLs of different frequencies from the South China Sea in the summer of 2014 are used as the observed data sets. The interpretation of the multidimensional PPD requires the calculation of its moments, such as the mean, covariance, and marginal distributions, which provide parameter estimates and uncertainties. Considering that the sensitivities of shallow- zone TLs vary for different frequencies of the bottom parameters in the deep ocean, this research obtains bottom parameters at varying frequencies. Then, the inversion results are compared with the sampling data and the correlations between bottom parameters are determined. Furthermore, we show the inversion results for multi- frequency combined inversion. The inversion results are verified by the experimental TLs and the numerical results, which are calculated using the inverted bottom parameters for different source depths and receiver depths at the corresponding frequency.展开更多
The forward design of trajectory planning strategies requires preset trajectory optimization functions,resulting in poor adaptability of the strategy and an inability to accurately generate obstacle avoidance trajecto...The forward design of trajectory planning strategies requires preset trajectory optimization functions,resulting in poor adaptability of the strategy and an inability to accurately generate obstacle avoidance trajectories that conform to real driver behavior habits.In addition,owing to the strong time-varying dynamic characteristics of obstacle avoidance scenarios,it is necessary to design numerous trajectory optimization functions and adjust the corresponding parameters.Therefore,an anthropomorphic obstacle-avoidance trajectory planning strategy for adaptive driving scenarios is proposed.First,numerous expert-demonstrated trajectories are extracted from the HighD natural driving dataset.Subsequently,a trajectory expectation feature-matching algorithm is proposed that uses maximum entropy inverse reinforcement learning theory to learn the extracted expert-demonstrated trajectories and achieve automatic acquisition of the optimization function of the expert-demonstrated trajectory.Furthermore,a mapping model is constructed by combining the key driving scenario information that affects vehicle obstacle avoidance with the weight of the optimization function,and an anthropomorphic obstacle avoidance trajectory planning strategy for adaptive driving scenarios is proposed.Finally,the proposed strategy is verified based on real driving scenarios.The results show that the strategy can adjust the weight distribution of the trajectory optimization function in real time according to the“emergency degree”of obstacle avoidance and the state of the vehicle.Moreover,this strategy can generate anthropomorphic trajectories that are similar to expert-demonstrated trajectories,effectively improving the adaptability and acceptability of trajectories in driving scenarios.展开更多
This paper first establishes the prior globe dynamical model by geophysics,which is a solid earth elastic deformation model.Then,the parameters of the globe dynamic model can be obtained by inverting the geodetic data...This paper first establishes the prior globe dynamical model by geophysics,which is a solid earth elastic deformation model.Then,the parameters of the globe dynamic model can be obtained by inverting the geodetic data.The inverse method can be used in seismology and geology,and to make earthquake prediction.展开更多
Analysis of a four-dimensional displacement vector on the fabric of space-time in the special or general case into two Four-dimensional vectors, according to specific conditions leads to the splitting of the total fab...Analysis of a four-dimensional displacement vector on the fabric of space-time in the special or general case into two Four-dimensional vectors, according to specific conditions leads to the splitting of the total fabric of space-time into a positive subspace-time that represents the space of causality and a negative subspace-time which represents a space without causality, thus, in the special case, we have new transformations for the coordinates of space and time modified from Lorentz transformations specific to each subspace, where the contraction of length disappears and the speed of light is no longer a universal constant. In the general case, we have new types of matric tensor, one for positive subspace-time and the other for negative subspace-time. We also find that the speed of the photon decreases in positive subspace-time until it reaches zero and increases in negative subspace-time until it reaches the speed of light when the photon reaches the Schwarzschild radius.展开更多
According to inverse heat transfer theory, the evolutions of synthetic surface heat transfer coefficient(SSHTC) of the quenching surface of 7B50 alloy during water-spray quenching were simulated by the Pro CAST soft...According to inverse heat transfer theory, the evolutions of synthetic surface heat transfer coefficient(SSHTC) of the quenching surface of 7B50 alloy during water-spray quenching were simulated by the Pro CAST software based on accurate cooling curves measured by the modified Jominy specimen and temperature-dependent thermo-physical properties of 7 B50 alloy calculated using the JMat Pro software. Results show that the average cooling rate at 6 mm from the quenching surface and 420-230 ℃(quench sensitive temperature range) is 45.78℃/s. The peak-value of the SSHTC is 69 kW/(m^2·K) obtained at spray quenching for 0.4 s and the corresponding temperature of the quenching surface is 160 ℃. In the initial stage of spray quenching, the phenomenon called "temperature plateau" appears on the cooling curve of the quenching surface. The temperature range of this plateau is 160-170℃ with the duration about 3 s. During the temperature plateau, heat transfer mechanism of the quenching surface transforms from nucleate boiling regime to single-phase convective regime.展开更多
We use magnetic field data observed by the Swarm mission from 2014 to 2020 to construct,for the first time,a two-dimensional(2 D)lithospheric magnetic anomaly model of Egypt and its surrounding area.Nighttime data dur...We use magnetic field data observed by the Swarm mission from 2014 to 2020 to construct,for the first time,a two-dimensional(2 D)lithospheric magnetic anomaly model of Egypt and its surrounding area.Nighttime data during quiet geomagnetic conditions has been expanded in terms of the Legendre polynomial in harmonic terms N=6-50.The damped least square method has been used to estimate the model coefficients based on the lithospheric magnetic data.Modeled data at two different altitudes(438-448 km and 503-511 km)were compared with the CHAOS model.Results exhibit that the 2 D model is superior to the CHAOS model in the capability of extracting more information about small-scale crustal anomaly field.At low altitudes(438-448 km),the strength of the anomaly field increases,but the noise of the external fields has greatly reduced at high altitudes(503-511 km).Besides,the magnetic anomaly field at low altitudes has illuminated short-scale anomalies that didn’t appear at high altitudes.Both the total and vertical magnetic anomaly vectors showed their ability to reveal tectonic structures compared with Moho depth map and the geological maps.展开更多
In this paper,a new earthquake location method based on the waveform inversion is proposed.As is known to all,the waveform misfit function under the L2 measure is suffering from the cycle skipping problem.This leads t...In this paper,a new earthquake location method based on the waveform inversion is proposed.As is known to all,the waveform misfit function under the L2 measure is suffering from the cycle skipping problem.This leads to a very small convergence domain of the conventional waveform based earthquake location methods.In present study,by introducing and solving two simple sub-optimization problems,we greatly expand the convergence domain of the waveform based earthquake location method.According to a large number of numerical experiments,the new method expands the range of convergence by several tens of times.This allows us to locate the earthquake accurately even from some relatively bad initial values.展开更多
A review of ten-year's practice in developing the improved simultaneous physical retrieval method(ISPRM)is given in the hope that some creative ideas can be drawn from it.The improvement upon the SPRM is associate...A review of ten-year's practice in developing the improved simultaneous physical retrieval method(ISPRM)is given in the hope that some creative ideas can be drawn from it.The improvement upon the SPRM is associated with the under-determinedness of this ill-posed inverse problem.In our experiment,the precondition is observed that prior information must be independent of the satellite measurements.The well-posed retrieval theory has told us that the forward process is fundamental for the retrieval,and it is the bridge between the input of satellite radiance and the output of retrievals.In order to obtain a better result from the forward process. the full advantage of every prior information available must be taken.It is necessary to turn the ill- posed inverse problem into the well-posed one.Then by using the Ridge regression or Bayes algorithm to find the optimal combination among the first guess,the theoretical analogue information and the satellite observations,the impact of the under-determinedness of this inverse problem on the numerical solution is minimized.展开更多
In the second paper on the inverse relativity model, we explained in the first paper [1] that analyzing the four-dimensional displacement vector on space-time according to a certain approach leads to the splitting of ...In the second paper on the inverse relativity model, we explained in the first paper [1] that analyzing the four-dimensional displacement vector on space-time according to a certain approach leads to the splitting of space-time into positive and negative subspace-time. Here, in the second paper, we continue to analyze each of the four-dimensional vectors of velocity, acceleration, momentum, and forces on the total space-time fabric. According to the approach followed in the first paper. As a result, in the special case, we obtain new transformations for each of the velocity, acceleration, momentum, energy, and forces specific to each subspace-time, which are subject to the positive and negative modified Lorentz transformations described in the first paper. According to these transformations, momentum remains a conserved quantity in the positive subspace and increases in the negative subspace, while the relativistic total energy decreases in the positive subspace and increases in the negative subspace. In the general case, we also have new types of energy-momentum tensor, one for positive subspace-time and the other for negative subspace-time, where the energy density decreases in positive subspace-time and increases in negative subspace-time, and we also obtain new gravitational field equations for each subspace-time.展开更多
In this paper,we apply the Wasserstein-Fisher-Rao(WFR)metric from the unbalanced optimal transport theory to the earthquake location problem.Compared with the quadratic Wasserstein(W2)metric from the classical optimal...In this paper,we apply the Wasserstein-Fisher-Rao(WFR)metric from the unbalanced optimal transport theory to the earthquake location problem.Compared with the quadratic Wasserstein(W2)metric from the classical optimal transport theory,the advantage of this method is that it retains the important amplitude information as a new constraint,which avoids the problem of the degeneration of the optimization objective function near the real earthquake hypocenter and origin time.As a result,the deviation of the global minimum of the optimization objective function based on the WFR metric from the true solution can be much smaller than the results based on the W2 metric when there exists strong data noise.Thus,we develop an accurate earthquake location method under strong data noise.Many numerical experiments verify our conclusions.展开更多
基金Supported by the National Natural Science Foundation of China under Grant No 11174235
文摘We develop a new approach to estimating bottom parameters based on the Bayesian theory in deep ocean. The solution in a Bayesian inversion is characterized by its posterior probability density (PPD), which combines prior information about the model with information from an observed data set. Bottom parameters are sensitive to the transmission loss (TL) data in shadow zones of deep ocean. In this study, TLs of different frequencies from the South China Sea in the summer of 2014 are used as the observed data sets. The interpretation of the multidimensional PPD requires the calculation of its moments, such as the mean, covariance, and marginal distributions, which provide parameter estimates and uncertainties. Considering that the sensitivities of shallow- zone TLs vary for different frequencies of the bottom parameters in the deep ocean, this research obtains bottom parameters at varying frequencies. Then, the inversion results are compared with the sampling data and the correlations between bottom parameters are determined. Furthermore, we show the inversion results for multi- frequency combined inversion. The inversion results are verified by the experimental TLs and the numerical results, which are calculated using the inverted bottom parameters for different source depths and receiver depths at the corresponding frequency.
基金supported by the National Natural Science Foundation of China(51875302)。
文摘The forward design of trajectory planning strategies requires preset trajectory optimization functions,resulting in poor adaptability of the strategy and an inability to accurately generate obstacle avoidance trajectories that conform to real driver behavior habits.In addition,owing to the strong time-varying dynamic characteristics of obstacle avoidance scenarios,it is necessary to design numerous trajectory optimization functions and adjust the corresponding parameters.Therefore,an anthropomorphic obstacle-avoidance trajectory planning strategy for adaptive driving scenarios is proposed.First,numerous expert-demonstrated trajectories are extracted from the HighD natural driving dataset.Subsequently,a trajectory expectation feature-matching algorithm is proposed that uses maximum entropy inverse reinforcement learning theory to learn the extracted expert-demonstrated trajectories and achieve automatic acquisition of the optimization function of the expert-demonstrated trajectory.Furthermore,a mapping model is constructed by combining the key driving scenario information that affects vehicle obstacle avoidance with the weight of the optimization function,and an anthropomorphic obstacle avoidance trajectory planning strategy for adaptive driving scenarios is proposed.Finally,the proposed strategy is verified based on real driving scenarios.The results show that the strategy can adjust the weight distribution of the trajectory optimization function in real time according to the“emergency degree”of obstacle avoidance and the state of the vehicle.Moreover,this strategy can generate anthropomorphic trajectories that are similar to expert-demonstrated trajectories,effectively improving the adaptability and acceptability of trajectories in driving scenarios.
文摘This paper first establishes the prior globe dynamical model by geophysics,which is a solid earth elastic deformation model.Then,the parameters of the globe dynamic model can be obtained by inverting the geodetic data.The inverse method can be used in seismology and geology,and to make earthquake prediction.
文摘Analysis of a four-dimensional displacement vector on the fabric of space-time in the special or general case into two Four-dimensional vectors, according to specific conditions leads to the splitting of the total fabric of space-time into a positive subspace-time that represents the space of causality and a negative subspace-time which represents a space without causality, thus, in the special case, we have new transformations for the coordinates of space and time modified from Lorentz transformations specific to each subspace, where the contraction of length disappears and the speed of light is no longer a universal constant. In the general case, we have new types of matric tensor, one for positive subspace-time and the other for negative subspace-time. We also find that the speed of the photon decreases in positive subspace-time until it reaches zero and increases in negative subspace-time until it reaches the speed of light when the photon reaches the Schwarzschild radius.
基金Project(2016YFB0300801)supported by the National Key Research and Development Program of ChinaProject(51371045)supported by the National Natural Science Foundation of China
文摘According to inverse heat transfer theory, the evolutions of synthetic surface heat transfer coefficient(SSHTC) of the quenching surface of 7B50 alloy during water-spray quenching were simulated by the Pro CAST software based on accurate cooling curves measured by the modified Jominy specimen and temperature-dependent thermo-physical properties of 7 B50 alloy calculated using the JMat Pro software. Results show that the average cooling rate at 6 mm from the quenching surface and 420-230 ℃(quench sensitive temperature range) is 45.78℃/s. The peak-value of the SSHTC is 69 kW/(m^2·K) obtained at spray quenching for 0.4 s and the corresponding temperature of the quenching surface is 160 ℃. In the initial stage of spray quenching, the phenomenon called "temperature plateau" appears on the cooling curve of the quenching surface. The temperature range of this plateau is 160-170℃ with the duration about 3 s. During the temperature plateau, heat transfer mechanism of the quenching surface transforms from nucleate boiling regime to single-phase convective regime.
文摘We use magnetic field data observed by the Swarm mission from 2014 to 2020 to construct,for the first time,a two-dimensional(2 D)lithospheric magnetic anomaly model of Egypt and its surrounding area.Nighttime data during quiet geomagnetic conditions has been expanded in terms of the Legendre polynomial in harmonic terms N=6-50.The damped least square method has been used to estimate the model coefficients based on the lithospheric magnetic data.Modeled data at two different altitudes(438-448 km and 503-511 km)were compared with the CHAOS model.Results exhibit that the 2 D model is superior to the CHAOS model in the capability of extracting more information about small-scale crustal anomaly field.At low altitudes(438-448 km),the strength of the anomaly field increases,but the noise of the external fields has greatly reduced at high altitudes(503-511 km).Besides,the magnetic anomaly field at low altitudes has illuminated short-scale anomalies that didn’t appear at high altitudes.Both the total and vertical magnetic anomaly vectors showed their ability to reveal tectonic structures compared with Moho depth map and the geological maps.
基金This work was supported by the National Nature Science Foundation of China(Grant Nos.41230210,41390452)Hao Wu was also partially supported by the National Nature Science Foundation of China(Grant Nos.11101236,91330203)and SRF for ROCS,SEM.The authors are grateful to Prof.Shi Jin for his helpful suggestions and discussions that greatly improve the presentation.Hao Wu would like to thank Prof.Ping Tong for his valuable comments.The authors would also like to thank the referees for their valuable suggestions which helped to improve the content and presentation of this paper.
文摘In this paper,a new earthquake location method based on the waveform inversion is proposed.As is known to all,the waveform misfit function under the L2 measure is suffering from the cycle skipping problem.This leads to a very small convergence domain of the conventional waveform based earthquake location methods.In present study,by introducing and solving two simple sub-optimization problems,we greatly expand the convergence domain of the waveform based earthquake location method.According to a large number of numerical experiments,the new method expands the range of convergence by several tens of times.This allows us to locate the earthquake accurately even from some relatively bad initial values.
基金Supported by NNSF of China under Grant(49794030#)National"973"Program No.4 (G1998040909#).
文摘A review of ten-year's practice in developing the improved simultaneous physical retrieval method(ISPRM)is given in the hope that some creative ideas can be drawn from it.The improvement upon the SPRM is associated with the under-determinedness of this ill-posed inverse problem.In our experiment,the precondition is observed that prior information must be independent of the satellite measurements.The well-posed retrieval theory has told us that the forward process is fundamental for the retrieval,and it is the bridge between the input of satellite radiance and the output of retrievals.In order to obtain a better result from the forward process. the full advantage of every prior information available must be taken.It is necessary to turn the ill- posed inverse problem into the well-posed one.Then by using the Ridge regression or Bayes algorithm to find the optimal combination among the first guess,the theoretical analogue information and the satellite observations,the impact of the under-determinedness of this inverse problem on the numerical solution is minimized.
文摘In the second paper on the inverse relativity model, we explained in the first paper [1] that analyzing the four-dimensional displacement vector on space-time according to a certain approach leads to the splitting of space-time into positive and negative subspace-time. Here, in the second paper, we continue to analyze each of the four-dimensional vectors of velocity, acceleration, momentum, and forces on the total space-time fabric. According to the approach followed in the first paper. As a result, in the special case, we obtain new transformations for each of the velocity, acceleration, momentum, energy, and forces specific to each subspace-time, which are subject to the positive and negative modified Lorentz transformations described in the first paper. According to these transformations, momentum remains a conserved quantity in the positive subspace and increases in the negative subspace, while the relativistic total energy decreases in the positive subspace and increases in the negative subspace. In the general case, we also have new types of energy-momentum tensor, one for positive subspace-time and the other for negative subspace-time, where the energy density decreases in positive subspace-time and increases in negative subspace-time, and we also obtain new gravitational field equations for each subspace-time.
基金National Natural Science Foundation of China(Grant Nos.11871297,11971258,U1839206)National Key Research and Development Program of China on Monitoring,Early Warning and Prevention of Major Natural Disaster(Grant No.2017YFC1500301)Tsinghua University Initiative Scientific Research Program.
文摘In this paper,we apply the Wasserstein-Fisher-Rao(WFR)metric from the unbalanced optimal transport theory to the earthquake location problem.Compared with the quadratic Wasserstein(W2)metric from the classical optimal transport theory,the advantage of this method is that it retains the important amplitude information as a new constraint,which avoids the problem of the degeneration of the optimization objective function near the real earthquake hypocenter and origin time.As a result,the deviation of the global minimum of the optimization objective function based on the WFR metric from the true solution can be much smaller than the results based on the W2 metric when there exists strong data noise.Thus,we develop an accurate earthquake location method under strong data noise.Many numerical experiments verify our conclusions.