大陆地壳岩石结构重建：方法与实例

The composition of the continental crust has long been a subject of interest to earth scientists as it can provide key information about the crustal growth and evolution of the continents. In this paper we make a comparative study on the lithological discrimination schemes featuring with the use of different seismic attributes, such as P-wave velocity, P- to S-wave velocity ratio, acoustic or elastic impedances, Lame impedances and high-sensitive identification factors. The results demonstrate that Lame impedances have more powerful constrains than other seismic attributes. In order to fully take the advantage of make the best of the different seismic response of crustal rock, we firstly use seismic attribute that have weak distinguish power to construct loose constrained lithological model, then use seismic attributes that have stronger distinguish power to tighten the constrains of lithological discrimination. We propose a joint scheme (chain constrain technique) by combing all available constrains to reduce the non-uniqueness in mapping rock distribution. We adopt chain constrain technique to construct lithological model beneath Tunxi-Wenzhou transect, Southeastern China, Manzhouli-Suifenhe transect, Northeastern China, and geophysical profile in Bohai Bay Basin, North China. The results can be suumarized as the follows: (1) We compare the sensitivity of different seismic factor constraints on rock types, and conclude that Lame impedances have tighter constrains than seismic velocity, Vp/Vs, density. (2) We propose chain constrains to construct lithological model from integrated geophysical data, and reduce the non-uniqueness in mapping rock distribution. (3) We reconstruct crustal lithological model beneath Tunxi-Wenzhou transect, Southeastern China. The results suggested that Jiangshan-Shaoxing fault is a crust-scale, and it is the boundary between Cathaysia and Yanthze blocks. (4) We construct crustal lithological model beneath Manzhouli-Suifenhe transect, Northeastern China. (5) We map the petrologic distribution along a geophysical profile in Bohai Bay Basin, North China, and construct a three-layered petrology model from the depth 2 km to about 10 km, consisted of basalt (the first layer), pelitic siltstone (the second layer), and silty mudstone and fine sandstone (the third layer).

AVO分析与地震资料叠前储层预测技术

In this paper we base on the anisotropic theory and Zoeppritz function of the transmission theory and the law of amplitude versus offset simplify seismic reflection coefficient of different media, analyze the characteristic of the gas or oil saturated stratum or the VTI and HTI models. Discuss the P wave reflection relationship and the meanings of the different parameters. We use measured parameters of a reservoir to simulate the characteristic of the reservoir, study the different effects of stratum saturated with gas or oil and analyze the characteristic of the seismic response of different models which change with different incident angles and different azimuths. Using the field data of logs ,analyze the rock property parameters, build the relationship of logs and parameters by Gassmann theory or empirical function. Calculate the density and the shear modulus and bulk modulus, reconstruct the log curves, calculate shear wave logs and correlate the logs affected by mud and other environmental factors. Finally perform the relationship of the seismic data log of saturated stratum and enhance the ability and reliability in reservoir prediction. Our aim is by the prestack seismic processing to get high solution and amplitude preserved seismic data. Because in incident angle gathers or azimuthal gathers, the low signal to noise ratio and low different covers affect the result of the prestack reservoir prediction. We apply prestack noise erase, cell regularization process and relatively amplitude preservation in the high solution seismic process routine to preserve the characteristic of stratum response, and erase the effects of the noise. In this paper we finished prestack invertion in the BYT survey and fractured reservoir depiction in MB survey. By the invertion and multiple attributes crossplot. we can get the stratum profiles and oil indicator profiles which can predict the distribution of the reservoir and oil. In the MB survey, we get orientation and density of fractured reservoir by the azimuthal seismic amplitude and depict the potential oil and gas reservoir. Prestak invertion works better in distinguishing oil and reservoir.

离散裂缝诱导各向异性粘弹性地震波模拟与特征分析

