Bézier曲线间最近距离的计算方法

针对Bzier曲线间最近距离计算问题,提出一种简捷、可靠的计算方法.该方法以Bernstein多项式算术运算为工具,建立Bzier曲线间最近距离的计算模型;然后充分利用Bzier曲面的凸包性质和de Casteljau分割算法进行求解.该方法几何意义明确,能有效地避免迭代初始值的选择和非线性方程组的求解,并可进一步推广应用于计算Bzier曲线/曲面间的最近距离.实验结果表明,该方法简捷、可靠且容易实现,与Newton-Raphson方法的融合可进一步提高该方法的运行速度.

基于广义正交多项式迭积微分算子的地震波场数值模拟研究

Seismic Numerical Modeling is one of bases of the Exploratory Seismology and Academic Seismology, also is a research field in great demand. Essence of seismic numerical modeling is to assume that structure and parameters of the underground media model are known, simulate the wave-field and calculate the numerical seismic record that should be observed. Seismic numerical modeling is not only a means to know the seismic wave-field in complex inhomogeneous media, but also a test to the application effect by all kinds of methods. There are many seismic numerical modeling methods, each method has its own merits and drawbacks. During the forward modeling, the computation precision and the efficiency are two pivotal questions to evaluate the validity and superiority of the method. The target of my dissertation is to find a new method to possibly improve the computation precision and efficiency, and apply the new forward method to modeling the wave-field in the complex inhomogeneous media. Convolutional Forsyte polynomial differentiator (CFPD) approach developed in this dissertation is robust and efficient, it shares some of the advantages of the high precision of generalized orthogonal polynomial and the high speed of the short operator finite-difference. By adjusting the operator length and optimizing the operator coefficient, the method can involve whole and local information of the wave-field. One of main tasks of the dissertation is to develop a creative, generalized and high precision method. The author introduce convolutional Forsyte polynomial differentiator to calculate the spatial derivative of seismic wave equation, and apply the time staggered grid finite-difference which can better meet the high precision of the convolutional differentiator to substitute the conventional finite-difference to calculate the time derivative of seismic wave equation, then creating a new forward method to modeling the wave-field in complex inhomogeneous media. Comparing with Fourier pseudo-spectral method, Chebyshev pseudo-spectral method, staggered- grid finite difference method and finite element method, convolutional Forsyte polynomial differentiator (CFPD) method has many advantages: 1. Comparing with Fourier pseudo-spectral method. Fourier pseudo-spectral method (FPS) is a local operator, its results have Gibbs effects when the media parameters change, then arose great errors. Therefore, Fourier pseudo-spectral method can not deal with special complex and random heterogeneous media. But convolutional Forsyte polynomial differentiator method can cover global and local information. So for complex inhomogeneous media, CFPD is more efficient. 2. Comparing with staggered-grid high-order finite-difference method, CFPD takes less dots than FD at single wave length, and the number does not increase with the widening of the studying area. 3. Comparing with Chebyshev pseudo-spectral method (CPS). The calculation region of Chebyshev pseudo-spectral method is fixed in , under the condition of unchangeable precision, the augmentation of calculation is unacceptable. Thus Chebyshev pseudo-spectral method is inapplicable to large area. CFPD method is more applicable to large area. 4. Comparing with finite element method (FE), CFPD can use lager grids. The other task of this dissertation is to study 2.5 dimension (2.5D) seismic wave-field. The author reviews the development and present situation of 2.5D problem, expatiates the essentiality of studying the 2.5D problem, apply CFPD method to simulate the seismic wave-field in 2.5D inhomogeneous media. The results indicate that 2.5D numerical modeling is efficient to simulate one of the sections of 3D media, 2.5D calculation is much less time-consuming than 3D calculation, and the wave dispersion of 2.5D modeling is obviously less than that of 3D modeling. Question on applying time staggered-grid convolutional differentiator based on CFPD to modeling 2.5D complex inhomogeneous media was not studied by any geophysicists before, it is a fire-new creation absolutely. The theory and practices prove that the new method can efficiently model the seismic wave-field in complex media. Proposing and developing this new method can provide more choices to study the seismic wave-field modeling, seismic wave migration, seismic inversion, and seismic wave imaging.

