Shock propagation and MPT noise from a transonic rotor in non-uniform flow

One of the major challenges in hig4h-speed fan stages used in compact, embedded propulsion systems is inlet distortion noise. A body-force-based approach for the prediction of multiple-pure-tone (MPT) noise was previously introduced and validated. In this paper, it is employed with the objective of quantifying the effects of non-uniform flow on the generation and propagation of MPT noise. First-of-their-kind back-to-back coupled aero-acoustic computations were carried out using the new approach for conventional and serpentine inlets. Both inlets delivered flow to the same NASA/GE R4 fan rotor at equal corrected mass flow rates. Although the source strength at the fan is increased by 45 dB in sound power level due to the non-uniform inflow, farfield noise for the serpentine inlet duct is increased on average by only 3.1 dBA overall sound pressure level in the forward arc. This is due to the redistribution of acoustic energy to frequencies below 11 times the shaft frequency and the apparent cut-off of tones at higher frequencies including blade-passing tones. The circumferential extent of the inlet swirl distortion at the fan was found to be 2 blade pitches, or 1/11th of the circumference, suggesting a relationship between the circumferential extent of the inlet distortion and the apparent cut-off frequency perceived in the far field. A first-principles-based model of the generation of shock waves from a transonic rotor in non-uniform flow showed that the effects of non-uniform flow on acoustic wave propagation, which cannot be captured by the simplified model, are more dominant than those of inlet flow distortion on source noise. It demonstrated that non-linear, coupled aerodynamic and aeroacoustic computations, such as those presented in this paper, are necessary to assess the propagation through non-uniform mean flow. A parametric study of serpentine inlet designs is underway to quantify these propagation effects. Copyright © 2011 by ASME.

Gaussian Approximation Potential: an interatomic potential derived from first principles Quantum Mechanics

Simulation of materials at the atomistic level is an important tool in studying microscopic structure and processes. The atomic interactions necessary for the simulation are correctly described by Quantum Mechanics. However, the computational resources required to solve the quantum mechanical equations limits the use of Quantum Mechanics at most to a few hundreds of atoms and only to a small fraction of the available configurational space. This thesis presents the results of my research on the development of a new interatomic potential generation scheme, which we refer to as Gaussian Approximation Potentials. In our framework, the quantum mechanical potential energy surface is interpolated between a set of predetermined values at different points in atomic configurational space by a non-linear, non-parametric regression method, the Gaussian Process. To perform the fitting, we represent the atomic environments by the bispectrum, which is invariant to permutations of the atoms in the neighbourhood and to global rotations. The result is a general scheme, that allows one to generate interatomic potentials based on arbitrary quantum mechanical data. We built a series of Gaussian Approximation Potentials using data obtained from Density Functional Theory and tested the capabilities of the method. We showed that our models reproduce the quantum mechanical potential energy surface remarkably well for the group IV semiconductors, iron and gallium nitride. Our potentials, while maintaining quantum mechanical accuracy, are several orders of magnitude faster than Quantum Mechanical methods.

Strategies for modeling turbulent flows in electronics

Hybrid methods based on the Reynolds Averaged Navier Stokes (RANS) equations and the Large Eddy Simulation (LES) formulation are investigated to try and improve the accuracy of heat transfer and surface temperature predictions for electronics systems and components. Two relatively low Reynolds number flows are studied using hybrid RANS-LES, RANS-Implicit-LES (RANS-ILES) and non-linear LES models. Predictions using these methods are in good agreement with each other, even using different grid resolutions. © 2008 IEEE.

Predictive techniques for turbulent oscillatory flows

The prediction of turbulent oscillatory flow at around transitional Reynolds numbers is considered for an idealized electronics system. To assess the accuracy of turbulence models, comparison is made with measurements. A stochastic procedure is used to recover instantaneous velocity time traces from predictions. This procedure enables more direct comparison with turbulence intensity measurements which have not been filtered to remove the oscillatory flow component. Normal wall distances, required in some turbulence models, are evaluated using a modified Poisson equation based technique. A range of zero, one and two equation turbulence models are tested, including zonal and a non-linear eddy viscosity models. The non-linear and zonal models showed potential for accuracy improvements.

Vibration of a high-altitude tethered balloon with application to climate engineering

