Modelling and simulation of the second-order Doppler error of a laser dual-frequency interferometer

Only the first- order Doppler frequency shift is considered in current laser dual- frequency interferometers; however; the second- order Doppler frequency shift should be considered when the measurement corner cube ( MCC) moves at high velocity or variable velocity because it can cause considerable error. The influence of the second- order Doppler frequency shift on interferometer error is studied in this paper, and a model of the second- order Doppler error is put forward. Moreover, the model has been simulated with both high velocity and variable velocity motion. The simulated results show that the second- order Doppler error is proportional to the velocity of the MCC when it moves with uniform motion and the measured displacement is certain. When the MCC moves with variable motion, the second- order Doppler error concerns not only velocity but also acceleration. When muzzle velocity is zero the second- order Doppler error caused by an acceleration of 0.6g can be up to 2.5 nm in 0.4 s, which is not negligible in nanometric measurement. Moreover, when the muzzle velocity is nonzero, the accelerated motion may result in a greater error and decelerated motion may result in a smaller error.

Optimal H2 order-one reduction by solving eigenproblems for polynomial equations

A method is given for solving an optimal H2 approximation problem for SISO linear time-invariant stable systems. The method, based on constructive algebra, guarantees that the global optimum is found; it does not involve any gradient-based search, and hence avoids the usual problems of local minima. We examine mostly the case when the model order is reduced by one, and when the original system has distinct poles. This case exhibits special structure which allows us to provide a complete solution. The problem is converted into linear algebra by exhibiting a finite-dimensional basis for a certain space, and can then be solved by eigenvalue calculations, following the methods developed by Stetter and Moeller. The use of Buchberger's algorithm is avoided by writing the first-order optimality conditions in a special form, from which a Groebner basis is immediately available. Compared with our previous work the method presented here has much smaller time and memory requirements, and can therefore be applied to systems of significantly higher McMillan degree. In addition, some hypotheses which were required in the previous work have been removed. Some examples are included.

Semi-implicit second order finite volume jet stream computations

Semi-implicit, second order temporal and spatial finite volume computations of the flow in a differentially heated rotating annulus are presented. For the regime considered, three cyclones and anticyclones separated by a relatively fast moving jet of fluid or "jet stream" are predicted. Two second order methods are compared with, first order spatial predictions, and experimental measurements. Velocity vector plots are used to illustrate the predicted flow structure. Computations made using second order central differences are shown to agree best with experimental measurements, and to be stable for integrations over long time periods (> 1000s). No periodic smoothing is required to prevent divergence.

Lift and the leading-edge vortex

Flapping wings often feature a leading-edge vortex (LEV) that is thought to enhance the lift generated by the wing. Here the lift on a wing featuring a leading-edge vortex is considered by performing experiments on a translating flat-plate aerofoil that is accelerated from rest in a water towing tank at a fixed angle of attack of 15°. The unsteady flow is investigated with dye flow visualization, particle image velocimetry (PIV) and force measurements. Leading-and trailing-edge vortex circulation and position are calculated directly from the velocity vectors obtained using PIV. In order to determine the most appropriate value of bound circulation, a two-dimensional potential flow model is employed and flow fields are calculated for a range of values of bound circulation. In this way, the value of bound circulation is selected to give the best fit between the experimental velocity field and the potential flow field. Early in the trajectory, the value of bound circulation calculated using this potential flow method is in accordance with Kelvin's circulation theorem, but differs from the values predicted by Wagner's growth of bound circulation and the Kutta condition. Later the Kutta condition is established but the bound circulation remains small; most of the circulation is contained instead in the LEVs. The growth of wake circulation can be approximated by Wagner's circulation curve. Superimposing the non-circulatory lift, approximated from the potential flow model, and Wagner's lift curve gives a first-order approximation of the measured lift. Lift is generated by inertial effects and the slow buildup of circulation, which is contained in shed vortices rather than bound circulation. © 2013 Cambridge University Press.

亚硝化—厌氧氨氧化处理猪场废水研究