The real media always attenuate and distort seismic waves as they propagate in the earth. This behavior can be modeled with a viscoelastic and anisotropic wave equation. The real media can be described as fractured media. In this thesis, we present a high-order staggered grid finite-difference scheme for 2-D viscoelastic wave propagation in a medium containing a large number of small finite length fractures. We use the effective medium approach to compute the anisotropic parameters in each grid cell. By comparing our synthetic seismogram by staggered-grid finite-difference with that by complex-ray parameter ray tracing method, we conclude that the high-order staggered-grid finite-difference technique can effectively used to simulate seismic propagation in viscoelastic-anisotropic media. Synthetic seismograms demonstrate that strong attenuation and significant frequency dispersion due to viscosity are important factors of reducing amplitude and delaying arrival time varying with incidence angle or offset. On the other hand, the amount of scattered energy not only provides an indicator of orientation of fracture sets, but can also provide information about the fracture spacing. Analysis of synthetic seismograms from dry- and fluid-filled fractures indicates that dry-filled fractures show more significant scattering on seismic wavefields than fluid-filled ones, and offset-variations in P-wave amplitude are observable. We also analyze seismic response of an anticlinal trap model that includes a gas-filled fractured reservoir with high attenuation, which attenuates and distorts the so-called bright spot.

基于工程地质力学的岩土工程加固理论、技术与实践