A balloon tethered at an altitude of 20 km could deliver a particulate cloud leading to global cooling. Tethering a balloon at this altitude poses significant problems with respect to vibration and stability, especially in regions of high wind. No-one has ever proposed, yet alone launched, a balloon at an altitude of 20 km tethered to the ground. Owing to wind, the tether needs to be 23 km in length and is to be fixed to a ship at sea or on land in equatorial regions. Whilst the balloon at 20 km is subject to relatively modest wind conditions, at jet stream altitudes (10km) the tether will experience much higher wind loadings, not only because of the high wind speeds of up to 300 km / hr but also because of the high air density. A tether of circular cross section in these high winds will be subject to horizontal and downward drag forces that would bring the aerostat down. For this reason it is advantageous to consider a self-aligning tether of an aerodynamic cross section whereby it is possible to reduce the drag substantially. One disadvantage of a non-circular tether is the possibility of flutter and galloping instabilities. It is reasonably straightforward to model these phenomena for short lengths of aerofoil, but the situation becomes more complex for a 20 km tensioned tether with large deflection and curvature, variable wind speed, variable air density and variable tension. Analysis using models of infinite length are used to establish the stability at a local scale where the tension, aerodynamic and geometric properties are considered constant. Dispersion curve analysis is useful here. But for dynamics on a long-wavelength scale (several km) then a full non-linear analysis is required. This non-linear model can be used to establish the local values of tension appropriate for the dispersion analysis. This keynote presentation will give some insight into these issues.

An efficient approach of modelling the flexural cracking behaviour of un-notched plain concrete prisms subject to monotonic and cyclic loading

In the field of vibration-based damage detection of concrete structures efficient damage models are needed to better understand changes in the vibration properties of cracked structures. These models should quantitatively replicate the damage mechanisms in concrete and easily be used as damage detection tools. In this paper, the flexural cracking behaviour of plain concrete prisms subject to monotonic and cyclic loading regimes under displacement control is tested experimentally and modelled numerically. Four-point bending tests on simply supported un-notched prisms are conducted, where the cracking process is monitored using a digital image correlation system. A numerical model, with a single crack at midspan, is presented where the cracked zone is modelled using the fictitious crack approach and parts outside that zone are treated in a linear-elastic manner. The model considers crack initiation, growth and closure by adopting cyclic constitutive laws. A multi-variate Newton-Raphson iterative solver is used to solve the non-linear equations to ensure equilibrium and compatibility at the interface of the cracked zone. The numerical results agree well with the experiments for both loading scenarios. The model shows good predictions of the degradation of stiffness with increasing load. It also approximates the crack-mouth-opening-displacement when compared with the experimental data of the digital image correlation system. The model is found to be computationally efficient as it runs full analysis for cyclic loading in less than 2. min, and it can therefore be used within the damage detection process. © 2013 Elsevier Ltd.

Assessing boost-assist options for turbocharged engines using 1-D engine simulation and model predictive control

Delivering acceptable low end torque and good transient response is a significant challenge for all turbocharged engines. As downsized gasoline engines and Diesel engines make up a larger and larger proportion of the light-duty engines entering the market, the issue takes on greater significance. Several schemes have been proposed to improve torque response in highly boosted engines, including the use of electrical assist turbochargers and compressed air assist. In this paper we examine these methods with respect to their effectiveness in improving transient response and their relative performance along with some of the practical considerations for real world application. Results shown in this paper are from 1-D simulations using the Ricardo WAVE software package. The simulation model is based on a production light-duty Diesel engine modified to allow the introduction of compressed air at various points in the air-path as well as direct torque application to the turbocharger shaft (such as might be available from an electrical assist turbocharger). Whilst the 1-D simulation software provides a suitable environment for investigating the various boost assistance options, the overall air path performance also depends upon the control system. The introduction of boost assistance complicates the control in two significant ways: the system may run into constraints (such as compressor surge) that are not encountered in normal operation and the assistance introduces an additional control input. Production engine controllers are usually based on gain-scheduled PID control and extensive calibration. For this study, the non-linear nature of the engine together with the multiple configurations considered and the slower than real-time execution of 1-D models makes such an approach time consuming. Moreover, an ad-hoc approach would leave some doubt as to the fairness of comparisons between the different boost-assist options. Model Predictive Control has been shown to offer a convenient approach to controlling the 1-D simulations in a close to optimal manner for a typical Diesel VGT-EGR air path configuration. We show that the same technique can be applied to all the considered assistance methods with only modest calibration effort required. Copyright © 2012 SAE International.