猪场废水COD浓度高、氨氮浓度高、悬浮物浓度高，已成为农村面源污染的主要来源，并严重威胁到农村饮用水安全。猪场废水氨氮浓度高、处理难度大，如何采用经济高效的方法，去除氨氮使其达到排放标准，一直是猪场废水处理中面临的重要难题。 厌氧氨氧化是近年受到国内外水处理研究者广泛关注的新型生物脱氮技术，具有不需要外加有机碳源、节省供氧量、降低能耗等优点。虽然国内外研究者对厌氧氨氧化过程的脱氮机理、厌氧氨氧化菌的生理生化特性等进行了多方面的研究，但已有的报道大多以模拟废水为研究对象，以猪场废水为研究对象的报道，在国内外文献中极少有报导。 本论文以猪场废水为主要研究对象，考察了猪场废水的亚硝化过程、厌氧氨氧化的启动过程，并对亚硝化和厌氧氨氧化联合用于猪场废水脱氮进行了探索。 1．论文首先研究了猪场废水的亚硝化过程，考察了废水水质和主要运行条件对亚硝化过程的影响。实验表明：（1）亚硝化阶段反应时间为8到10h时，出水中氨氮和亚硝酸盐浓度比可达到1：1～1：1.23，满足厌氧氨氧化反应对二者比例的要求；达到前述要求时，氨氮去除率达到58.3～65.6 %，亚硝化率在整个过程均保持在97 %以上，COD去除率在59.2～68.6 %；（2）曝气量（溶解氧）对亚硝化过程影响显著，随着曝气量增大，达到厌氧氨氧化要求的氨氮与亚硝酸盐氮浓度比例所需水力停留时间τ越短，pH出现明显下降的时间越短；（3）τ对应的pH在7.8～8.1之间，无需进行pH调节即可满足厌氧氨氧化反应对pH的要求；（4）氨氮和COD降解过程遵循一级反应动力学，氨氮和COD降解的速率常数分别为0.0656～0.0724 1/h和0.0491～0.0664 1/h。 2．在进行亚硝化过程研究的同时，以模拟废水为试验对象，进行厌氧氨氧化启动研究。以反硝化污泥和养殖厂储水池厌氧底泥的混合污泥作为接种污泥，历时大约100天，培育出具有厌氧氨氧化活性的污泥，氨氮和亚硝酸盐氮最高进水浓度分别为223.8 mg/L和171.4 mg/L，去除率最高分别达48%和41.5%，此时二者消耗比例为1.33：1。 3．在猪场废水的亚硝化研究完成和厌氧氨氧化过程初步启动成功后，在模拟废水中逐步加入猪场废水的亚硝化处理出水，逐步实现亚硝化和厌氧氨氧化的组合。亚硝化出水添加到厌氧反应器后，厌氧氨氧化反应仍可继续进行，且去除效率逐步提高。研究发现添加的亚硝化出水中携带的亚硝化细菌在厌氧氨氧化菌膜外层生长并累积，增加了厌氧氨氧化反应基质的传质阻力，妨碍了厌氧氨氧化效率的提高。 4．亚硝化－厌氧氨氧化实际工程应用探索中，生物接触氧化池可在有效去除废水中的有机物的同时实现亚硝化，出水中氨氮和亚硝酸盐比例平均为1.10，可满足后续厌氧氨氧化的要求；在适宜的进水浓度和温度下，ABR池出现了厌氧氨氧化启动的迹象；研究同时发现，水质的波动和气温的变化是工程中影响厌氧氨氧化菌活性的重要因素。 论文的主要创新点在于：（1）以猪场废水为研究对象，以实现厌氧氨氧化为目标，对亚硝化过程进行了比较详细的考察，获得了亚硝化出水满足厌氧氨氧化要求的工艺条件，通过对其COD和氨氮降解过程的考察，得出亚硝化阶段COD降解和氨氮去除的动力学模型；（2）对亚硝化－厌氧氨氧化处理猪场废水进行了探索，确立了影响其污染物去除率稳定的重要因素。 论文的上述研究成果，为厌氧氨氧化技术的实用性研究提供理论依据。 Piggery wastewater, which is characterized by high concentration of COD、ammonium and suspend substance, has become a most important source of non-point source pollution and also severely threats drinking water security in rural area. How to discharge piggery wastewater with the ammonium concentration meeting standard by economical and effective method? This is the most urgent problem in piggery wastewater treatment. As a new biological nitrogen removal technology, Anammox process has been paid more and more attention by researchers all over the world. Anammox has advantages of no need of organic carbon addition, low oxygen consumption and energy consumption. Plenty of investigations have been carried out to the mechanism, physiological and biochemical characteristic of bacteria about Anammox. Most of researches focused on synthetic wastewater, there is rare report about its application in piggery wastewater. In this paper，experimental studies were performed to investigate Sharon process in treatment of piggery wastewater，the start up process of Annammox using synthetic wastewater were studied, the feasibility of applying Sharon-Anammox process in the nitrogen removal of piggery wastewater was evaluated. 1. Sharon process of piggery wastewater was firstly investigated to analyze the effects of water quality and main running parameters, which meet the NH4+-N to NO2--N ratio requirement of successive Anammox. Results showed: (1)During Sharon Process，after 8～10 hours’ reaction the NH4+-N to NO2--N ratio in effluent reached 1:1.0～1:1.23, when the removal percentage of NH4+-N was 58.3～65.6 %, a semi-nitration rate of above 97 % was achieved during the process; meanwhile 59.2～68.6 % of the COD was also removed. (2)The aeration rate（oxygen） had obvious effect on the hydraulic retention time(τ) which met the NH4+-N to NO2--N ratio requirement of Anammox. As aeration rate increased, the hydraulic retention time(τ) was shortened. (3) The pH corresponding to τ was between 7.8 and 8.1, thus it needed no artificial adjustment. (4) The reduction of ammonia and COD followed the first-order reaction kinetics. The velocity constants of ammonia and COD were 0.0656～0.0724 1/h and 0.0491～0.0664 1/h, respectively. 2. The startup of Anammox process using the artificial wastewater was performed simultaneously with Sharon. The aim was to investigate the running parameters of Anammox and make foundation for the combination stage. By using the mixture of denitrifying sludge and anaerobic sludge in tank of the breeding factory, sludge of Anammox activity was cultivated in UASB after 100 days. The removal percentage of NH4+-N and NO2-N were up to 48% and 41.5%, respectively, when the NH4+-N and NO2-N influent concentration were 223.8 mg/L and 171.4 mg/L, respectively, the NH4+-N and NO2-N removal rate was 1.33:1. 3. After investigation of Sharon and startup of Anammox, effluent of Sharon process was added into the synthetic wastewater to combine Sharon and Anammox step by step. It took some time after the addition of Sharon effluent that Anammox reaction continued and the removal rate kept increasing. It indicated that nitrifying bacteria were carried by the Sharon effluent cumulated in the outer layer of Anammox. This enhanced transfer resistance of Anammox reaction and the increasing removal rate was restrained. 4. In the bio-contact oxidation pond of practical project, Sharon process were carried out successfully and organic compounds were removed effectively. An average NO2-N/ NH4+-N rate of 1:1.0 was achieved in the effluent, which met the requirement of successive Anammox. Under condition of suitable influent concentration and temperature, there was evidence that Anammox could start up in ABR. The variety of wastewater and temperature had great affects on Anammox activity in practical engineering. Innovation of this paper: (1) The Sharon process for treating piggery wastewater was discussed in details. Technological parameters that met requirement of Anammox were obtained. The dynamic models of COD and ammonium removal in the process were educed. (2) Sharon-Ananmmox for treatment of piggery wastewater was investigated, and the primary influencing factors was studied. This paper could be a theoretical consult for research of Anammox utility.

