延安地区全新世黄土孢粉分析及其植被演化

采用孢粉分析和 1 4 C测年法 ,通过花粉植被化模拟技术 ,研究了延安地区全新世的孢粉组合特征及其古植被的演替历史。全新世时期延安地区的气候变化在全球和中国区域气候的格局下 ,也呈现出早、中、晚三个气候阶段 ,中期 (75 0 0～ 45 0 0 a)为气候最佳期 ,与之相对应的植被演替也分为三个阶段 :中期植被为松、桦、栎、漆为主具亚热带成分的暖温带落叶阔叶林 ,早期和晚期的气候较干冷 ,乔木花粉的含量减少 ,呈现以松、栎为主的暖温带针阔叶混交林

基于主动建模的无人直升机增强LQR控制

为了解决无人直升机控制问题,通过把主动建模与LQR(Linear Quadratic Regulator)控制相结合,提出一种能补偿模型差的控制方法。该方法在悬停状态下,采用简化模型设计LQR控制器,并通过UKF(Un-scented-Kalman-Filter)在线估计简化模型与全状态模型的模型差,使用模型差作为补偿项对LQR控制增强。针对实际直升机动力学模型进行仿真,验证了基于UKF的估计和增强LQR控制的有效性。仿真实验结果证明,基于UKF的主动建模技术能够快速估计状态和参数变化,并且增强LQR控制能够使系统适应模型不确定性。

地震波模拟：非规则网格离散的数值解法、吸收边界及其实际应用

This dissertation presents a series of irregular-grid based numerical technique for modeling seismic wave propagation in heterogeneous media. The study involves the generation of the irregular numerical mesh corresponding to the irregular grid scheme, the discretized version of motion equations under the unstructured mesh, and irregular-grid absorbing boundary conditions. The resulting numerical technique has been used in generating the synthetic data sets on the realistic complex geologic models that can examine the migration schemes. The motion equation discretization and modeling are based on Grid Method. The key idea is to use the integral equilibrium principle to replace the operator at each grid in Finite Difference scheme and variational formulation in Finite Element Method. The irregular grids of complex geologic model is generated by the Paving Method, which allow varying grid spacing according to meshing constraints. The grids have great quality at domain boundaries and contain equal quantities of nodes at interfaces, which avoids the interpolation of parameters and variables. The irregular grid absorbing boundary conditions is developed by extending the Perfectly Matched Layer method to the rotated local coordinates. The splitted PML equations of the first-order system is derived by using integral equilibrium principle. The proposed scheme can build PML boundary of arbitrary geometry in the computational domain, avoiding the special treatment at corners in a standard PML method and saving considerable memory and computation cost. The numerical implementation demonstrates the desired qualities of irregular grid based modeling technique. In particular, (1) smaller memory requirements and computational time are needed by changing the grid spacing according to local velocity; (2) Arbitrary surfaces and interface topographies are described accurately, thus removing the artificial reflection resulting from the stair approximation of the curved or dipping interfaces; (3) computational domain is significantly reduced by flexibly building the curved artificial boundaries using the irregular-grid absorbing boundary conditions. The proposed irregular grid approach is apply to reverse time migration as the extrapolation algorithm. It can discretize the smoothed velocity model by irregular grid of variable scale, which contributes to reduce the computation cost. The topography. It can also handle data set of arbitrary topography and no field correction is needed.

电离层参量M(3000)F2和hmF2的经验正交函数方法建模研究