First principles interatomic potential for tungsten based on Gaussian process regression

An accurate description of atomic interactions, such as that provided by first principles quantum mechanics, is fundamental to realistic prediction of the properties that govern plasticity, fracture or crack propagation in metals. However, the computational complexity associated with modern schemes explicitly based on quantum mechanics limits their applications to systems of a few hundreds of atoms at most. This thesis investigates the application of the Gaussian Approximation Potential (GAP) scheme to atomistic modelling of tungsten - a bcc transition metal which exhibits a brittle-to-ductile transition and whose plasticity behaviour is controlled by the properties of $\frac{1}{2} \langle 111 \rangle$ screw dislocations. We apply Gaussian process regression to interpolate the quantum-mechanical (QM) potential energy surface from a set of points in atomic configuration space. Our training data is based on QM information that is computed directly using density functional theory (DFT). To perform the fitting, we represent atomic environments using a set of rotationally, permutationally and reflection invariant parameters which act as the independent variables in our equations of non-parametric, non-linear regression. We develop a protocol for generating GAP models capable of describing lattice defects in metals by building a series of interatomic potentials for tungsten. We then demonstrate that a GAP potential based on a Smooth Overlap of Atomic Positions (SOAP) covariance function provides a description of the $\frac{1}{2} \langle 111 \rangle$ screw dislocation that is in agreement with the DFT model. We use this potential to simulate the mobility of $\frac{1}{2} \langle 111 \rangle$ screw dislocations by computing the Peierls barrier and model dislocation-vacancy interactions to QM accuracy in a system containing more than 100,000 atoms.

中国卧龙自然保护区不同海拔川滇高山栎(Quercus aquifolioides)群体的遗传变异