功能菌剂的混合发酵研究

随着化工行业的发展，大量有毒有害难降解有机物随工业废水的排放进入环境，这些物质能够在环境中长期存在、积累和扩散，通过食物链对动植物的生存及人类的健康造成不良影响。本文以苯酚、对氯硝基苯、氯苯和十六烷为模拟污染物，以前期研制的功能菌剂为对象，经过紫外线线诱变筛选出优于出发菌株的功能菌，对诱变后功能菌的理化性能进行了研究，对菌种进行了鉴定，在此基础上，就其相互之间的微生态关系进行研究，为混合发酵提供理论基础，并就其最佳发酵条件及发酵参数进行了研究，最后对发酵产品的性能进行了检测。目前，国内外有关功能菌剂的研究还存在多方面的不足，主要包括：①由于多菌种混合发酵过程较为复杂，各菌之间存在复杂的相互作用，影响因素较多，关于菌种之间的相互关系研究得很少，环境功能菌剂的发酵方法大多采用单独发酵后混合的方式。单独发酵对原材料、设备和能源的利用率较低，对于多菌种制剂发酵，在设备、能源和原材料的方面造成的浪费更大，将会大幅增加菌剂的生产成本，影响多菌种功能菌剂的发展；②功能菌剂生产过程的质量控制方面研究得较少；③功能菌剂产品的稳定性、抗冲击性能研究得较少，对环境微生物制剂的研究主要集中在菌种选育和培养条件优化方面。 通过本论文研究，得到以下主要结论。 （1）在紫外线诱变处理中，用紫外线对发生一定程度退化的出发菌株进行诱变处理后，六株具有高效降解性能的菌株被筛选出来，诱变筛选出的菌株形态和ERIC-PCR指纹图谱与出发菌株相比发生了明显改变；而且诱变后的菌株对目标难降解底物的降解能力均得到改善，其中，FPN、FCB、F14、FEm对目标底物的降解率提高了20%以上；诱变后菌株经过7次连续传代接种后，对目标难降解底物的降解率无显著变化，具有一定的遗传稳定性。并对诱变后的功能菌进行了初步的鉴定，这6株菌都分别是芽孢杆菌。 （2）对诱变后的功能菌相互之间的微生态关系进行了研究，通过抑菌实验、生长量以及基质消耗量的比较，确定它们之间的生长关系是无害共栖关系，可以进行混合发酵。 （3）对该功能菌剂进行发酵培养条件研究，结果表明发酵培养基的最佳成分（g/L）：葡萄糖 31.0g/L、玉米粉10.0g/L、磷酸氢二钾1.0g/L、硫酸铵1.1g/L、硫酸镁0.55g/L。通过研究不同的培养条件对菌体生长和降解性能的影响，确定了最佳培养条件：培养基初始pH7.5；最适温度32℃；培养基装液量125mL（250 mL三角瓶），以及培养时间对降解性能的影响，培养20 h的产物对降解最为有利。通过研究添加不同目标污染物对菌体生长和降解性能的影响，确定了添加目标污染物的最佳量以及最佳时间：苯酚投加量：1.125 g/L，对氯硝基苯投加量：0.1 g/L；最佳投加时间为发酵培养开始后4 h。 （4）以摇瓶分批发酵最优条件为基础，对FPN、F10、FCB、FNa、F14 和 FEm进行了摇瓶分批发酵试验。以摇瓶分批发酵试验数据为依据，对功能菌剂分批发酵动力学进行了研究，建立了菌体生长和基质消耗的动力学模型，拟合模型能较好的反映功能菌剂分批发酵过程。 （5）功能菌剂和活性污泥协同作用，可以提高系统的生物降解能力，功能菌剂投加量为2%，新鲜活性污泥3500 mg/L，降解24 h条件下，功能菌剂和活性污泥的协同作用对COD的去除率和对照组相比，最多的提高了36.8%。功能菌剂和活性污泥协同作用以及活性污泥的单独作用，其生物降解过程均符合一级反应动力学过程，功能菌剂和活性污泥协同作用的生物降解动力学方程为：，相关系数97%。采用SBR运行方式，引入功能菌剂的SBR系统明显能够改善和提高生物降解的效率。与仅有活性污泥的系统相比，系统对COD的平均去除率可以提高27.1%，同时，系统的耐负荷冲击以及耐毒害冲击的性能比仅有活性污泥的SBR系统强，特别是负荷冲击对引入功能菌剂的SBR系统影响很小。仅有活性污泥的SBR系统经过负荷冲击和毒害冲击之后，不能恢复到冲击之前的水平，而且系统有效作用时间的周期比引入功能菌剂的SBR系统相比大大缩短，而引入功能菌剂的SBR系统处理效果较为稳定，恢复能力很强。 Along with the development of industries, many recalcitrant organic chemicals have been discharged into natural environments together with wastewaters and can exist in waters, soil and sediments for a long time without degradation. These haz-ardous substances, their byporducts and metabolizabilities can be highly toxic, mu-tagenic and carcinogenic, thereby threatening animals, plants and human health through food chain. Consequently the removal of these compounds is of significant interest in the area of wastewater treatment. In this dissertation, the phenol, hydro-quinone, chlorobenzene and hexadecane treated as the model pollutants, the func-tional microorganism agent was used as the starting strains, they treated with ultra-violet light, and then the mutant strains with high degradation ability were screened out and identified primarily, the relationship between these stains were studied, the medium composition and fermentation conditions were optimized, the degradation ability of the fermented production was tested. The literature survey indicates that the study of the microorganism agent is far from complete and more information is re-quired on following problems. 1, Because of the complexity of relationship in mixed fermentation and the complicated factors, the study is hardly to process.2, There is a lack of information on the quality control of the producing process .3, And there is a lack of information on the stability about the microorganism agent. In this dissertation, the main results of the present study could be summarized as follows: (1)The degenerate starting strains were treated with the ultraviolet light, and six mutant strains with high biodegradation ability were screened out by using the me-dium with selective pressure of model pollutants. The mutant strains had great changes in colonialmorphology and ERIC-PCR fingerprinting. And the mutant strains got obvious advantages over the starting strains in degradation ability and over 20% improvement of removal rates was achieved for FPN、FCB、F14 and FEm. The de-gradation ability of the mutant strains was stable after seven generations. After that, the mutant strains were primarily identified as bacillus respectively. (2) The relationship between these mutant strains was studied. By the compari-son of antibiosis effect, biomass and consumption of substrate, the relationships were neutralism and they could be mixed fermented. (3) The optimized cultivation conditions were as follows: glucose 31.0 g/L, corn power 10 g/L, K2HPO4 1.0 g/L, (NH4)2SO4 1.1 g/L, MgSO4 0.55 g/L, initial pH7.5, temperature 32℃, working volume 125 mL/250 mL, and cultivation time 20h (con-sidering the time effect on degradation ability), adding pollutants phenol (1.125 g/L) and hydroquinone (0.1 g/L) into the broth at 4 h after cultivation. (4) Based on the above optimum condition, the batch fermentation was per-formed with strains FPN, F10, FCB, FNa, F14 and FEm in shake flask. The batch fermentation kinetics was studied based on the experimental data. Two kinetic models were constructed which could reflect the regularity of growth and substrate consump-tion in the process of batch fermentation. (5) The co-operation of functional microorganism agent and activated sludge could raise biodegradation of system by adding some microorganism agent and 3500 mg/L fresh activated sludge. Bioaugumentation by the addition of high effective deg-radation culture enhanced the treatment effect of SBR system and the COD removal rate was increased by 20%-36.8%. Its biodegradation matched first-order dynamical reaction equation, and the reaction equation was ln0.2327.391ct=−+. The micro-organism agent had the effect of optimization to activated sludge micro-ecosystem. The SBR system adding 2% microorganism agent, the average COD removal rate of that was increased by 27.1% and stronger anti-shock ability to load and toxicant were achieved (compared with SBR system just adding activated sludge). Especially the load-shock has barely effect to the SBR system adding microorganism agent. After the load and toxicant shock, the SBR system just adding activated sludge couldn’t come back to original level and the activated sludge micro-ecosystem was frustrated. The applying of microorganism agent increased biological activity and system’s re-sistance ability to load shock and toxicant shock.