In the engineering reinforcement of-rock and soil mass, engineers must consider how to obtain better reinforcing effect at the cost of less reinforcing expense, which, in fact, is the aim of reinforcement design. In order to accomplish the purpose, they require not only researching the material used to reinforce and its structure, but also taking into account of several important geological factors, such as the structure and property of rock and soil mass. How to improve the reinforcing effect according to engineering geomechanical principle at the respect of the reinforcement of engineering soil and rock mass is studied and discussed in this paper. The author studies the theory, technology and practice of geotechnical reinforcement based on engineering geomechanics, taking example for the soil treatment of Zhengzhou Airport, the effect analysis of reinforcement to the slope on the left bank of Wuqiangxi Hydropower Station and the reinforcing design of the No. 102 Landslide and unique sand-slide slope on the Sichuan-Tibet Highway. The paper is comprised of two parts for the convenience of discussion. In the first part, from the first chapter to the fifth chapter, trying to perform the relevant research and application at the viewpoint of soil mass engineering geomechanics, the author mainly discusses the study of reinforcing soft ground soil through dynamical consolidation and its application. Then, in the second part, from the sixth chapter to the eleventh chapter, the study of new technologies in the rock slope reinforcement and their application are discussed. The author finds that not only better reinforcing effect can be gained in the research where the principle and method of rock mass engineering geomechanics is adopted, but also new reinforcing technologies can be put forward. Zhengzhou Airport is an important one in central plains. It lies on Yellow River alluvial deposit and the structure of stratum is complex and heterogeneous. The area of airport is very large, which can result in differential settlement easily, damage of airport and aircraft accident, whereas, there are no similar experiences to dispose the foundation, so the foundation treatment become a principal problem. During the process of treatment, the method of dynamic compaction was adopted after compared with other methods using the theory of synthetic integration. Dynamic compaction is an important method to consolidate foundation, which was successfully used in the foundation of Zhengzhou Airport. For fill foundation, controlling the thickness of fill so as to make the foundation treatment can reach the design demand and optimum thickness of the fill is a difficult problem. Considering this problem, the author proposed a calculation method to evaluate the thickness of fill. The method can consider not only the self-settlement of fill but also the settlement of the ground surface under applied load so as to ensure the settlement occurred during the using period can satisfy the design demand. It is proved that the method is correct after using it to choose reasonable energy of dynamic compaction to treat foundation. At the same time, in order to examine the effect of dynamic compaction, many monitor methods were adopted in the test such as static loading test, modulus of resilience test, deep pore pressure -test, static cone penetration test and the variation of the pore volume measurement. Through the tests, the author summarized the discipline of the accumulation and dissipation of pore pressure in Yellow River alluvial deposit under the action of dynamic compaction, gave a correct division of the property change of silt and clay under dynamic compaction, determined the bearing capacity of foundation after treatment and weighted the reinforcing effect of dynamic consolidation from the variation of the soil particle in microcosmic and the parameter of soil mass' density. It can be considered that the compactness of soil is in proportion to the energy of dynamic compaction. This conclusion provided a reference to the research of the "Problem of Soil Structure-the Central Problem of Soil Mechanics in 21 Century ". It is also important to strengthen rock mass for water conservancy and electric power engineering. Slip-resistance pile and anchoring adit full of reinforced concrete are usually adopted in engineering experience to strengthen rock mass and very important for engineering. But there also some deficiency such as the weakest section can't be highlighted, the monitor is inconvenient and the diameter of pile and adit is very large etc. The author and his supervisor professor Yangzhifa invented prestressed slip-resistance pile and prestressed anchoring adit full of reinforced concrete, utilizing the advantage that the prestressed structure has better anti-tensile characteristic (this invention is to be published). These inventions overcome the disadvantages of general slip-resistance pile and anchoring adit full of reinforced concrete and have the functions of engineering prospecting, strengthening, drainage and monitor simultaneous, so they have better strengthened effect and be more convenient for monitor and more economical than traditional methods. Drainage is an important factor in treatments of rock mass and slop. In view of the traditional drainage method that drainage pore often be clogged so as to resulted in incident, professor Yangzhifa invented the method and setting of guide penetration by fiber bundle. It would take good effect to use it in prestressed slip-resistance pile and anchoring adit full of reinforced concrete. In this paper, the author took example for anchoring adit full of reinforced concrete used to strengthen Wuqiangxi left bank to simulate the strengthened effect after consolidated by prestressed slip-resistance pile, took example for 102 landslide occurred along Sichuan-Tibet highway to simulate the application of slip-resistance pile and the new technology of drainage. At the same time the author proposed the treatment method of flowing sand in Sichuan-Tibet highway, which will benefit the study on strengthening similar engineering. There are five novelties in the paper with the author's theoretical study and engineering practice: 1. Summarizing the role of pore water pressure accumulation and dissipation of the Yellow River alluvial and diluvial soil under the action of dynamical consolidation, which has instructive significance in the engineering construction under the analogical engineering geological conditions in the future. It has not been researched by the predecessors. 2. Putting forward the concept of density D in microcosmic based on the microcosmical structure study of the soil sample. Adopting D to weight the reinforcing effect of dynamic consolidation is considered to be appropriate by the means of comparing the D values of Zhengzhou Airport's ground soil before with after dynamically consolidating reinforcement, so a more convenient balancing method can be provided for engineering practice. 3. According to the deep research into the soil mass engineering geology, engineering rock and soil science, soil mechanics, as well as considerable field experiments, improving the consolidating method in airport construction, from the conventional method, which is dynamically compactmg original ground surface firstly, then filling soil and dynamically layer-consolidating or layer-compacting at last to the upgraded method, which is performing dynamical consolidation after filling soil to place totally at the extent of the certain earth-filling depth. The result of the dynamical consolidation not only complies with the specifications, but also reduces the soil treatment investment by 10 million RMB. 4. Proposing the method for calculating the height of the filled soil by the means of estimating the potential displacement produced in the original ground surface and the filled earth soil under the possible load, selecting the appropriate dynamically-compacting power and determining the virtual height of the filled earth soil. The method is proved to be effective and scientific. 5. According to the thought of Engineering Geomechanics Metal-Synthetic Methodology (EGMS), patenting two inventions (to the stage of roclamation, with Professor Yang Zhi-fa, the cooperative tutor, and etc.) in which multi-functions, engineering geological investigation, reinforcement, drainage and strength remedy, are integrated all over in one body at the viewpoint of the breakage mechanism of the rock slope.

接收函数方法与中国及邻区上地幔间断面的研究