海拔梯度造成的环境异质性，如崎岖的地形、复杂的植被结构以及花期延迟等可能会极大地影响到物种的形态和遗传变异格局。理解物种形态和遗传变异的海拔格局对于物种多样性的管理和保护是非常重要的。尽管植物群体遗传学是一个飞速发展的研究领域，然而与海拔相关的形态变异、遗传变异及群体间遗传差异的研究却很少。到目前为止，还不清楚遗传变异与海拔之间是否必然的相关性。 川滇高山栎是一种重要的生态和经济型树种，广泛分布于中国西南的四川、西藏、贵州和云南省的高海拔地区，在保持水土、调节气候方面起着十分重要的作用。尽管主要受阳光限制而仅分布于阳坡，但其海拔梯度范围较大，表明川滇高山栎对不同的环境具有很强的适应性。本文通过叶型及生理响应、微卫星分子标记和扩增性片段长度多态性方法，试图探索川滇高山栎叶沿海拔梯度的形态和生理响应及其沿海拔梯度的遗传变异格局，为川滇高山栎的保护和利用提供进一步的遗传学理论依据和技术指导。 对叶形、含氮量及碳同位素的试验结果表明，平均比叶面积、气孔密度、气孔长度和气孔指数等气孔参数随海拔的升高呈非线性变化。在海拔大于2800 m时，川滇高山栎的比叶面积、气孔长度和气孔指数都随海拔升高而降低，但是在海拔小于2800 m时，这些指标都随海拔的升高而增大。相对而言，单位叶面积的含氮量和碳同位素则表现出相反的变化模式。另外，比叶面积是决定碳同位素沿海拔梯度变化的最重要参数。本研究结果表明，海拔2800 m附近是川滇高山栎生长和发育的最适地带，在这里生长的植物叶片厚度更薄、气孔更大、叶碳同位素值更小。 利用六对微卫星引物对五个不同海拔川滇高山栎群体遗传多样性进行研究，结果表明，群体内表现出较高的遗传多样性，平均每位点等位基因数11.33个，平均期望杂合度达0.820。群体间差异较小，分化仅为6.6％。聚类分析也并没有显示出明显的海拔格局。然而低频率等位基因却与海拔呈显著性正相关(R2=0.97, P < 0.01)，表明在高海拔处，川滇高山栎以更多的稀有基因来适应恶劣的环境条件。本试验结果表明由海拔梯度形成的选择性压力对川滇高山栎群体的遗传变异影响并不明显。 为了进一步探讨川滇高山栎群体遗传变异与海拔之间的相互关系，我们还对其进行了扩增性片段长度多态性分析。结果表明：(1)随海拔的升高(从群体WL2到群体WL5)，群体内遗传变异降低，而群体间遗传差异增加；(2)低海拔群体WL1表现出最低的遗传变异性(HE = 0.181)，同时与其余四个群体间呈现出最大的遗传差异性(平均FST = 0.0596)；(3)在除去低海拔群体WL1后，Mantel检测表明群体间遗传距离与海拔距离之间表现出正相关性。另外，研究结果还表明，遗传变异受生境条件(过度的湿热环境)及人为干扰(火烧、砍伐和放牧)的影响，这一点至少在低海拔群体WL1上发生了作用。 通过叶形态、生理及DNA分子水平的研究，结果表明叶形态特征和碳同位素与海拔紧密相关，与海拔之间呈非线性变化，海拔2,800 m附近是川滇高山栎生长和发育的最适地带。海拔梯度在一定程度上会影响到川滇高山栎群体的遗传变异结构，但在这样一个狭窄的地理分布区域里，这种影响并不足以导致群体间较大的遗传分化。同时生境条件及人为干扰也是影响遗传变异的限制性因子，不容忽视。 Altitudinal gradients impose heterogeneous environmental conditions, such as rugged topography, a complex pattern of vegetation and flowering delay, and they likely furthermore markedly affect the morphological and genetic variation pattern of a species. Understanding altitudinal pattern of morphological and genetic variation at a species is important for the management and conservation of species diversity. Although plant population genetics is a fast growing field of research, there are only few recent investigations, which analyzed the genetic differentiation and changes of intra-population variation along altitudinal gradients. At present, it is still unclear whether there are some common patterns of morphological and genetic variation with altitude. Quercus aquifolioides Rehder & E.H. Wilson, which is an important ecological and economical endemic woody plant species, is widely distributed in the Yunnan and Sichuan provinces, Southwest China. Its large range of habitat across different altitudes implies strong adaptation to different environments, although it is mainly restricted to sunny, south facing slopes. It plays a very important role in preventing soil erosion, soil water loss and regulating climate, as well as in retaining ecological stability. In this paper, we tried to understand the altitudinal pattern of morphological and genetic variation along altitudinal gradients through the experiments of leaf morphological and physiological responses, microsatellite analysis and AFLP markers. In leaf morphological and physiological responses experiment, we measured leaf morphology, nitrogen content and carbon isotope composition (as an indicator of water use efficiency) of Q. aquifolioides along an altitudinal gradient. We found that these leaf morphological and physiological responses to altitudinal gradients were non-linear with increasing altitude. Specific leaf area, stomatal length and index increased with increasing altitude below 2,800 m, but decreased with increasing altitude above 2,800 m. In contrast, leaf nitrogen content per unit area and carbon isotope composition showed opposite change patterns. Specific leaf area seemed to be the most important parameter that determined the carbon isotope composition along the altitudinal gradient. Our results suggest that near 2,800 m in altitude could be the optimum zone for growth and development of Q. aquifolioides, and highlight the importance of the influence of altitude in research on plant physiological ecology. Genetic variation and differentiation were investigated among five natural populations of Q. aquifolioides occurring along an altitudinal gradient that varied from 2,000 to 3,600 m above sea level in the Wolong Natural Reserve of China, by analyzing variation at six microsatellite loci. The results showed that the populations were characterized by relatively high intra-population variation with the average number of alleles equaling 11.33 per locus and the average expected heterozygosity (HE) being 0.779. The amount of genetic variation varied only little among populations, which suggests that the influence of altitude factors on microsatellite variation is limited. However, there is a significantly positive correlation between altitude and the number of low-frequency alleles (R2=0.97, P < 0.01), which indicates that Q. aquifolioides from high altitudes has more unique variation, possibly enabling adaptation to severe conditions. F statistics showed the presence of a slight deficiency of heterozygosity (FIS=0.136) and a low level of differentiation among populations (FST=0.066). The result of the cluster analysis demonstrates that the grouping of populations does not correspond to the altitude of the populations. Based on the available data, it is likely that the selective forces related to altitude are not strong enough to significantly differentiate the populations of Q. aquifolioides in terms of microsatellite variation. To further elucidate genetic variation pattern of Q. aquifolioides populations under sub-alpine environments, genetic variation and differentiation were investigated along altitudinal gradients using AFLP markers. The altitudinal populations with an average altitude interval of 400 m, i.e. WL1, WL2, WL3, WL4 and WL5, correspond to the altitudes 2,000, 2,400, 2,800, 3,200 and 3,600 m, respectively. Our results were as follows: (i) decreasing genetic variation (ranging from 0.253 to 0.210) and increasing genetic differentiation with altitude were obtained from the WL2 to the WL5 population; (ii) the WL1 population showed the lowest genetic variation (HE = 0.181) and the highest genetic differentiation (average FST = 0.0596) with the other four populations; (iii) the positive correlation was obtained using Mantel tests between genetic and altitude distances except for the WL1 population. Our results suggest that altitudinal gradients may have influenced the genetic variation pattern of Q. aquifolioides populations to some extent. In addition, habitat environments (unfavorable wet and hot conditions) and human disturbances (burning, grazing and felling) were possible influencing factors, especially to the low-altitude WL1 population. The present study shows that there were close correlations between morphological features and carbon isotope composition in our data. This indicates that a coordinated plant response modified these parameters simultaneously across different altitudes. Around 2,800 m altitude there seems to be an optimum zone for growth and development of Q. aquifolioides, as indicated by thinner leaves, larger stomata and more negative d13C values. All available evidence indicates altitudinal gradients may have influenced the genetic variation pattern of Q. aquifolioides to some extent. Decreasing genetic variation and increasing genetic differentiation with altitude was obtained except for the WL1 population. And the environment of habitats and human disturbances were also contributing factors, which impact genetic variation pattern, especially to the low-altitude WL1 population.