Chemisorption and reaction characteristics of methyl radicals on Cu(110)

Methyl radicals are generated by pyrolysis of azomethane, and the condition for achieving neat adsorption on Cu(110) is described for studying their chemisorption and reaction characteristics. The radical-surface system is examined by X-ray photoemission spectroscopy, ultraviolet photoemission spectroscopy, temperature-programmed desorption, low-energy electron diffraction (LEED), and high-resolution electron energy loss spectroscopy under ultrahigh vacuum conditions. It is observed that a small fraction of impinging CH3 radicals decompose into methylene possibly on surface defect sites. This type of CH2 radical has no apparent effect on CH3(ads) surface chemistry initiated by dehydrogenation to form active CH2(ads) followed by chain reactions to yield high-mass alkyl products. All thermal desorption products, such as H-2, CH4, C2H4, C2H6, and C3H6, are detected with a single desorption peak near 475 K. The product yields increase with surface coverage until saturation corresponding to 0.50 monolayer of CH3(ads). The mass distribution is, however, invariant with initial CH3(ads) coverage, and all desorbed species exhibit first-order reaction kinetics. LEED measurement reveals a c(2 x 2) adsorbate structure independent of the amount of gaseous exposure. This strongly suggests that the radicals aggregate into close-packed two-dimensional islands at any exposure. The islanding behavior can be correlated with the reaction kinetics and is deemed to be essential for the chain propagation reactions. Some relevant aspects of the CH3/Cu(111) system are also presented. The new results are compared with those of prior studies employing methyl halides as radical sources. Major differences are found in the product distribution and desorption kinetics, and these are attributed to the influence of surface halogen atoms present in those earlier investigations.