复杂条件地层滤波预测的理论与关键技术研究

The theory researches of prediction about stratigraphic filtering in complex condition are carried out, and three key techniques are put forward in this dissertation. Theoretical aspects: The prediction equations for both slant incidence in horizontally layered medium and that in laterally variant velocity medium are expressed appropriately. Solving the equations, the linear prediction operator of overlaid layers, then corresponding reflection/transmission operators, can be obtained. The properties of linear prediction operator are elucidated followed by putting forward the event model for generalized Goupillaud layers. Key technique 1: Spectral factorization is introduced to solve the prediction equations in complex condition and numerical results are illustrated. Key technique 2: So-called large-step wavefield extrapolation of one-way wave under laterally variant velocity circumstance is studied. Based on Lie algebraic integral and structure preserving algorithm, large-step wavefield depth extrapolation scheme is set forth. In this method, the complex phase of wavefield extrapolation operator’s symbol is expressed as a linear combination of wavenumbers with the coefficients of this linear combination in the form of the integral of interval velocity and its derivatives over depth. The exponential transform of the complex phase is implemented through phase shifting, BCH splitting and orthogonal polynomial expansion. The results of numerical test show that large-step scheme takes on a great number of advantages as low accumulating error, cheapness, well adaptability to laterally variant velocity, small dispersive, etc. Key technique 3: Utilizing large-step wavefield extrapolation scheme and based on the idea of local harmonic decomposition, the technique generating angle gathers for 2D case is generalized to 3D case so as to solve the problems generating and storing 3D prestack angle gathers. Shot domain parallel scheme is adopted by which main duty for servant-nodes is to compute trigonometric expansion coefficients, while that for host-node is to reclaim them with which object-oriented angle gathers yield. In theoretical research, many efforts have been made in probing into the traits of uncertainties within macro-dynamic procedures.

电离层行进式扰动的GPS台网监测研究