土壤侵蚀高维指数非线性模型及其参数辨识

模型验证和数据库组建基础上,用WinEPIC模型定量模拟研究了黄土高原半湿润区长武、半干旱区固原和半干旱偏旱区海原20～30年内苜蓿草地水分生产潜力、10m土层土壤有效含水量和土壤湿度剖面分布特征的动态变化.结果表明:长武、固原和海原苜蓿草地水分生产潜力模拟值随降水量变化而呈现波动性降低趋势,其平均值分别为8.81、3.83和2.48t.hm-2;长武、固原和海原苜蓿草地10m土层逐月土壤有效含水量模拟值均呈现明显的波动性降低趋势,模拟初期,4～8年生苜蓿草地土壤干燥化趋势十分强烈,此后,随降水量变化长期在较低水平上波动;随着苜蓿生长年限的延长,苜蓿草地土壤干层逐年加深、加厚,长武、固原和海原土壤干层分布深度达到10m所需时间依次为6、6和4年,此后苜蓿草地降水渗深以下土层长期维持较为稳定的干燥化状态;苜蓿草地水分持续利用的合理年限为半湿润区8～10年,半干旱区6～8年,半干旱偏旱区4～6年.

地球介质非线性与各向异性的地震波响应及其成像效应

The practice of geophysical prospecting shows us the complex interior earth. The studies of the complexity play an important role and practical guide for the subsurface structure. At present, the complexity of the earth mainly means lateral and vertical homogeneity, anisotropy and non-linear quality. And the anisotropy and non-linear media studies become the frontier in seismology and exploration seismology. This paper summarizes the development of complexities and presents the forward and inverse in the non-linear and anisotropic media. Firstly, the paper introduces the theory of seismic wave propagation in the non-linear and anisotropic media, the theoretical basis for simulation and inversion research. Secondly, high quality numerical simulation method with little dispersion has been developed to investigate the influence of complexity including anisotropy and non-linear multi-component seismograms. Because most real data in seismology have a single component, we developed two aspects work on anisotropic multi-component imaging. One is prestack reflection migration. The result show that distorted images are obtained if data from anisotropic media are migrated using isotropic extrapolation. Moreover, image quality will be improved greatly after considering anisotropy in subsurface layers. The other one is the we take advantage of multi-component data to inversion of the anisotropic parameters jointly seimic reflection travel time and polarization information. Based on these research works, we get the following results: 1.Combing numerical simulation, systematical studies indicate that anisotropy and non-linear seismograms characters are significant to detect cracked belts in the earth and to understand deformation field and mechanism. 2.Based on anisotropic media models, we developed an efficient prestack migration method for subsurface structure and different observation methods seismic data, which improving the imaging quality with VSP, seismograms and real data. 3.Jointly seismic inversion combining seismic anisotropic reflection traveltimes and polarizations data show that the complete wrong inversion and the following explanation will be resulted by ignoring anisotropy.