The ionospheric parameter M(3000)F2 (the so-called transmission factor or the propagation factor) is important not only in practical applications such as frequency planning for radio-communication but also in ionospheric modeling. This parameter is strongly anti-correlated with the ionospheric F2-layer peak height hmF2，a parameter often used as a key anchor point in some widely used empirical models of the ionospheric electron density profile (e.g., in IRI and NeQuick models). Since hmF2 is not easy to obtain from measurements and M(3000)F2 can be routinely scaled from ionograms recorded by ionosonde/digisonde stations distributed globally and its data has been accumulated for a long history, usually the value of hmF2 is calculated from M(3000)F2 using the empirical formula connecting them. In practice, CCIR M(3000)F2 model is widely used to obtain M(3000)F2 value. However, recently some authors found that the CCIR M(3000)F2 model has remarkable discrepancies with the measured M(3000)F2, especially in low-latitude and equatorial regions. For this reason, the International Reference Ionosphere (IRI) research community proposes to improve or update the currently used CCIR M(3000)F2 model. Any efforts toward the improvement and updating of the current M(3000)F2 model or newly development of a global hmF2 model are encouraged. In this dissertation, an effort is made to construct the empirical models of M(3000)F2 and hmF2 based on the empirical orthogonal function (EOF) analysis combined with regression analysis method. The main results are as follows: 1. A single station model is constructed using monthly median hourly values of M(3000)F2 data observed at Wuhan Ionospheric Observatory during the years of 1957–1991 and compared with the IRI model. The result shows that EOF method is possible to use only a few orders of EOF components to represent most of the variance of the original data set. It is a powerful method for ionospheric modeling. 2. Using the values of M(3000)F2 observed by ionosondes distributed globally, data at grids uniformly distributed globally were obtained by using the Kriging interpolation method. Then the gridded data were decomposed into EOF components using two different coordinates: (1) geographical longitude and latitude; (2) modified dip (Modip) and local time. Based on the EOF decompositions of the gridded data under these two coordinates systems, two types of the global M(3000)F2 model are constructed. Statistical analysis showed that the two types of the constructed M(3000)F2 model have better agreement with the observational M(3000)F2 than the M(3000)F2 model currently used by IRI. The constructed models can represent the global variations of M(3000)F2 better. 3. The hmF2 data used to construct the hmF2 model were converted from the observed M(3000)F2 based on the empirical formula connecting them. We also constructed two types of the global hmF2 model using the similar method of modeling M(3000)F2. Statistical analysis showed that the prediction of our models is more accurate than the model of IRI. This demonstrated that using EOF analysis method to construct global model of hmF2 directly is feasible. The results in this thesis indicate that the modeling technique based on EOF expansion combined with regression analysis is very promising when used to construct the global models of M(3000)F2 and hmF2. It is worthwhile to investigate further and has the potential to be used to the global modeling of other ionospheric parameters.

准噶尔盆地南缘台13井区侏罗系油藏表征与建模

Junggar Basin has a large amount of recoverable reserves, However, due to the unfavorable factors, such as bad seismic data quality, complex structure with many faults and less wells, the exploration of oil and gas is still relatively limited, so advanced theoretical guidance and effective technical supports are desirable. Based on the theories of sedimentology, as well as comprehensive studies of outcrops, seismic data, drilling data and setting of this area, the paper establishes the isochronous correlation framework, and analyzes the sedimentary facies types and provenance direction, and obtains the profile and plain maps of the sedimentary facies combined with the logging constrained inversion. Then the paper analyzes the reservoir controlling factors, reservoir lithology attribute, 4-property relationship and sensibility based on the sedimentary facies research, and sets up a 3D geological model using facies controlled modeling. Finally, the paper optimizes some target areas with the conclusions of reservoir, structure and reservoir formation.Firstly, the paper establishs the isochronous correlation framework by the seismic data, drilling data and setting of this area. The sedimentary facies in Tai13 well block are braided river and meandering river according to the analysis of the lithology attribute, logging facies and sedimentary structure attribute of outcrop. The concept of “wetland” is put forward for the first time. The provenance direction of Badaowan and Qigu formation is obtained by the geology setting, sedimentary setting and paleocurrent direction. The paper obtains the profile and plain maps of the sedimentary facies from the sand value of the wells and the sand thickness maps from the logging constrained inversion. Then, this paper takes characteristics and control factors of the Jurassic reservoirs analysis on thin section observation, scanning transmission electron microscope observation and find out the petrology characteristics of reservoir, space types of reservoir and lithofacies division. In this area, primary pores dominate in the reservoir pores, which believed that sedimentation played the most important roles of the reservoir quality and diagenesis is the minor factor influencing secondary porosity. Using stochastic modeling technique，the paper builds quantitative 3-D reservoir Parameter. Finally, combined the study of structure and reservoir formation, the reservoir distribution regularity is concluded: (a) structures control the reservoir formation and accumulation. (b) Locating in the favorable sedimentary facies belt. And the area which meets these conditions mentioned above is a good destination for exploration.