In this dissertation, we investigated two types of traveling ionospheric disturbances (TIDs)/gravity waves (GWs) triggered separately by auroral energy input during super geomagnetic storms and solar terminator (ST) under quiet geomagnetic conditions (kp<3+) using TEC measurements from the global network of GPS receivers. Research into the generation and propagation of TIDs/GWs during storms greatly enhance our understandings on the evolution processes of energy transportation from the high-latitude’s magnetosphere to the low-latitude ionosphere and the conjugated effect of TIDs propagation between the northern and southern hemispheres. Our results revealed that the conjugacy of propagation direction between the northern and southern hemispheres was subject to the influence of Coriolis force. We also figure out the evolution processes of ionospheric disturbances at the global scale. These are important topics that had not been well addressed previously. In addition, we also obtained thee wave structures of medium scale TIDs excited by the solar terminator (ST) moving over the northern America and physical mechanisms involved. Our observations confirm that the ST is a stable and repetitive source of ionospheric wave disturbances and the evidence of solar terminator generated disturbances has been demonstrated experimentally via the GPS TEC measurement. The main researches and results of this dissertation are as follows. First, the global traveling ionospheric disturbances (TIDs) during the drastic magnetic storms of October 29–31, 2003 were analyzed using the Global Position System (GPS) total electron content (TEC) data observed in the Asian-Australian, European and North American sectors. We collected the most comprehensive set of the TEC data from more than 900 GPS stations on the International GNSS Services (IGS) website and introduce here a strategy that combines polynomial fitting and multi-channel maximum entropy spectral analysis to obtain TID parameters. Moreover, in collaboration with my thesis advisor, I have developed an imaging technique of 2-dimensional map of TIDs structures to obtain spatial and temporal maps of large scale traveling ionospheric disturbances (LSTIDs). The clear structures of TEC perturbations map during the passage of TIDs were displayed. The results of our study are summarized as follows: (1) Large-scale TIDs (LSTIDs) and medium-scale TIDs (MSTIDs) were detected in all three sectors after the sudden commencement (SC) of the magnetic storm, and their features showed longitudinal and latitudinal dependences. The duration of TIDs was longer at higher latitudes than at middle latitudes, with a maximum of about 16 h. The TEC variation amplitude of LSTIDs was larger in the North American sector than in the two other sectors. At the lower latitudes, the ionospheric perturbations were more complicated, and their duration and amplitude were relatively longer and larger. (2) The periods and phase speeds of TIDs were different in these three sectors. In Europe, the TIDs propagated southward; in North America and Asia, the TIDs propagated southwestward; in the near-equator region, the disturbances propagated with the azimuth (the angle of the propagation direction of the LSTIDs measured clockwise from due north with 0°) of 210° showing the influence of Coriolis force; in the Southern Hemisphere, the LSTIDs propagated conjugatedly northwestward. Both the southwestward and northeastward propagating LSTIDs are found in the equatorial region. These results mean that the Coriolis effect cannot be ignored for the wave propagation of LSTIDs and that the propagation direction is correlated with the polar magnetic activity. (3) The day (day of year: 301) before the SC (sudden commencement) of magnetic storm, we observed a sudden TEC skip disturbances (±10 TECU). It should be a response for the high flux of proton during the solar flare event, but not the magnetic storms. Next, the most comprehensive and dense GPS network’s data from North-America region were used in this paper to analyze the medium scale traveling ionospheric disturbances (MSTIDs) which were generated by the moving solar terminator during the quiet days in 2005. We applied the multi-channel maximum entropy spectral analysis to calculated TID parameters, and found that the occurrence of ST-MSTIDs depends on the seasonal variations. The results of our study are summarized as follows: (1) MSTIDs stimulated by the moving ST (ST-MSTIDs) are detected at mid-latitudes after the passage of the solar terminator with the life time of 2～3 hours and the variation amplitude of 0.2～0.8 TECU. Spectral analysis indicated that the horizontal wavelength, average period, horizontal phase velocity of the MSTIDs are around 300±150 km，150±80 m/s and 25±15 min, respectively. In addition, ST-MSTIDs have wave fronts elongating the moving ST direction and almost parallel to ST. (2) The statistical results demonstrate that the dusk MSTIDs stimulated by ST is more obvious than the dawn MSTIDs in summer. On the contrary, the more-pronounced dawn MSTIDs occurs in winter. (3) Further analysis indicates that the seasonal variations of ST-MSTIDs occurrence frequency are most probably related to the seasonal differences of the variations of EUV flux in the ionosphere region and recombination process during sunrise and sunset period at mid-latitudes. Statistical study of occurrence characteristics of TIDs using the GPS network in North-American and European during solar maximum, In conclusion, statistical studies of the propagation characteristics of TIDs, which excited by the two common origins including geomagnetic storms and moving solar terminator, were involved with global GPS TEC databasein this thesis. We employed the multichannel maximum entropy spectral analysis method to diagnose the characteristics of propagation and evolvement of ionospheric disturbances, also, the characteristics of their regional distribution and climatological variations were revealed by the statistic analysis. The results of these studies can improve our knowledge about the energy transfer in the solar-terrestrial system and the coupling process between upper and lower atmosphere (thermosphere-ionosphere-mesosphere). On the other hand, our results of the investigation on TIDs generated by particular linear origin such as ST are important for developing ionospheric irregularity physics and modeling the transionosphere radio wave propagation. Besides, the GPS TEC representation of the ST-generated ionospheric structure suggests a better possibility for investigating this phenomenon. Subsequently, there are scientific meaning of the result of this dissertation to deeply discuss the energy transfer and coupling in the ionosphere, as well as realistic value to space weather forecast in the ionosphere region.

复杂介质中地震波场数值模拟研究——基于广义正交多项式迭积微分算子法