层析静校正技术研究和应用

Static correction is one of the indispensable steps in the conventional onshore seismic data processing, particularly in the western part of China; it is theoretically and practically significant to resolve the issue of static correction. Conventional refraction static correction is put forward under the assumption that layered medium is horizontal and evenly distributed. The complicated nature of the near surface from western part of China is far from the assumption. Therefore, the essential way to resolve the static correction problem from the complex area is to develop a new theory. In this paper, a high-precision non-linear first arrival tomography is applied to solve the problem, it moved beyond the conventional refraction algorithm based on the layered medium and can be used to modeling the complex near surface. Some of the new and creative work done is as follows: One. In the process of first arrival tomographic image modeling, a fast high-order step algorithm is used to calculate the travel time for first arrival and ray path and various factors concerning the fast step ray tracing algorithm is analyzed. Then the second-order and third-order differential format is applied to the step algorithm which greatly increased the calculation precision of the ray tracing and there is no constraint to the velocity distribution from the complex areas. This method has very strong adaptability and it can meet the needs of great velocity variation from the complicated areas. Based on the numerical calculation, a fast high-order step is a fast, non-conditional and stable high-precision tomographic modeling algorithm. Two, in the tomographic inversion, due to the uneven fold coverage and insufficient information, the inversion result is unstable and less reliable. In the paper, wavelet transform is applied to the tomographic inversion which has achieved a good result. Based on the result of the inversion from the real data, wavelet tomographic inversion has increased the reliability and stability of the inversion. Three. Apply the constrained high-precision wavelet tomographic image to the static correction processing from the complex area. During tomographic imaging, by using uphole survey, refraction shooting or other weathering layer method, weathering layer can be identified before the image. Because the group interval for the shot first arrival is relatively big, there is a lack of precision for the near surface inversion. In this paper, an inversion method of the layer constraint and well constraint is put forward, which can be used to compensate the shallow velocity of the inversion for the shot first arrival and increase the precision of the tomographic inversion. Key words: Tomography ，Fast marching method，Wavelet transform, Static corrections, First break

青藏高原及邻区面波频散反演--岩石圈-软流层体系结构研究