Second-order Stokes solutions for internal waves in three-layer density-stratified fluid

In this paper, internal waves in three-layer stratified fluid are investigated by using a perturbation method, and the second-order asymptotic solutions of the velocity potentials and the second-order Stokes solutions of the associated elevations of the interfacial waves are presented based on the small amplitude wave theory. As expected, the first-order solutions are consistent with ordinary linear theoretical results, and the second-order solutions describe the second-order modification on the linear theory and the interactions between the two interfacial waves. Both the first-order and second-order solutions derived depend on the depths and densities of the three-layer fluid. It is also noted that the solutions obtained from the present work include the theoretical results derived by Umeyama as special cases.

Representation of microscale ocean wavenumber spectrum correct to the second order

In this paper, the analytical representations of four wave source functions in high-frequency spectrum range are given on the basis of ocean wave theory and dimensional analysis, and the perturbation method is used to solve the governing equations of ocean wave high-frequency spectrum on the basis of the temporally stationary and locally homogeneous scale relations of microscale wave. The microscale ocean wavenumber spectrum correct to the second order has an explicit structure, its first order part represents the equilibrium between different source functions, and its second order part represents the contribution of microscale wave propagation.

Second-order Stokes wave solutions for intefacial waves in three-layer stratified fluid with background current

In this paper, interfacial waves in three-layer stratified fluid with background current are investigated using a perturbation method, and the second-order asymptotic solutions of the velocity potentials and the second-order Stokes wave solutions of the associated elevations of the interfacial waves are presented based on the small amplitude wave theory, and the Kelvin-Helmholtz instability of interfacial waves is studied. As expected, for three-layer stratified fluid with background current, the first-order asymptotic solutions (linear wave solutions), dispersion relation and the second-order asymptotic solutions derived depend on not only the depths and densities of the three-layer fluid but also the background current of the fluids, and the second-order Stokes wave solutions of the associated elevations of the interfacial waves describe not only the second-order nonlinear wave-wave interactions between the interfacial waves but also the second-order nonlinear interactions between the interfacial waves and currents. It is also noted that the solutions obtained from the present work include the theoretical results derived by Chen et al (2005) as a special case. It also shows that with the given wave number k (real number) the interfacial waves may show Kelvin-Helmholtz instability.

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

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.

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

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