The receiver function method applied in researching the discontinuities in upper mantle was systematically studied in this paper. Using the theoretical receiver functions, the characteristics of P410S and P660S phases were analyzed, and the influencing factors for detection of these phases were discussed. The stability of receiver function was studied, and a new computational method of receiver function, RFSSMS (Receiver Function of Stack and Smooth of Multi seismic-records at a Single station), was put forward. We built initial reference velocity model for the media beneath each of 18 seismic stations respectively; then estimated the buried depths of 410-km and 660-km discontinuities(simply marked as '410' and '660') under the stations by using the arrive time differences of P410S and P660S with P. We developed a new receiver function inversion method -PGARFI (Peeling-Genetic Algorithm of Receiver Function Inversion), to obtain the whole crust and upper mantle velocity structure and the depths of discontinuities beneath a station. The major works and results could be summarized as follows: (1) By analysis of the theoretical receiver functions with different velocity models and different ray parameters, we obtain the knowledge: The amplitudes of P410S and P660S phases are decreasing with the increasing of epicentral distance A , and the arrival time differences of these phases with P are shorter as A is longer. The multiple refracted and/or reflected waves yielded on Moho and the discontinuities in the crust interfere the identification of P410S. If existing LVZ under the lithosphere, some multiple waves caused by LVZ will interfere the identification of P410S. The multiple waves produced by discontinuity lied near 120km depth will mix with P410s phase in some range of epicentral distance; and the multiple waves concerned with the discontinuity lied near 210km depth will interfere the identification of P660S. The epicentral distance for P4i0s identification is limited, the upper limit is 80° . The identification of P660S is not restricted by the epicenter distance obviously. The identification of P410S and P6gos in the theoretical receiver functions is interfered weakly from the seismic wave attenuation caused by the media absorption if the Q value in a reasonable range. (2) The stability of receiver function was studied by using synthetic seismograms with different kind of noise. The results show that on the condition of high signal-noise-ratio of seismic records, the high frequency background noise and the low frequency microseism noise do not influence the calculating result of receiver function. But the media "scattering noise" influence the stability of receiver function. When the scattering effect reach some level, the identification of P4iOs and P66os is difficult in single receiver function which is yielded from only one seismic record. We provided a new method to calculate receiver function, that is, with a group of earthquake records, stacking the R and Z components respectively in the frequency domain, and weighted smooth the stacked Z component, then compute the complex spectrum ratio of R to Z. This method can improve the stability of receiver function and protrude the P4i0s and P66os in the receiver function curves. (3) 263 receiver functions were provided from 1364 three component broadband seismograms recorded at 18 stations in China and adjacent areas for the tele-earthquakes. The observed arrival time differences of P410S and P660S with P were obtained in these receiver functions. The initial velocity model for every station was built according to the prior research results. The buried depths of '410' and '660' under a station were acquired by the way of adjusting the depths of these two discontinuities in the initial velocity model until the theoretical arrival time differences of P410S and P660S with P well conformed to the observed. The results show an obvious lateral heterogeneity of buried depths of ' 410' and (660' . The depth of '410' is shallower beneath BJI, XAN, LZH and ENH, but deeper under QIZ and CHTO, and the average is 403km . The average depth of '660' is 663km, deeper under MDJ and MAJO, but shallower under QIZ and HYB. (4) For inversing the whole crust and upper mantle velocity structure, a new inversion method -PGARFI (Peeling-Genetic Algorithm of Receiver Function Inversion) has- been developed here. The media beneath a station is divided into segments, then the velocity structure is inversed from receiver function from surface to deep successively. Using PGARFI, the multi reflection / refraction phases of shallower discontinuities are isolated from the first order refraction transform phase of deep discontinuity. The genetic algorithm with floating-point coding was used hi the inversion of every segment, and arithmetical crossover and non-uniform mutation technologies were employed in the genetic optimization. 10 independent inversions are completed for every segment, and 50 most excellent velocity models are selected according to the priority of fitness from all models produced in the inversion process. The final velocity structure of every segment is obtained from the weighted average of these 50 models. Before inversion, a wide range of velocity variation with depth and depth range of the main discontinuities are given according to priori knowledge. PGARFI was verified with numerical test and applied in the inversion of the velocity structure beneath HIA station down to 700km depth.

双聚焦偏移速度估计与可分表示叠前弹性反演

In exploration seismology, the geologic target of oil and gas reservoir in complex medium request the high accuracy image of the structure and lithology of the medium. So the study of the prestack image and the elastic inversion of seismic wave in the complex medium come to the leading edge. The seismic response measured at the surface carries two fundamental pieces of information: the propagation effects of the medium and the reflections from the different layer boundaries in the medium. The propagation represent the low-wavenumber component of the medium, it is so-called the trend or macro layering, whereas the reflections represent the high-wavenumber component of the medium, it is called the detailed or fine layering. The result of migration velocity analysis is the resolution of the low-wavenumber component of the medium, but the prestack elastic inversion provided the resolution of the high-wavvenumber component the medium. In the dissertation, the two aspects about the migration velocity estimation and the elastic inversion have been studied.Firstly, any migration velocity analysis methods must include two basic elements: the criterion that tell us how to know whether the model parameters are correct and the updating that tell us how to update the model parameters when they are incorrect, which are effected on the properties and efficiency of the velocity estimation method. In the dissertation, a migration velocity analysis method based on the CFP technology has been presented in which the strategy of the top-down layer stripping approach are adapted to avoid the difficult of the selecting reduce .The proposed method has a advantage that the travel time errors obtained from the DTS panel are defined directly in time which is the difference with the method based on common image gather in which the residual curvature measured in depth should be converted to travel time errors.In the proposed migration velocity analysis method, the four aspects have been improved as follow:? The new parameterization of velocity model is provided in which the boundaries of layers are interpolated with the cubic spline of the control location and the velocity with a layer may change along with lateral position but the value is calculated as a segmented linear function of the velocity of the lateral control points. The proposed parameterization is suitable to updating procedure.? The analytical formulas to represent the travel time errors and the model parameters updates in the t-p domain are derived under local lateral homogeneous. The velocity estimations are iteratively computed as parametric inversion. The zero differential time shift in the DTS panel for each layer show the convergence of the velocity estimation.? The method of building initial model using the priori information is provided to improve the efficiency of velocity analysis. In the proposed method, Picking interesting events in the stacked section to define the boundaries of the layers and the results of conventional velocity analysis are used to define the velocity value of the layers? An interactive integrate software environment with the migration velocity analysis and prestack migration is built.The proposed method is firstly used to the synthetic data. The results of velocity estimation show both properties and efficiency of the velocity estimation are very good.The proposed method is also used to the field data which is the marine data set. In this example, the prestack and poststack depth migration of the data are completed using the different velocity models built with different method. The comparison between them shows that the model from the proposed method is better and improves obviously the quality of migration.In terms of the theoretical method of expressing a multi-variable function by products of single-variable functions which is suggested by Song Jian (2001), the separable expression of one-way wave operator has been studied. A optimization approximation with separable expression of the one-way wave operator is presented which easily deal with the lateral change of velocity in space and wave number domain respectively and has good approach accuracy. A new prestack depth migration algorithm based on the optimization approximation separable expression is developed and used to testing the results of velocity estimation.Secondly, according to the theory of the seismic wave reflection and transmission, the change of the amplitude via the incident angle is related to the elasticity of medium in the subsurface two-side. In the conventional inversion with poststack datum, only the information of the reflection operator at the zero incident angles can be used. If the more robust resolutions are requested, the amplitudes of all incident angles should be used.A natural separable expression of the reflection/transmission operator is represented, which is the sum of the products of two group functions. One group function vary with phase space whereas other group function is related to elastic parameters of the medium and geological structure.By employing the natural separable expression of the reflection/transmission operator, the method of seismic wave modeling with the one-way wave equation is developed to model the primary reflected waves, it is adapt to a certain extent heterogeneous media and confirms the accuracy of AVA of the reflections when the incident angle is less than 45'. The computational efficiency of the scheme is greatly high.The natural separable expression of the reflection/transmission operator is also used to construct prestack elastic inversion algorithm. Being different from the AVO analysis and inversion in which the angle gathers formed during the prstack migration are used, the proposed algorithm construct a linear equations during the prestack migration by the separable expression of the reflection/transmission operator. The unknowns of the linear equations are related to the elasticity of the medium, so the resolutions of them provided the elastic information of the medium.The proposed method of inversion is the same as AVO inversion in , the difference between them is only the method processing the amplitude via the incident angle and computational domain.

Transonic Aeroelastic Numerical Simulation in Aeronautical Engineering

A lower-upper symmetric Gauss-Seidel (LU-SGS) subiteration scheme is constructed for time-marching of the fluid equations. The Harten-Lax-van Leer-Einfeldt-Wada (HLLEW) scheme is used for the spatial discretization. The same subiteration formulation is applied directly to the structural equations of motion in generalized coordinates. Through subiteration between the fluid and structural equations, a fully implicit aeroelastic solver is obtained for the numerical simulation of fluid/structure interaction. To improve the ability for application to complex configurations, a multiblock grid is used for the flow field calculation and transfinite interpolation (TFI) is employed for the adaptive moving grid deformation. The infinite plate spline (IPS) and the principal of virtual work are utilized for the data transformation between the fluid and structure. The developed code was first validated through the comparison of experimental and computational results for the AGARD 445.6 standard aeroelastic wing. Then, the flutter character of a tail wing with control surface was analyzed. Finally, flutter boundaries of a complex aircraft configuration were predicted.

Nonlinear analysis of pile foundation jacket platform

In this paper, the nonlinear collapse of the BOHAI-8 pile foundation jacket platform has been analyzed. The ultimate load and collapse process of two computational models of the structure are given. One model is of fixed support whose length is eight times the pile leg diameter and the other considers the nonlinearity of the soil-pile interaction.

A new dynamic model for study of dislocation pattern formation

Based on the principle given in nonlinear diffusion-reaction dynamics, a new dynamic model for dislocation patterning is proposed by introducing a relaxation time to the relation between dislocation density and dislocation flux. The so-called chemical potential like quantities, which appear in the model can be derived from variation principle for free energy functional of dislocated media, where the free energy density function is expressed in terms of not only the dislocation density itself but also their spatial gradients. The Linear stability analysis on the governing equations of a simple dislocation density shows that there exists an intrinsic wave number leading to bifurcation of space structure of dislocation density. At the same time, the numerical results also demonstrate the coexistence and transition between different dislocation patterns.

Transonic Aeroelastic Numerical Simulation in Aeronautical Engineering

A lower-upper symmetric Gauss-Seidel (LU-SGS) subiteration scheme is constructed for time-marching of the fluid equations. The Harten-Lax-van Leer-Einfeldt-Wada (HLLEW) scheme is used for the spatial discretization. The same subiteration formulation is applied directly to the structural equations of motion in generalized coordinates. Through subiteration between the fluid and structural equations, a fully implicit aeroelastic solver is obtained for the numerical simulation of fluid/structure interaction. To improve the ability for application to complex configurations, a multiblock grid is used for the flow field calculation and transfinite interpolation (TFI) is employed for the adaptive moving grid deformation. The infinite plate spline (IPS) and the principal of virtual work are utilized for the data transformation between the fluid and structure. The developed code was first validated through the comparison of experimental and computational results for the AGARD 445.6 standard aeroelastic wing. Then, the flutter character of a tail wing with control surface was analyzed. Finally, flutter boundaries of a complex aircraft configuration were predicted.

Flow-induced vibrations of long circular cylinders modeled by coupled nonlinear oscillators

The dynamics of long slender cylinders undergoing vortex-induced vibrations (VIV) is studied in this work. Long slender cylinders such as risers or tension legs are widely used in the field of ocean engineering. When the sea current flows past a cylinder, it will be excited due to vortex shedding. A three-dimensional time domain model is formulated to describe the response of the cylinder, in which the in-line (IL) and cross-flow (CF) deflections are coupled. The wake dynamics, including in-line and cross-flow vibrations, is represented using a pair of non-linear oscillators distributed along the cylinder. The wake oscillators are coupled to the dynamics of the long cylinder with the acceleration coupling term. A non-linear fluid force model is accounted for to reflect the relative motion of cylinder to current. The model is validated against the published data from a tank experiment with the free span riser. The comparisons show that some aspects due to VIV of long flexible cylinders can be reproduced by the proposed model, such as vibrating frequency, dominant mode number, occurrence and transition of the standing or traveling waves. In the case study, the simulations show that the IL curvature is not smaller than CF curvature, which indicates that both IL and CF vibrations are important for the structural fatigue damage.

Coherence measurement of a Ti:sapphire femtosecond pulsed laser

The interference patterns produced by Gaussian-shaped broad-bandwidth femtosecond pulsed laser sources are derived. The interference pattern contains both spatial and temporal properties of laser beam. Interference intensity dependent on the bandwidth of femtosecond laser are given. We demonstrate experimentally both the spatial and the temporal coherence properties of a Ti:sapphire femtosecond pulse laser, as well as its power spectrum by using a pinhole pair.

Influence of temperature and salinity fluctuations on propagation behaviour of partially coherent beams in oceanic turbulence

A theoretical study of the behaviour of partially coherent beams propagating through oceanic turbulence has been performed. Based on the previously developed knowledge of beam spreading of a partially coherent beam in the atmosphere and the spatial power spectrum of the refractive index of ocean water, we study the normalized root-mean-square width of a partially coherent beam on propagation through oceanic turbulence and its turbulence distance which may be a measure of turbulence resistance. Our analysis indicates that the behaviour of partially coherent beams on propagation may be described by the rate of dissipation of the mean-squared temperature chi(T) and that of salinity chi(S). In terms of a quantity w that defines the contributions of the temperature and salinity distributions to the distribution of the refractive index, chi(S) could be written as a function of chi(T) and w. Therefore, the behaviour of partially coherent beams on propagation can be characterized only by chi(T) for a given w. The results are shown for curved surfaces, from which one can see that partially coherent beams exhibit robust turbulence resistance when the water volume has a smaller chi(T).

用于钕玻璃啁啾脉冲放大系统光谱整形的多层介质膜反射镜

在千焦拍瓦高功率啁啾脉冲放大系统设计中,为了尽量消除增益窄化和增益饱和效应的影响,同时尽可能提高高功率激光输出脉冲信噪比参数,激光脉冲时空和光谱的整形问题备受关注.提出一种光谱整形新方法,利用特定结构的多层介质膜反射镜,可实现对大能量高功率啁啾脉冲钕玻璃放大系统中啁啾脉冲的光谱整形.研究结果表明,只要合理选择多层介质膜系的结构参数,可有效地控制其反射率分布,且在保证反射相位基本不变的条件下其调制度可超过60%.针对钕玻璃1053nm波长设计而成的光谱整形反射镜,反射带宽可达到196nm,色分辨率约为0.1nm,在几十纳米波长范围所对应的相位偏差小于12mrad,相当于λ/524,能够满足"神光Ⅱ"千焦拍瓦改造的技术指标要求.

冬小麦抗寒性和光合作用对增强UV-B辐射与温度交互作用的响应

人类活动产生的氯氟烃化合物破坏了大气臭氧层，导致了到达地球表面的UV-B辐射大幅度增加。UV-B辐射增强可以影响到植物的生长、形态与发育等各个方面，因此有关增强UV-B辐射对植物的影响，及其与许多环境因子复合作用的研究都已经广泛开展。但是增强UV-B辐射与温度，特别是与低温的相互作用的研究报道很少。在北半球的晚秋至早春这段时期里，一些越冬生长的植物将面临着UV-B辐射增强和低温的双重胁迫，因此，迫切需要进行UV-B辐射和低温生长环境下植物的响应及其机制的研究。 以人工气候生长室中生长的冬小麦(Triticum aestivum)幼苗为试验材料，研究了低剂量（4.2 kJ m-2 d-1 UV-BBE，LUVB）和较高剂量（7.0 kJ m-2 d-1 UV-BBE，HUVB）UV-B辐射处理对20/16℃条件下幼苗抗寒力的交叉适应性及其抗氧化系统的反应;同时还研究了在两种生长温度（25/20℃和10/5℃）条件下，低剂量（4.2 kJ m-2 d-1 UV-BBE，LUVB）和超高剂量（10.3 kJ m-2 d-1 UV-BBE，SHUVB）UV-B辐射处理幼苗的生长速率、光合与荧光参数、叶黄素循环色素、抗氧化系统、以及抗寒性和酚类物质等生理反应，以期阐明不同温度条件下生长的冬小麦对UV-B辐射的生长、光合作用以及抗寒性响应与适应机制。主要结果如下： 1.在LUVB辐射处理下，在20/16℃和25/20℃条件下生长的冬小麦幼苗LT50值都显著降低，HUVB辐射处理对在20/16℃条件下生长的幼苗LT50值也可以显著降低，而SHUVB辐射对25/20℃条件下生长的幼苗LT50值没有显著影响。但是，LUVB和SHUVB辐射处理都导致了10/5℃条件下生长的幼苗LT50值的显著增加。表明适当的UV-B辐射能增强较高温度（20/16℃或25/20℃）条件下冬小麦幼苗的抗寒力，即表现出对冷冻低温的交叉适应性，但低温（10/5℃）生长条件却削弱了UV-B辐射下冬小麦的抗寒能力。 2.在20/16℃条件下接受UV-B辐射预处理的幼苗在-6℃条件下冷冻胁迫6 h再缓慢恢复6 h后，与未进行UV-B辐射处理的对照相比，其叶片过氧化氢酶（CAT）、愈创木酚过氧化物酶（GPX）、谷胱甘肽还原酶（GR）活性，谷胱甘肽氧化还原比例（GSH/GSSG）都显著提高，而由硫代巴比妥酸反应物质（TBARS）代表的膜质过氧化程度显著低于对照。此外，UV-B辐射期间处理幼苗的H2O2含量较对照显著增加，而冷冻恢复以后却明显低于对照。表明UV-B辐射诱导的抗寒力的提高应该与冷冻恢复后植株体内抗氧化系统的上调表达有关，H2O2可能参与了UV-B辐射对低温的交叉适应的信号传导。 3.除25/20℃生长条件下的LUVB处理的小麦幼苗外，UV-B辐射显著降低幼苗的相对生长速率（RGR）、净光合速率（Pn）、光系统II最大量子产量（Fv/Fm）、光系统II实际量子产量（(F΄m−Fs)/F΄m）以及光化学淬灭（qP），但是UV-B辐射并不影响叶片胞间CO2浓度（Ci），而且冬小麦幼苗生长和光合作用的抑制被增加的UV-B辐射剂量和降低的温度加强。UV-B辐射引起的光抑制由非气孔限制所导致，而且主要与PS II光化学效率降低有关。 4.UV-B辐射显著增加了两个温度条件（20/16℃或25/20℃）下生长的冬小麦幼苗叶黄素循环过程中紫黄素（V）的合成，但抑制了V向玉米黄质（Z）的转化，从而造成了对照与LUVB辐射处理幼苗之间的叶片中脱环氧化比例（DEPS）和NPQ无显著性差异，但SHUVB辐射处理幼苗叶片中DEPS和NPQ显著降低。因此，在本试验条件下，增强UV-B辐射处理的冬小麦可能并不通过热耗散形式形成光保护机制，光抑制形成的过剩激发能的耗散可能更多地通过代谢途径来实现。 5.UV-B辐射处理提高了在25/20℃条件下幼苗的超氧化物歧化酶（SOD）、抗坏血酸过氧化物酶（APX）和GR等活性，以及抗坏血酸氧化还原比例（AsA/DHA）和GSH/GSSG;但是在10/5℃下，UV-B辐射除了导致SOD和CAT活性升高之外，对APX活性和AsA/DHA并不产生明显影响，但GPX和GSH/GSSG则显著降低。说明UV-B辐射幼苗的抗氧化系统在较高生长温度下显著地增强，而在低温10/5℃下被严重地削弱或降低，即低温阻止了代谢途径的光保护机制的正常运转。 6.多酚物质在UV-B辐射或低温10/5℃条件下都能显著地累积，且在UV-B辐射和低温复合作用下增加尤其显著，表明多酚物质在两个温度生长条件下特别是低温条件下都参与了对UV-B辐射幼苗的保护。 7.在高温条件下仅仅SHUVB处理的幼苗TBARS含量显著增加，而低温10/5℃条件下两个UV-B辐射处理都非常显著地上升，说明与高温生长条件相比较，低温加重了UV-B辐射引起的氧化胁迫，低温10/5℃条件下幼苗多酚的增加以及抗氧化系统的部分增强都没有能阻止UV-B辐射对幼苗的伤害。