The Qinghai-Tibet Plateau lies in the place of the continent-continent collision between Indian and Eurasian plates. Because of their interaction the shallow and deep structures are very complicated. The force system forming the tectonic patterns and driving tectonic movements is effected together by the deep part of the lithosphere and the asthenosphere. It is important to study the 3-D velocity structures, the spheres and layers structures, material properties and states of the lithosphere and the asthenosphere for getting knowledge of their formation and evolution, dynamic process, layers coupling and exchange of material and energy. Based on the Rayleigh wave dispersion theory, we study the 3-D velocity structures, the depths of interfaces and thicknesses of different layers, including the crust, the lithosphere and the asthenosphere, the lithosphere-asthenosphere system in the Qinghai-Tibet Plateau and its adjacent areas. The following tasks include: (1)The digital seismic records of 221 seismic events have been collected, whose magnitudes are larger than 5.0 over the Qinghai-Tibet Plateau and its adjacent areas. These records come from 31 digital seismic stations of GSN , CDSN、NCDSN and part of Indian stations. After making instrument response calibration and filtering, group velocities of fundamental mode of Rayleigh waves are measured using the frequency-time analysis (FTAN) to get the observed dispersions. Furthermore, we strike cluster average for those similar ray paths. Finally, 819 dispersion curves (8-150s) are ready for dispersion inversion. (2)From these dispersion curves, pure dispersion data in 2°×2° cells of the areas (18°N-42°N, 70°E-106°E) are calculated by using function expansion method, proposed by Yanovskaya. The average initial model has been constructed by taking account of global AK135 model along with geodetic, geological, geophysical, receiving function and wide-angle reflection data. Then, initial S-wave velocity structures of the crust and upper mantle in the research areas have been obtained by using linear inversion (SVD) method. (3)Taking the results of the linear inversion as the initial model, we simultaneously invert the S wave velocities and thicknesses by using non-linear inversion (improved Simulated Annealing algorithm). Moreover, during the temperature dropping the variable-scale models are used. Comparing with the linear results, the spheres and layers by the non-linear inversion can be recognized better from the velocity value and offset. (4)The Moho discontinuity and top interface of the asthenosphere are recognized from the velocity value and offset of the layers. The thicknesses of the crust, lithosphere and asthenosphere are gained. These thicknesses are helpful to studying the structural differentia between the Qinghai-Tibet Plateau and its adjacent areas and among geologic units of the plateau. The results of the inversion will provide deep geophysical evidences for studying deep dynamical mechanism and exploring metal mineral resource and oil and gas resources. The following conclusions are reached by the distributions of the S wave velocities and thicknesses of the crust, lithosphere and asthenosphere, combining with previous researches. (1)The crust is very thick in the Qinghai-Tibet Plateau, varying from 60 km to 80 km. The lithospheric thickness in the Qinghai-Tibet Plateau is thinner (130-160 km) than its adjacent areas. Its asthenosphere is relatively thicker, varies from 150 km to 230 km, and the thickest area lies in the western Qiangtang. India located in south of Main Boundary thrust has a thinner crust (32-38 km), a thicker lithosphere of about 190 km and a rather thin asthenosphere of only 60 km. Sichuan and Tarim basins have the crust thickness less than 50km. Their lithospheres are thicker than the Qinghai-Tibet Plateau, and their asthenospheres are thinner. (2)The S-wave velocity variation pattern in the lithosphere-asthenosphere system has band-belted distribution along east-westward. These variations correlate with geology structures sketched by sutures and major faults. These sutures include Main Boundary thrust (MBT), Yarlung-Zangbo River suture (YZS), Bangong Lake-Nujiang suture (BNS), Jinshajiang suture (JSJS), Kunlun edge suture (KL). In the velocity maps of the upper and middle crust, these sutures can be sketched. In velocity maps of 250-300 km depth, MBT, BNS and JSJS can be sketched. In maps of the crustal thickness, the lithospheric thickness and the asthenospheric thickness, these sutures can be still sketched. In particular, MBT can be obviously resolved in these velocity maps and thickness maps. (3)Since the collision between India and Eurasian plate, the “loss” of surface material arising from crustal shortening is caused not only by crustal thickening but also by lateral extrusion material. The source of lateral extrusion lies in the Qiangtang block. These materials extrude along the JSJS and BNS with both rotation and dispersion in Daguaiwan. Finally, it extends toward southeast direction. (4)There is the crust-mantle transition zone of no distinct velocity jump in the lithosphere beneath the Qiangtang Terrane. It has thinner lithosphere and developed thicker asthenosphere. It implies that the crust-mantle transition zone of partial melting is connected with the developed asthenosphere. The underplating of asthenosphere may thin the lithosphere. This buoyancy might be the main mechanism and deep dynamics of the uplift of the Qinghai-Tibet hinterland. At the same time, the transport of hot material with low velocity intrudes into the upper mantle and the lower crust along cracks and faults forming the crust-mantle transition zone.

地震数据波动方程反演方法研究