断层封闭性及其与油气成藏的关系研究-以沾化凹陷孤岛地区为例

Godao area, located in the east of the Zhanhua depression of Jirang sag in Bohai Bay Basin, is the studying area in my dissertation. It is first time that fault sealing properties and the related relationship with the pool forming are studied in Gudao area. On the base of the analysis of the regional tectonics, the author has studied the tertiary structural evolution of the Gudao area and distinguished the fault's level and put forward the distinguishing principle. The geometrical feature, mechanical characters, developmental mechanism and history of the boundary faults in the tectonic unit of this area are all studied and emphasized especially. The buried history of oil-generating depression (that is Gudao depression) and the history of oil and gas migration simultaneously are discussed, the juxtaposition relationship between boundary fault evolutionary history and oil and gas migrated history are expatiated. To the geological condition of the Gudao area, three level faults sealing properties of this area were discussed in detail. Their characteristics of behavior and the intrinsic relationship between their sealing and oil and gas migrated reservoir are elucidated. The pool-forming models related to fault seal are exposed. The author has studied the lithologies of different order of faults, the relationship of occurrence assemblage analysis, normal stress of fault plane in different depth and shale smear factor faults. Then analysis their role in the various faults sealing and confirms fault sealing marks of three different orders faults and exposes the mechanism of fault sealing. Shale smear zone formed by first order fault in lasting activities is one different type of fault breccia and mainly controlled factor to its entrapment of petroleum. Effective sealing threshold value and fault displacement is ascertained. Mainly controlled factor of second fault sealing is bigger compressive stress loaded on fault plane. According to this, quantitatively evaluated index is given. Shale smear zone is necessary condition for second fault stress entrapment. The juxtaposed relationship between the different lithologies within third order faults is most important controlled factor for its sealing. Based on various order of fault sealing features and mechanism in Gudao area, the author proposed three orders of fault sealing models. Shale smear zone sealing model, normal stress sealing model and lithologies juxtaposed sealing model are suggested to first, second and third order fault respectively. The conclusion of this studying has not only the very important theoretical significance and practical value in Gudo area but also the very important guiding role for other areas of related aspects.

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

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.