In this paper, we propose a new numerical modeling method – Convolutional Forsyte Polynomial Differentiator (CFPD), aimed at simulating seismic wave propagation in complex media with high efficiency and accuracy individually owned by short-scheme finite differentiator and general convolutional polynomial method. By adjusting the operator length and optimizing the operator coefficient, both global and local informations can be easily incorporated into the wavefield which is important to invert the undersurface geological structure. The key issue in this paper is to introduce the convolutional differentiator based on Forsyte generalized orthogonal polynomial in mathematics into the spatial differentiation of the first velocity-stress equation. To match the high accuracy of the spatial differentiator, this method in the time coordinate adopts staggered grid finite difference instead of conventional finite difference to model seismic wave propagation in heterogeneous media. To attenuate the reflection artifacts caused by artificial boundary, Perfectly Matched Layer (PML) absorbing boundary is also being considered in the method to deal with boundary problem due to its advantage of automatically handling large-angle emission. The PML formula for acoustic equation and first-order velocity-stress equation are also derived in this paper. There is little difference to implement the PML boundary condition in all kind of wave equations, but in Biot media, special attenuation factors should be taken. Numerical results demonstrate that the PML boundary condition is better than Cerjan absorbing boundary condition which makes it more suitable to hand the artificial boundary reflection. Based on the theories of anisotropy, Biot two-phase media and viscous-elasticity, this paper constructs the constitutive relationship for viscous-elastic and two-phase media, and further derives the first-order velocity-stress equation for 3D viscous-elastic and two-phase media. Numerical modeling using CFPD method is carried out in the above-mentioned media. The results modeled in the viscous-elastic media and the anisotropic pore elastic media can better explain wave phenomena of the true earth media, and can also prove that CFPD is a useful numerical tool to study the wave propagation in complex media.

扩散波勘探方法基础及其在石油油气藏开发中的应用研究

As we know, the essence of exploration is objective body determined by getting the information. Such as seismic、electrical and electromagnetic prospecting, they are the common methods of the exploration. Therefore, They have a complete set of theory now. In fact, the effective information can also be got by the diffusion way, it is called diffusion prospecting. The diffusion way prospecting is necessary and important. The way of diffusion prospecting is studied in the paper and main works include below: (1) On the basis of studying basic law of the diffusion, the paper gives the idea of diffusion wave and the formulas of computing diffusion wave function. (2) The paper studies the way of the diffusion prospecting and the methods of data processing. At the same time, it also expounds the characteristics and the applied foreground of the diffusion prospecting. (3) The paper gives the tomography idea and the basic method of diffusion CT. Meanwhile, it also expounds the foreground that the diffusion CT is applied in oil development prospecting. (4) As the inversion of the diffusion equation is a part of the diffusion prospecting way, the methods of diffusion equation inversion are studied and the two formulas are deduced --Laplace transform and polynomial fitting inversion formulas. As the other important result of diffusion equation inversion, the inversion can offer a new analysis method for well Testing in oil development. In order to show a set of methods in the paper feasible, forward、inversion and CT numerical simulation are done in the paper.

中国北方典型生态区表土孢粉组合特征及其与现代植被和气候的关系