地质体科学计算可视化研究及其应用

Describing visually space-time properties of geological phenomena consists of one of the most important parts in geology research. Such visual images are of usually helpful for analyzing geological phenomena and for discovering the regulations behind geological phenomena. This report studies mainly three application problems of scientific visualization in geology: (Dvisualizing geological body A new geometric modeling technique with trimmed surface patches has been eveloped to visualize geological body. Constructional surfaces are represented as trimmed surfaces and a constructional solid is represented by the upper and lower surface composed of trimmed surface patches from constructional surfaces. The technique can completely and definitely represent the structure of geological body. It has been applied in visualization for the coal deposit in Huolinhe, the aquifer thermal energy storage in Tianjin and the structure of meteorite impact in Cangshan et al. (2)visualizing geological space field Efficient visualization methods have been discussed. Marching-Cube algorithm used has been improved and is used to extract iso~surface from 3D data set, iso-line from 2D data set and iso-point from ID data set. The improved method has been used to visualize distribution and evolution of the abnormal pressures in Zhungaer Basin. (3)visualizing porous space a novel way was proposed to define distance from any point to a convex set. Thus a convex set skeleton-based implicit surface modeling technique is developed and used to construct a simplified porous space model. A Buoyancy Percolation numerical simulation platform has been developed to simulate the process of migration of oil in the porous media saturated with water.

Modeling Ammonothermal Growth of GaN Single Crystals: The Role of Transport

Single crystal gallium nitride (GaN) is an important technological material used primarily for the manufacture of blue light lasers. An important area of contemporary research is developing a viable growth technique. The ammonothermal technique is an important candidate among many others with promise of commercially viable growth rates and material quality. The GaN growth rates are a complicated function of dissolution kinetics, transport by thermal convection and crystallization kinetics. A complete modeling effort for the growth would involve modeling each of these phenomena and also the coupling between these. As a first step, the crystallization and dissolution kinetics were idealized and the growth rates as determined purely by transport were investigated. The growth rates thus obtained were termed ‘transport determined growth rates’ and in principle are the maximum growth rates that can be obtained for a given configuration of the system. Using this concept, a parametric study was conducted primarily on the geometric and the thermal boundary conditions of the system to optimize the ‘transport determined growth rate’ and determine conditions when transport might be a bottleneck.

Modeling, numerical methods, and simulation for particle-fluid two-phase flow problems

In this paper, we study the issues of modeling, numerical methods, and simulation with comparison to experimental data for the particle-fluid two-phase flow problem involving a solid-liquid mixed medium. The physical situation being considered is a pulsed liquid fluidized bed. The mathematical model is based on the assumption of one-dimensional flows, incompressible in both particle and fluid phases, equal particle diameters, and the wall friction force on both phases being ignored. The model consists of a set of coupled differential equations describing the conservation of mass and momentum in both phases with coupling and interaction between the two phases. We demonstrate conditions under which the system is either mathematically well posed or ill posed. We consider the general model with additional physical viscosities and/or additional virtual mass forces, both of which stabilize the system. Two numerical methods, one of them is first-order accurate and the other fifth-order accurate, are used to solve the models. A change of variable technique effectively handles the changing domain and boundary conditions. The numerical methods are demonstrated to be stable and convergent through careful numerical experiments. Simulation results for realistic pulsed liquid fluidized bed are provided and compared with experimental data. (C) 2004 Elsevier Ltd. All rights reserved.

A simple method to simulate shrinkage-induced cracking in cement-based composites by lattice-type modeling

A new numerical procedure is proposed to investigate cracking behaviors induced by mismatch between the matrix phase and aggregates due to matrix shrinkage in cement-based composites. This kind of failure processes is simplified in this investigation as a purely spontaneous mechanical problem, therefore, one main difficulty during simulating the phenomenon lies that no explicit external load serves as the drive to propel development of this physical process. As a result, it is different from classical mechanical problems and seems hard to be solved by using directly the classical finite element method (FEM), a typical kind of "load -> medium -> response" procedures. As a solution, the actual mismatch deformation field is decomposed into two virtual fields, both of which can be obtained by the classical FEM. Then the actual response is obtained by adding together the two virtual displacement fields based on the principle of superposition. Then, critical elements are detected successively by the event-by-event technique. The micro-structure of composites is implemented by employing the generalized beam (GB) lattice model. Numerical examples are given to show the effectiveness of the method, and detailed discussions are conducted on influences of material properties.

Physical Modeling And Parametric Study On Two-Degree-Of-Freedom Viv Of A Cylinder Near Rigid Wall

Unlike most previous studies on the transverse vortex-induced vibration(VIV) of a cylinder mainly under the wallfree condition (Williamson & Govardhan,2004),this paper experimentally investigates the vortex-induced vibration of a cylinder with two degrees of freedom near a rigid wall exposed to steady flow.The amplitude and frequency responses of the cylinder are discussed.The lee wake flow patterns of the cylinder undergoing VIV were visualized by employing the hydrogen bubble technique.The effects of the gap-to-diameter ratio (e0/D) and the mass ratio on the vibration amplitude and frequency are analyzed.Comparisons of VIV response of the cylinder are made between one degree (only transverse) and two degrees of freedom (streamwise and transverse) and those between the present study and previous ones.The experimental observation indicates that there are two types of streamwise vibration,i.e.the first streamwise vibration (FSV) with small amplitude and the second streamwise vibration (SSV) which coexists with transverse vibration.The vortex shedding pattem for the FSV is approximately symmetric and that for the SSV is alternate.The first streamwise vibration tends to disappear with the decrease of e0/D.For the case of large gap-to-diameter ratios (e.g.e0/D = 0.54～1.58),the maximum amplitudes of the second streamwise vibration and transverse one increase with the increasing gapto-diameter ratio.But for the case of small gap-to-diameter ratios (e.g.e0/D = 0.16,0.23),the vibration amplitude of the cylinder increases slowly at the initial stage (i.e.at small reduced velocity V,),and across the maximum amplitude it decreases quickly at the last stage (i.e.at large Vr).Within the range ofthe examined small mass ratio (m＜4),both streamwise and transverse vibration amplitude of the cylinder decrease with the increase of mass ratio for the fixed value of V,.The vibration range (in terms of Vr ) tends to widen with the decrease of the mass ratio.In the second streamwise vibration region,the vibration frequency of the cylinder with a small mass ratio (e.g.mx = 1.44) undergoes a jump at a certain Vr,.The maximum amplitudes of the transverse vibration for two-degree-of-freedom case is larger than that for one-degree-of-freedom case,but the transverse vibration frequency of the cylinder with two degrees of freedom is lower than that with one degree of freedom (transverse).

Analysis of mode characteristics for deformed square resonators by FDTD technique

The mode frequencies and quality factors (Q-factors) in two-dimensional (2-D) deformed square resonators are analyzed by finite-difference time-domain (FDTD) technique. The results show that the deformed square cavities with circular and cut corners have larger Q-factors than the perfect ones at certain conditions. For a square cavity with side length of 2 mu m and refractive index of 3.2, the mode Q-factor can increase 13 times as the perfect corners are replaced by a quarter of circle with radius of 0.3 pm. Furthermore the blue shift with the increasing deformations is found as a result of the reduction in effective resonator area. In square cavities with periodic roughness at sidewalls which maintains the symmetry of the square, the Q-factors of the whisperin gallery (WG)-like modes are still one order of magnitude larger that those of non-WG-like modes. However, the Q-tactors of these two types of modes are of the same order in the square cavity with random roughness. We also find that the rectangular and rhombic deformation largely reduce the Q-factors with the increasing offset and cause the splitting of the doubly degenerate modes due to the breaking of certain symmetry properties.

Fast monostatic GPR modeling

We propose the exploding-reflector method to simulate a monostatic survey with a single simulation. The exploding reflector, used in seismic modeling, is adapted for ground-penetrating radar (GPR) modeling by using the analogy between acoustic and electromagnetic waves. The method can be used with ray tracing to obtain the location of the interfaces and estimate the properties of the medium on the basis of the traveltimes and reflection amplitudes. In particular, these can provide a better estimation of the conductivity and geometrical details. The modeling methodology is complemented with the use of the plane-wave method. The technique is illustrated with GPR data from an excavated tomb of the nineteenth century.

Modeling gross primary production of alpine ecosystems in the Tibetan Plateau using MODIS images and climate data

The eddy covariance technique provides measurements of net ecosystem exchange (NEE) Of CO2 between the atmosphere and terrestrial ecosystems, which is widely used to estimate ecosystem respiration and gross primary production (GPP) at a number Of CO2 eddy flux tower sites. In this paper, canopy-level maximum light use efficiency, a key parameter in the satellite-based Vegetation Photosynthesis Model (VPM), was estimated by using the observed CO2 flux data and photosynthetically active radiation (PAR) data from eddy flux tower sites in an alpine swamp ecosystem, an alpine shrub ecosystem and an alpine meadow ecosystem in Qinghai-Tibetan Plateau, China. The VPM model uses two improved vegetation indices (Enhanced Vegetation Index (EVI), Land Surface Water Index (LSWI)) derived from the Moderate Resolution Imaging Spectral radiometer (MODIS) data and climate data at the flux tower sites, and estimated the seasonal dynamics of GPP of the three alpine grassland ecosystems in Qinghai-Tibetan Plateau. The seasonal dynamics of GPP predicted by the VPM model agreed well with estimated GPP from eddy flux towers. These results demonstrated the potential of the satellite-driven VPM model for scaling-up GPP of alpine grassland ecosystems, a key component for the study of the carbon cycle at regional and global scales. (c) 2006 Elsevier Inc. All rights reserved.

On the evaluation of diffusivities in gels using the diffusion cell technique

Modeling of the gel-immobilized cell system requires accurate measurement of diffusion coefficients. Three methods of the quasi-steady-state (QSS) method, the time-lag (TL) method and a variant quasi-steady-state (VQSS) method were critically assessed and compared for the evaluation of diffusivities using the diffusion cell technique. Experimental data from our laboratory were used for the analysis of the influence of crucial theoretical assumptions not being fulfilled in each method. The results highlighted a risk in obtaining highly variable diffusion coefficients by not validating the QSS and the accuracy of the measurements. In the TL method, the estimation of diffusivities based on the plot intercept that was mostly used in the literature, results in a many fold lower value when compared to that based on the plot slope. The comparison with the QSS and VQSS methods confirmed similar diffusivity obtained by the TL method based on the plot slope. It thus suggested that the correct estimation of diffusivities by the TL method could be based on the plot slope only. Furthermore, the errors associated with the solute mass in the gel, the sample withdrawal and the non-negligible concentration changes in the chambers were also discussed. It is concluded that diffusion cell technique has to be employed cautiously for a correct evaluation of diffusivities. (C) 2001 Elsevier Science B.V. All rights reserved.