Oceanic controls on the primary production of the northwest European continental shelf: model experiments under recent past conditions and a potential future scenario

In this paper we clearly demonstrate that changes in oceanic nutrients are a first order factor in determining changes in the primary production of the northwest European continental shelf on time scales of 5–10 yr. We present a series of coupled hydrodynamic ecosystem modelling simulations, using the POLCOMS-ERSEM system. These are forced by both reanalysis data and a single example of a coupled ocean-atmosphere general circulation model (OA-GCM) representative of possible conditions in 2080–2100 under an SRES A1B emissions scenario, along with the corresponding present day control. The OA-GCM forced simulations show a substantial reduction in surface nutrients in the open-ocean regions of the model domain, comparing future and present day time-slices. This arises from a large increase in oceanic stratification. Tracer transport experiments identify a substantial fraction of on-shelf water originates from the open-ocean region to the south of the domain, where this increase is largest, and indeed the on-shelf nutrient and primary production are reduced as this water is transported on-shelf. This relationship is confirmed quantitatively by comparing changes in winter nitrate with total annual nitrate uptake. The reduction in primary production by the reduced nutrient transport is mitigated by on-shelf processes relating to temperature, stratification (length of growing season) and recycling. Regions less exposed to ocean-shelf exchange in this model (Celtic Sea, Irish Sea, English Channel, and Southern North Sea) show a modest increase in primary production (of 5–10%) compared with a decrease of 0–20% in the outer shelf, Central and Northern North Sea. These findings are backed up by a boundary condition perturbation experiment and a simple mixing model.