Geophysical inversion is a theory that transforms the observation data into corresponding geophysical models. The goal of seismic inversion is not only wave velocity models, but also the fine structures and dynamic process of interior of the earth, expanding to more parameters such as density, aeolotropism, viscosity and so on. As is known to all, Inversion theory is divided to linear and non-linear inversion theories. In rencent 40 years linear inversion theory has formed into a complete and systematic theory and found extensive applications in practice. While there are still many urgent problems to be solved in non-linear inversion theory and practice. Based on wave equation, this dissertation has been mainly involved in the theoretical research of several non-linear inversion methods: waveform inversion, traveltime inversion and the joint inversion about two methods. The objective of gradient waveform inversion is to find a geologic model, thus synthetic seismograms generated by this geologic model are best fitted to observed seismograms. Contrasting with other inverse methods, waveform inversion uses all characteristics of waveform and has high resolution capacity. But waveform inversion is an interface by interface method. An artificial parameter limit should be provided in each inversion iteration. In addition, waveform information will tend to get stuck in local minima if the starting model is too far from the actual model. Based on velocity scanning in traditional seismic data processing, a layer-by-layer waveform inversion method is developed in this dissertation to deal with weaknesses of waveform inversion. Wave equation is used to calculate the traveltime and derivative (perturbation of traveltime with respect to velocity) in wave-equation traveltime inversion (WT). Unlike traditional ray-based travetime inversion, WT has many advantages. No ray tracing or traveltime picking and no high frequency assumption is necessary and good result can be got while starting model is far from real model. But, comparing with waveform inversion, WT has low resolution. Waveform inversion and WT have complementary advantages and similar algorithm, which proves that the joint inversion is a better inversion method. And another key point which this dissertation emphasizes is how to give fullest play to their complementary advantages on the premise of no increase of storage spaces and amount of calculation. Numerical tests are implemented to prove the feasibility of inversion methods mentioned above in this dissertation. Especially for gradient waveform inversion, field data are inversed. This field data are acquired by our group in Wali park and Shunyi district. Real data processing shows there are many problems for waveform inversion to deal with real data. The matching of synthetic seismograms with observed seismograms and noise cancellation are two primary problems. In conclusion, on the foundation of the former experiences, this dissertation has implemented waveform inversions on the basis of acoustic wave equation and elastic wave equation, traveltime inversion on the basis of acoustic wave equation and traditional combined waveform traveltime inversion. Besides the traditional analysis of inversion theory, there are two innovations: layer by layer inversion of seimic reflection data inversion and rapid method for acoustic wave-equation joint inversion.

基于随机介质理论的多尺度储层非均质性定量描述

Application of long-term exploration for oil and gas shows that the reservoir technology of prediction is one of the most valuable methods. Quantitative analysis of reservoir complexity is also a key technology of reservoir prediction. The current reservoir technologies of prediction are based on the linear assumption of various physical relationships. Therefore, these technologies cannot handle complex reservoirs with thin sands, high heterogeneities in lithological composition and strong varieties in petrophysical properties. Based on the above-mentioned complex reservoir, this paper conducts a series of researches. Both the comprehending and the quantitative analysis of reservoir heterogeneities have been implemented using statistical and non-linear theories of geophysics. At the beginning, the research of random media theories about reservoir heterogeneities was researched in this thesis. One-dimensional (1-D) and two-dimensional (2-D) random medium models were constructed. The autocorrelation lengths of random medium described the mean scale of heterogeneous anomaly in horizontal and deep directions, respectively. The characteristic of random medium models were analyzed. We also studied the corresponding relationship between the reservoir heterogeneities and autocorrelation lengths. Because heterogeneity of reservoir has fractal nature, we described heterogeneity of reservoir by fractal theory based on analyzing of the one-dimensional (1-D) and two-dimensional (2-D) random medium models. We simulated two-dimensional (2-D) random fluctuation medium in different parameters. From the simulated results, we can know that the main features of the two-dimensional (2-D) random medium mode. With autocorrelation lengths becoming larger, scales of heterogeneous geologic bodies in models became bigger. In addition, with the autocorrelation lengths becoming very larger, the layer characteristic of the models is very obvious. It would be difficult to identify sandstone such as gritstone, clay, dense sandstone and gas sandstone and so on in the reservoir with traditional impedance inversion. According to the obvious difference between different lithologic and petrophysical impedance, we studied multi-scale reservoir heterogeneities and developed new technologies. The distribution features of reservoir lithological and petrophysical heterogeneities along vertical and transverse directions were described quantitatively using multi-scale power spectrum and heterogeneity spectrum methods in this paper. Power spectrum (P spectrum) describes the manner of the vertical distribution of reservoir lithologic and petrophysical parameters and the large-scale and small-scale heterogeneities along vertical direction. Heterogeneity spectrum (H spectrum) describes the structure of the reservoir lithologic and petrophysical parameters mainly, that is to say, proportional composition of each lithological and petrophysical heterogeneities are calculated in this formation. The method is more reasonable to describe the degree of transverse multi-scale heterogeneities in reservoir lithological and petrophysical parameters. Using information of sonic logs in Sulige oil field, two spectral methods have been applied to the oil field, and good analytic results have been obtained. In order to contrast the former researches, the last part is the multi-scale character analysis of reservoir based on the transmission character of wave using the wavelet transform. We discussed the method applied to demarcate sequence stratigraphy and also analyzed the reservoir interlayer heterogeneity.