3D traveltime calculation is widely used in seismic exploration technologies such as seismic migration and tomography. The fast marching method (FMM) is useful for calculating 3D traveltime and has proven to be effi...3D traveltime calculation is widely used in seismic exploration technologies such as seismic migration and tomography. The fast marching method (FMM) is useful for calculating 3D traveltime and has proven to be efficient and stable. However, it has low calculation accuracy near the source, which thus gives it low overall accuracy. This paper proposes a joint traveltime calculation method to solve this problem. The method firstly employs the wavefront construction method (WFC), which has a higher calculation accuracy than FMM in calculating traveltime in the small area near the source, and secondly adopts FMM to calculate traveltime for the remaining grid nodes. Due to the increase in calculation precision of grid nodes near the source, this new algorithm is shown to have good calculation precision while maintaining the high calculation efficiency of FMM, which is employed in most of the computational area. Results are verified using various numerical models.展开更多
An initial alignment technique for the strapdown inertial navigation system (SINS) of vehicles in the moving state is researched. By selecting an odometer as the system’s external sensor, the mathematical model for t...An initial alignment technique for the strapdown inertial navigation system (SINS) of vehicles in the moving state is researched. By selecting an odometer as the system’s external sensor, the mathematical model for the alignment in the moving state is established and the observability of the system is analyzed. The results show that the SINS can successfully achieve the precision alignment in 10 min when the vehicle is moving toward the prearranged place after its staying for several seconds to perform the coarse alignment. The precision of alignment can also be improved in the moving state compared with that in the static state.展开更多
The conventional acoustic logging interpretation method, which is based on vertical wells that penetrate isotropic formations, is not suitable for horizontal and deviated wells penetrating anisotropic formations. This...The conventional acoustic logging interpretation method, which is based on vertical wells that penetrate isotropic formations, is not suitable for horizontal and deviated wells penetrating anisotropic formations. This unsuitability is because during horizontal and deviated well drilling, cuttings will splash on the well wall or fall into the borehole bottom and form a thin bed of cuttings. In addition, the high velocity layers at different depths and intrinsic anisotropy may affect acoustic logging measurements. In this study, we examine how these factors affect the acoustic wave slowness measured in horizontal and deviated wells that are surrounded by an anisotropic medium using numerical simulation. We use the staggered-grid finite difference method in time domain (FDTD) combined with hybrid-PML. First, we acquire the acoustic slowness using a simulated array logging system, and then, we analyze how various factors affect acoustic slowness measurements and the differences between the effects of these factors. The factors considered are high-velocity layers, thin beds of cuttings, dipping angle, formation thickness, and anisotropy. The simulation results show that these factors affect acoustic wave slowness measurements differently. We observe that when the wavelength is much smaller than the distance between the borehole wall and high velocity layer, the true slowness of the formation could be acquired. When the wavelengths are of the same order (i.e., in the near-field scenarios), the geometrical acoustics theory is no longer applicable. Furthermore, when a thin bed of cuttings exists at the bottom of the borehole, Fermat's principle is still applicable, and true slowness can be acquired. In anisotropic formations, the measured slowness changes with increments in the dipping angle. Finally, for a measurement system with specific spacing, the slowness of a thin target layer can be acquired when the distance covered by the logging tool is sufficiently long. Based on systematical simulations with different dipping angles and anisotropy in homogenous TI media, slowness estimation charts are established to quantitatively determine the slowness at any dipping angle and for any value of the anisotropic ratio. Synthetic examples with different acoustic logging tools and different elastic parameters demonstrate that the acoustic slowness estimation method can be conveniently applied to horizontal and deviated wells in TI formations with high accuracy.展开更多
The basic principles of sea surface temperature (SST) remote sensing using infrared and microwave radiometers are introduced, and the differences between two sensors for retrieving sea surface temperature are invest...The basic principles of sea surface temperature (SST) remote sensing using infrared and microwave radiometers are introduced, and the differences between two sensors for retrieving sea surface temperature are investigated. The ground resolution, atmospheric effect, sea surface wind, skin depth and so on have important influence on precision of sea surface temperature retrieved by two sensors. The better understanding of the advantage and disadvantage of sea surface temperature detected by infrared and microwave radiometers would help us to imply SST remote sensing data more effectively and correctly.展开更多
A method for designing an X-ray flatness filter for medical electron linac is developed. It is used in the optimization process in the electron beam radiation system. Monte Carlo simulation method is used and two exam...A method for designing an X-ray flatness filter for medical electron linac is developed. It is used in the optimization process in the electron beam radiation system. Monte Carlo simulation method is used and two examples of real radiation system optimization processes for China-made medical electron linac are provided: 15 MV X- ray system of BJ-20 linac, and 12 MeV electron system of BJ-14. Results are verified by using the traditional method.展开更多
The slip-sweep technique is one of the high-efficiency, high-fidelity, and environmental vibroseis seismic prospecting techniques which consists of a vibrator group sweeping without waiting for the previous group's s...The slip-sweep technique is one of the high-efficiency, high-fidelity, and environmental vibroseis seismic prospecting techniques which consists of a vibrator group sweeping without waiting for the previous group's sweep to terminate. The cycle time can be reduced drastically and hence the production efficiency can be increased significantly but harmonic distortion of one sweep will leak into the record of the other sweep. In this paper, we propose an anti-correlation method for removing harmonic distortion in vibroseis data. This method is based on decomposition of the ground force signal into fundamental and harmonic components. Then the corresponding anti-correlation operator can be computed to estimate the energy of each harmonic after correlating the vibroseis data with the corresponding harmonic component. Finally, the vibroseis harmonic noise to be removed can be obtained by subtracting the extracted harmonic noise from the traces of the previous group's sweep. The advantage of the proposed method is that it can process both uncorrelated and correlated vibroseis seismic data. Moreover, the algorithm is simple, stable, and computationally fast. Especially, the significant contribution of this method is a considerable reduction in the harmonic without any alteration of the desired signals. The method was tested on both synthetic and field data sets to validate the good harmonic noise suppression results.展开更多
This paper deals with the research of the GPS/INS integrated navigation system applying Extended Kalman Filter, which involves integrated principles, scheme and technology of combining with real INS and GPS receiver d...This paper deals with the research of the GPS/INS integrated navigation system applying Extended Kalman Filter, which involves integrated principles, scheme and technology of combining with real INS and GPS receiver data. Emphases are placed on the modeling of system errors and implementation of the integrated system. Both loose and tightly coupled GPS/INS integrated in schemes are analyzed. On the basis of our experience accumulated in the research of GPS/INS for many years, the GPS/INS integrated navigation developing system is developed. It can be put into efficient and economic use in the study and design of integrated navigation system. It plays an important role in the aeronautical and astronautical fields in China. This system is not only a computer aided design software but also a semi physical simulation system by obtaining real INS and GPS receiver data. So the key software unit of the developing system could be conveniently transferred into practical engineering software in actual hardware integrated system. The application of this system shows that the design ideas and integrated scheme of this development system are successful, and can achieve good navigation result.展开更多
A fast U D factorization based Kalman filter for the 21 state integrated global positioning system and inertial navigation system (GPS/INS) is developed from the point of engineering implementation. The conventio...A fast U D factorization based Kalman filter for the 21 state integrated global positioning system and inertial navigation system (GPS/INS) is developed from the point of engineering implementation. The conventional Kalman filter is widely used for integration of GPS/INS, however, due to the model and numerical computation errors, the Kalman filter may diverge in engineering implementation. In order to solve this problem, an extended Kalman filter based on the U D factorization is proposed. Moreover, the high order integrated system suffers from the problem of long computation time, leading to difficulties in real time applications. An algorithmic approach is developed to improve the computational speed. A typical aircraft trajectory is simulated to compare the improvement in the computational speed and the navigation accuracy using the conventional Kalman filter and the fast Kalman filter based on the U D factorization. The results indicate that the methods proposed in this paper are very effective in overcoming these problems for the high dynamic integrated GPS/INS system.展开更多
A two-step method is proposed for detection and identification of invisible impact damage in composite structure under temperature changes using Lamb waves.First,a statistical outlier analysis is employed to distingui...A two-step method is proposed for detection and identification of invisible impact damage in composite structure under temperature changes using Lamb waves.First,a statistical outlier analysis is employed to distinguish whether the changes of Lamb wave signals are induced by damage within a monitoring area or are only affected by temperature changes.Damage indices are defined after the Lamb wave signals are processed by Fourier transform,and a Monte Carlo procedure is used to obtain the damage threshold value for the damage indices at the undamaged state.If the damage indices in the operation state exceed the threshold value,the presence of damage is determined.Then,a probabilistic damage imaging algorithm displaying probabilities of the presence of damage within the monitoring area is adopted to fuse information collected from multiple actuator-sensor paths to identify the location of damage.Damage indices under damaged state are used to generate the diagnostic image.Experimental study on a stiffened composite panel with random temperature changes is performed to demonstrate the effectiveness of the proposed method.展开更多
To simulate the actual flowfield at the exit of the supersonic/hypersonic inlet, a wind tunnel is designed to study the flow in the scramjet isolator under the asymmetric incoming flow. And compression fields in the i...To simulate the actual flowfield at the exit of the supersonic/hypersonic inlet, a wind tunnel is designed to study the flow in the scramjet isolator under the asymmetric incoming flow. And compression fields in the isolator are investigated using wall static and pitot pressure measurements. Three incoming Mach numbers are considered as 1.5, 1.8 and 2. Results show that the increase of the asymmetry of the flow at the isolator entrance leads to the increase of the shock train length in the isolator for a given pressure ratio. Based on the analysis of the flow asymmetry effect at the isolator entrance on the shock train length, a modified correlation is proposed to calculate the length of the shock train. Predicted results of the proposed correlation are in good agreement with the experimental data.展开更多
The problem of joint eigenvalue estimation for the non-defective commuting set of matrices A is addressed. A procedure revealing the joint eigenstructure by simultaneous diagonalization of. A with simultaneous Schur d...The problem of joint eigenvalue estimation for the non-defective commuting set of matrices A is addressed. A procedure revealing the joint eigenstructure by simultaneous diagonalization of. A with simultaneous Schur decomposition (SSD) and balance procedure alternately is proposed for performance considerations and also for overcoming the convergence difficulties of previous methods based only on simultaneous Schur form and unitary transformations, it is shown that the SSD procedure can be well incorporated with the balancing algorithm in a pingpong manner, i. e., each optimizes a cost function and at the same time serves as an acceleration procedure for the other. Under mild assumptions, the convergence of the two cost functions alternately optimized, i. e., the norm of A and the norm of the left-lower part of A is proved. Numerical experiments are conducted in a multi-dimensional harmonic retrieval application and suggest that the presented method converges considerably faster than the methods based on only unitary transformation for matrices which are not near to normality.展开更多
Ocean waves are the core environmental elements affecting the movements and structure design of ships. Statistical analysis of wave parameters is the basis for the establishment of long-term ship environmental adaptab...Ocean waves are the core environmental elements affecting the movements and structure design of ships. Statistical analysis of wave parameters is the basis for the establishment of long-term ship environmental adaptability prediction model. The observations from coastal stations, buoys, altimeters and volunteer ships that cover from 1993 to 2011 were interpolated into miller Ion-lat grids by using bilinear method and the analytical fields of ocean waves were given. By using optimal interpolation, the analysis wave fields were assimilated into the WAVEWATCH III (WW3) simulation results. From the assimilated results, the wave rose statistics, the wave height of muitiyear return period and the extreme 2-D wave spectrum are related to the ship seakeeping were calculated. Finally, the wave statistics in China offshore were analyzed in detail.展开更多
A speed sensorless vector control system of induction motor with estimated rotor speed and rotor flux using a new reduced order extended Kalman filter is proposed. With this method, two rotor flux components are sele...A speed sensorless vector control system of induction motor with estimated rotor speed and rotor flux using a new reduced order extended Kalman filter is proposed. With this method, two rotor flux components are selected as the state variables, and the rotor speed as an estimated parameter is regarded as an augmented state variable. The algorithm with reduced order decreases the computational complexity and makes the proposed estimator feasible to be implemented in real time. The simulation results show high accuracy of the estimation algorithm and good performance of speed control, and verify the usefulness of the proposed algorithm.展开更多
In order to measure the parameters of flight rocket by using radar,rocket impact point was estimated accurately for rocket trajectory correction.The Kalman filter with adaptive filter gain matrix was adopted.According...In order to measure the parameters of flight rocket by using radar,rocket impact point was estimated accurately for rocket trajectory correction.The Kalman filter with adaptive filter gain matrix was adopted.According to the particle trajectory model,the adaptive Kalman filter trajectory model was constructed for removing and filtering the outliers of the parameters during a section of flight detected by three-dimensional data radar and the rocket impact point was extrapolated.The results of numerical simulation show that the outliers and noise in trajectory measurement signal can be removed effectively by using the adaptive Kalman filter and the filter variance can converge in a short period of time.Based on the relation of filtering time and impact point estimation error,choosing the filtering time of 8-10 scan get the minimum estimation error of impact point.展开更多
基金supported by NSFC(Nos.41274120,41404085,and 41504084)
文摘3D traveltime calculation is widely used in seismic exploration technologies such as seismic migration and tomography. The fast marching method (FMM) is useful for calculating 3D traveltime and has proven to be efficient and stable. However, it has low calculation accuracy near the source, which thus gives it low overall accuracy. This paper proposes a joint traveltime calculation method to solve this problem. The method firstly employs the wavefront construction method (WFC), which has a higher calculation accuracy than FMM in calculating traveltime in the small area near the source, and secondly adopts FMM to calculate traveltime for the remaining grid nodes. Due to the increase in calculation precision of grid nodes near the source, this new algorithm is shown to have good calculation precision while maintaining the high calculation efficiency of FMM, which is employed in most of the computational area. Results are verified using various numerical models.
文摘An initial alignment technique for the strapdown inertial navigation system (SINS) of vehicles in the moving state is researched. By selecting an odometer as the system’s external sensor, the mathematical model for the alignment in the moving state is established and the observability of the system is analyzed. The results show that the SINS can successfully achieve the precision alignment in 10 min when the vehicle is moving toward the prearranged place after its staying for several seconds to perform the coarse alignment. The precision of alignment can also be improved in the moving state compared with that in the static state.
基金supported by National Natural Science Foundation of China(No.41204094)Science Foundation of China University of Petroleum,Beijing(No.2462015YQ0506)
文摘The conventional acoustic logging interpretation method, which is based on vertical wells that penetrate isotropic formations, is not suitable for horizontal and deviated wells penetrating anisotropic formations. This unsuitability is because during horizontal and deviated well drilling, cuttings will splash on the well wall or fall into the borehole bottom and form a thin bed of cuttings. In addition, the high velocity layers at different depths and intrinsic anisotropy may affect acoustic logging measurements. In this study, we examine how these factors affect the acoustic wave slowness measured in horizontal and deviated wells that are surrounded by an anisotropic medium using numerical simulation. We use the staggered-grid finite difference method in time domain (FDTD) combined with hybrid-PML. First, we acquire the acoustic slowness using a simulated array logging system, and then, we analyze how various factors affect acoustic slowness measurements and the differences between the effects of these factors. The factors considered are high-velocity layers, thin beds of cuttings, dipping angle, formation thickness, and anisotropy. The simulation results show that these factors affect acoustic wave slowness measurements differently. We observe that when the wavelength is much smaller than the distance between the borehole wall and high velocity layer, the true slowness of the formation could be acquired. When the wavelengths are of the same order (i.e., in the near-field scenarios), the geometrical acoustics theory is no longer applicable. Furthermore, when a thin bed of cuttings exists at the bottom of the borehole, Fermat's principle is still applicable, and true slowness can be acquired. In anisotropic formations, the measured slowness changes with increments in the dipping angle. Finally, for a measurement system with specific spacing, the slowness of a thin target layer can be acquired when the distance covered by the logging tool is sufficiently long. Based on systematical simulations with different dipping angles and anisotropy in homogenous TI media, slowness estimation charts are established to quantitatively determine the slowness at any dipping angle and for any value of the anisotropic ratio. Synthetic examples with different acoustic logging tools and different elastic parameters demonstrate that the acoustic slowness estimation method can be conveniently applied to horizontal and deviated wells in TI formations with high accuracy.
文摘The basic principles of sea surface temperature (SST) remote sensing using infrared and microwave radiometers are introduced, and the differences between two sensors for retrieving sea surface temperature are investigated. The ground resolution, atmospheric effect, sea surface wind, skin depth and so on have important influence on precision of sea surface temperature retrieved by two sensors. The better understanding of the advantage and disadvantage of sea surface temperature detected by infrared and microwave radiometers would help us to imply SST remote sensing data more effectively and correctly.
基金Supported by the National Natural Science Foundation of China (60672104,10675013)the Na-tional Basic Research Program of China ("973"Program)(2006CB705705)+1 种基金the 10th Five-Year Plan of the Ministry of Science and Technology of China(2001BA706B-05)the Joint Research Foundation of Beijing Municipal Commissionof Education~~
文摘A method for designing an X-ray flatness filter for medical electron linac is developed. It is used in the optimization process in the electron beam radiation system. Monte Carlo simulation method is used and two examples of real radiation system optimization processes for China-made medical electron linac are provided: 15 MV X- ray system of BJ-20 linac, and 12 MeV electron system of BJ-14. Results are verified by using the traditional method.
基金supported by the Sinopec Service Company and China National Petroleum Corporation
文摘The slip-sweep technique is one of the high-efficiency, high-fidelity, and environmental vibroseis seismic prospecting techniques which consists of a vibrator group sweeping without waiting for the previous group's sweep to terminate. The cycle time can be reduced drastically and hence the production efficiency can be increased significantly but harmonic distortion of one sweep will leak into the record of the other sweep. In this paper, we propose an anti-correlation method for removing harmonic distortion in vibroseis data. This method is based on decomposition of the ground force signal into fundamental and harmonic components. Then the corresponding anti-correlation operator can be computed to estimate the energy of each harmonic after correlating the vibroseis data with the corresponding harmonic component. Finally, the vibroseis harmonic noise to be removed can be obtained by subtracting the extracted harmonic noise from the traces of the previous group's sweep. The advantage of the proposed method is that it can process both uncorrelated and correlated vibroseis seismic data. Moreover, the algorithm is simple, stable, and computationally fast. Especially, the significant contribution of this method is a considerable reduction in the harmonic without any alteration of the desired signals. The method was tested on both synthetic and field data sets to validate the good harmonic noise suppression results.
文摘This paper deals with the research of the GPS/INS integrated navigation system applying Extended Kalman Filter, which involves integrated principles, scheme and technology of combining with real INS and GPS receiver data. Emphases are placed on the modeling of system errors and implementation of the integrated system. Both loose and tightly coupled GPS/INS integrated in schemes are analyzed. On the basis of our experience accumulated in the research of GPS/INS for many years, the GPS/INS integrated navigation developing system is developed. It can be put into efficient and economic use in the study and design of integrated navigation system. It plays an important role in the aeronautical and astronautical fields in China. This system is not only a computer aided design software but also a semi physical simulation system by obtaining real INS and GPS receiver data. So the key software unit of the developing system could be conveniently transferred into practical engineering software in actual hardware integrated system. The application of this system shows that the design ideas and integrated scheme of this development system are successful, and can achieve good navigation result.
文摘A fast U D factorization based Kalman filter for the 21 state integrated global positioning system and inertial navigation system (GPS/INS) is developed from the point of engineering implementation. The conventional Kalman filter is widely used for integration of GPS/INS, however, due to the model and numerical computation errors, the Kalman filter may diverge in engineering implementation. In order to solve this problem, an extended Kalman filter based on the U D factorization is proposed. Moreover, the high order integrated system suffers from the problem of long computation time, leading to difficulties in real time applications. An algorithmic approach is developed to improve the computational speed. A typical aircraft trajectory is simulated to compare the improvement in the computational speed and the navigation accuracy using the conventional Kalman filter and the fast Kalman filter based on the U D factorization. The results indicate that the methods proposed in this paper are very effective in overcoming these problems for the high dynamic integrated GPS/INS system.
基金Supported by the Aeronautical Science Foundation of China(2008ZA52012)the Six Kinds of Excellent Talent Project in Jiangsu Province of China(2010JZ004)the Research Foundation of Nanjing University of Aeronautics and Astronautics(NS2010027)~~
文摘A two-step method is proposed for detection and identification of invisible impact damage in composite structure under temperature changes using Lamb waves.First,a statistical outlier analysis is employed to distinguish whether the changes of Lamb wave signals are induced by damage within a monitoring area or are only affected by temperature changes.Damage indices are defined after the Lamb wave signals are processed by Fourier transform,and a Monte Carlo procedure is used to obtain the damage threshold value for the damage indices at the undamaged state.If the damage indices in the operation state exceed the threshold value,the presence of damage is determined.Then,a probabilistic damage imaging algorithm displaying probabilities of the presence of damage within the monitoring area is adopted to fuse information collected from multiple actuator-sensor paths to identify the location of damage.Damage indices under damaged state are used to generate the diagnostic image.Experimental study on a stiffened composite panel with random temperature changes is performed to demonstrate the effectiveness of the proposed method.
文摘To simulate the actual flowfield at the exit of the supersonic/hypersonic inlet, a wind tunnel is designed to study the flow in the scramjet isolator under the asymmetric incoming flow. And compression fields in the isolator are investigated using wall static and pitot pressure measurements. Three incoming Mach numbers are considered as 1.5, 1.8 and 2. Results show that the increase of the asymmetry of the flow at the isolator entrance leads to the increase of the shock train length in the isolator for a given pressure ratio. Based on the analysis of the flow asymmetry effect at the isolator entrance on the shock train length, a modified correlation is proposed to calculate the length of the shock train. Predicted results of the proposed correlation are in good agreement with the experimental data.
基金The National Natural Science Foundation of China(No.60572072,60496311),the National High Technology Researchand Development Program of China (863Program ) ( No.2003AA123310),the International Cooperation Project on Beyond 3G Mobile of China (No.2005DFA10360).
文摘The problem of joint eigenvalue estimation for the non-defective commuting set of matrices A is addressed. A procedure revealing the joint eigenstructure by simultaneous diagonalization of. A with simultaneous Schur decomposition (SSD) and balance procedure alternately is proposed for performance considerations and also for overcoming the convergence difficulties of previous methods based only on simultaneous Schur form and unitary transformations, it is shown that the SSD procedure can be well incorporated with the balancing algorithm in a pingpong manner, i. e., each optimizes a cost function and at the same time serves as an acceleration procedure for the other. Under mild assumptions, the convergence of the two cost functions alternately optimized, i. e., the norm of A and the norm of the left-lower part of A is proved. Numerical experiments are conducted in a multi-dimensional harmonic retrieval application and suggest that the presented method converges considerably faster than the methods based on only unitary transformation for matrices which are not near to normality.
基金supports from National Natural Science Foundation of China (No. 41406032 and No. 41376014)Open Fund of State Key Laboratory of Satellite Ocean Environment Dynamics (No. SOED1305)
文摘Ocean waves are the core environmental elements affecting the movements and structure design of ships. Statistical analysis of wave parameters is the basis for the establishment of long-term ship environmental adaptability prediction model. The observations from coastal stations, buoys, altimeters and volunteer ships that cover from 1993 to 2011 were interpolated into miller Ion-lat grids by using bilinear method and the analytical fields of ocean waves were given. By using optimal interpolation, the analysis wave fields were assimilated into the WAVEWATCH III (WW3) simulation results. From the assimilated results, the wave rose statistics, the wave height of muitiyear return period and the extreme 2-D wave spectrum are related to the ship seakeeping were calculated. Finally, the wave statistics in China offshore were analyzed in detail.
文摘A speed sensorless vector control system of induction motor with estimated rotor speed and rotor flux using a new reduced order extended Kalman filter is proposed. With this method, two rotor flux components are selected as the state variables, and the rotor speed as an estimated parameter is regarded as an augmented state variable. The algorithm with reduced order decreases the computational complexity and makes the proposed estimator feasible to be implemented in real time. The simulation results show high accuracy of the estimation algorithm and good performance of speed control, and verify the usefulness of the proposed algorithm.
文摘In order to measure the parameters of flight rocket by using radar,rocket impact point was estimated accurately for rocket trajectory correction.The Kalman filter with adaptive filter gain matrix was adopted.According to the particle trajectory model,the adaptive Kalman filter trajectory model was constructed for removing and filtering the outliers of the parameters during a section of flight detected by three-dimensional data radar and the rocket impact point was extrapolated.The results of numerical simulation show that the outliers and noise in trajectory measurement signal can be removed effectively by using the adaptive Kalman filter and the filter variance can converge in a short period of time.Based on the relation of filtering time and impact point estimation error,choosing the filtering time of 8-10 scan get the minimum estimation error of impact point.