Surface pollen assemblages and their relationhips with the modern vegetation and climate provide a foundation for investigating palaeo-environment conditions by fossil pollen analysis. A promising trend of palynology is to link pollen data more closely with ecology. In this study, I summarized the characteristics of surface pollen assemblages and their quantitative relation with the vegetation and climate of the typical ecological regions in northern China, based on surface pollen analysis of 205 sites and investigating of modern vegetation and climate. The primary conclusions are as follows:The differences in surface pollen assemblages for different vegetation regions are obvious. In the forest communities, the arboreal pollen percentages are more than 30%, herbs less than 50% and shrubs less than 10%; total pollen concentrations are more than 106 grains/g. In the steppe communities, arboreal pollen percentages are generally less than 5%; herb pollen percentages are more than 90%, and Artemisia and Chenopodiaceae are dominant in the pollen assemblages; total pollen concentrations range from 103 to 106 grains/g. In the desert communities, arboreal pollen percentages are less than 5%. Although Chenopodiaceae and Artemisia still dominate the pollen assemblages, Ephedra, Tamaricaceae and Nitraria are also significant important in the pollen assemblages; total pollen concentrations are mostly less than 104grains/g. In the sub-alpine or high and cold meadow communities, arboreal pollen percentages are less than 30%. and Cyperaceae is one of the most significant-taxa in the pollen assemblages. In the shrub communities, the pollen assemblages are consistent with the zonal vegetation; shrub pollen percentages are mostly less than 20%, except for Artemisia and Hippophae rhamnoides communities.There are obvious trends for the pollen percentage ratios of Artemisia to Chenopodiaceae (A/C), Pinus to Artemisia (P/A) and arbor to non-arbor (AP/NAP) in the different ecological regions. In the temperate deciduous broad-leaved forest region, the P/A ratios are generally higher than 0.1, the A/C ratios higher than 2 and the AP/NAP ratios higher than 0.3. In the temperate steppe regions, the P/A ratios are generally less than 0.1, the A/C ratios higher than 1 and the AP/NAP ratios less than 0.1. In the temperate desert regions, the P/A ratios are generally less than 0.1, the A/C ratios less than 1, and the AP/NAP ratios less than 0.1.The study on the representation and indication of pollen to vegetation shows that Pinus, Artemisia, Betula, Chenopodiaceae, Ephedra, Selaginella sinensis etc. are over-representative in the pollen assemblages and can only indicate the regional vegetation. Some pollen types, such as Quercus, Carpinus, Picea, Abies, Elaeagus, Larix, Salix, Pterocelis, Juglans, Ulmus, Gleditsia, Cotinus, Oleaceae, Spiraea, Corylus, Ostryopsis, Vites, Tetraena, Caragana, Tamaricaceae, Zygophyllum, Nitraria, Cyperaceae, Sanguisorba etc. are under-representative in the pollen assemblages, and can indicate the plant communities well. Populus, Rosaceae, Saxifranaceae, Gramineae, Leguminosae, Compositae, Caprifoliaceae etc. can not be used as significant indicators to the plants.The study on the relation of pollen percentages with plant covers shows that Pinus pollen percentages are more than 30% where pine trees exist in the surrounding region. The Picea+Abies pollen percentages are higher than 20% where the Picea+Abies trees are dominant in the communities, but less than 5% where the parent plants are sparse or absent. Larix pollen percentages vary from 5% to 20% where the Larix trees are dominant in the communities, but less than 5% where the parent plants are sparse or absent. Betula pollen percentages are higher than 40% where the Betula trees are dominant in the communities" but less than 5% where the parent plants are sparse or absent. Quercus pollen percentages are higher than 10% where the Quercus trees are dominant in the communities, but less than 1% where the parent plants sparse or absent. Carpinus pollen percentages vary from 5% to 15% where the Carpinus trees are dominant in the communities, but less than 1% where the parent plants are sparse or absent. Populus pollen percentages are about 0-5% at pure Populus communities, but cannot be recorded easily where the Populus plants mixed with other trees in the communities. Juglans pollen accounts for 25% to 35% in the forest of Juglans mandshurica, but less than 1% where the parent plants are sparse or absent. Pterocelis pollen percentages are less than 15% where the Pterocelis trees are dominant in the communities, but cannot be recorded easily where the parent plants are sparse or absent. Ulmus pollen percentages are more than 8% at Ulmus communities, but less than 1% where the Ulmus plants mixed with other trees in the communities. Vitex pollen percentages increase along with increasing of parent plant covers, but the maximum values are less than 10 %. Caragana pollen percentages are less than 20 % where the Caragana plant are dominant in the communities, and cannot be recorded easily where the parent plants are sparse or absent. Spiraea pollen percentages are less than 16 % where the Spiraea plant are dominant in the communities, and cannot be recorded easily where the parent plants are sparse or absent.The study on the relation of surface pollen assemblages with the modern climate shows that, in the axis 1 of DCA, surface samples scores have significant correlation with the average annual precipitations, and the highest determination coefficient (R2) is 0.8 for the fitting result of the third degree polynomial functions. In the axis 2 of DCA, the samples scores have significant correlation with the average annual temperatures, average July temperatures and average January temperatures, and the determination coefficient falls in 0.13-0.29 for the fitting result of the third degree polynomial functions with the highest determination coefficient for the average July temperature.The sensitivity of the different pollen taxa to climate change shows that some pollen taxa such as Pinus, Quercus, Carpinus, Juglans, Spiraea, Oleaceae, Gramineae, Tamariaceae and Ephedra are only sensitive to the change in precipitation.

Study on thermodynamic properties of polypyromellitimide molding powder

Polypyromellitimide molding powder has been prepared. In the 78-370 K range, the dependence of the specific heat capacity (c(p)) on the temperature (T) is given by the polynomial: c(p)=0.8163+0.4592X+0.02468X(2)+0.1192X(3)+0.05659X(4) (J K-1 g(-1)) where X=(T-225.5)/144.5. Thermal decomposition in air starts at 716 K, and is complete at 1034 K. The standard combustion enthalpy is Delta(c)H=-26.442 kJ g(-1). (C) 2000 Elsevier Science B.V. All rights reserved.

Low temperature heat capacity of (S)-ibuprofen

Molar heat capacities of ( S)-ibuprofen were precisely measured with a small sample precision automated adiabatic calorimeter over the temperature range from 80 to 370 K. Experimental heat capacities were fitted into a polynomial equation of heat capacities ( C-p,C- m) with reduced temperature ( X), [ X = f(T)]. The polynomial equations for ( S)-ibuprofen were C-p,C- m(s) = - 39.483 X-4 - 66. 649 X-3 + 95. 196 X-2 + 210. 84 X + 172. 98 in solid state and C-p,C- m(L) = 7. 191X(3) + 4. 2774 X-2 + 56. 365 X + 498. 5 in liquid state. The thermodynamic functions relative to the reference temperature of 298. 15 K, H-T - H-298.15 and S-T - S-298.15, were derived for the( S)-ibuprofen. A fusion transition at T-m = (324. 15 +/- 0. 02) K was found from the C-p - T curve. The molar enthalpy and entropy of the fusion transition were determined to be (18. 05 +/- 0. 31) kJ.mol(-1) and (55. 71 +/- 0. 95) J.mol(-1).K-1, respectively. The purity of the ( S)-ibuprofen was determined to be 99. 44% on the basis of the heat capacity measurement. Finally, the heat capacities of ( S)-ibuprofen and racemic ibuprofen were compared.

Thermodynamic properties of ethanol and gasoline blended fuel

The heat capacities (C-p) of five types of gasohol (50 wt % ethanol and 50 wt % unleaded gasoline 93(#) (E50), 60 wt % ethanol and 40 wt % unleaded gasoline 93(#) (E60), 70 wt % ethanol and 30 wt % unleaded gasoline 93(#) (E70), 80 wt % ethanol and 20 wt % unleaded gasoline 93(#) (E80), and 90 wt % ethanol and 10 wt % unleaded gasoline 93(#) (E90), where the "93" denotes the octane number) were measured by adiabatic calorimetry in the temperature range of 78-320 K. A glass transition was observed at 95.61, 96.14, 96.56, 96.84, and 97.08 K for samples from the E50, E60, E70, E80, and E90 systems, respectively. A liquid-solid phase transition and a solid-liquid phase transition were observed in the respective temperature ranges of 118-153 and 155-163 K for E50, 117-150 and 151-164 K for E60, 115-154 and 154-166 K for E70, 113-152 and 152-167 K for E80, and 112-151 and 1581-167 K for E90. The polynomial equations of Cp and the excess heat capacities (C-p(E)), with respect to the thermodynamic temperature, were established through least-squares fitting. Based on the thermodynamic relationship and the equations obtained, the thermodynamic functions and the excess thermodynamic functions of the five gasohol samples were derived.