Numerical simulation of wavefronts diffracted by apertures with circular symmetry

The numerical simulation of the wavefronts diffracted by apertures with circular symmetry is realized by a numerical method. It is based on the angular spectrum of plane waves, which ignored the vector nature of light. The on-axial irradiance distributions of plane wavefront and Gauss wavefront diffracted by the circular aperture have been calculated along the propagation direction. Comparisons of the simulation results with the analytical results and the experimental results tell us that it is a feasible method to calculate the diffraction of apertures. (c) 2006 Published by Elsevier GmbH.

羟基自由基降解琼胶及其产物的活性研究

琼胶是一种从石花菜等红藻中提取的，目前生产工艺和结构等方面研究比较成熟的海藻多糖，广泛应用于医药、仪器等行业。但是，海藻多糖因为具有分子量大，粘度大，溶解度较小的等特点，而使其应用范围受到限制。利用降解的手段对其进行修饰，降低分子量和粘度，改善溶解性，可以拓展其应用范围。并且根据文献报道，琼 胶寡糖具有一些特殊的生物活性，如抗氧化性，抗炎症等。因此，对琼胶降解的研究具有生要意义。本研究中，为了选择一种合适的降解方法，进行了几种水解方法的尝试，其中包括在不同湿度和酸度下盐酸水解，过氧化氢和醋酸催化水解，Fenton体系羟基自由基降解。对于酸水解和Fenton体系氧化还原降解方法，通过粘度法对反应的速度进行了比较，表明氧化还原降解反应中琼胶的粘度降低比较快，并且具有代表性和新意，确定为本实验的降解琼胶的方法并对氧化还原降解所得的产物进行了活性实验。通过模仿自然界普遍存在的氧化还原降解反应，利用Vc诱导的Fenton体系产生的羟基自由基氧化还原降解琼胶得到低分子量的琼胶。降解产物经过高速离心、60％乙醇沉淀，除去分子量比较大的降解产物和磷酸盐，得到可溶于60％乙醇的分子量估计小于3000的降争产物，其产率为85％。利用经Sephadex-G25凝胶色谱分离所香的不同分子量的级分进行分子量和α－葡萄糖苷酶抑制活性关系的实验。降解产物对α－葡萄糖苷酶的抑制率和各级分的浓度呈线性正相关，并且各级分的IC_(50)则随着分子量的降低而降低。另外，对所得的降解产物混合物进行了红外吸收光谱、质子去偶核磁共震碳谱和负离子基质辅助激光诱导－飞行时间质谱结构分析。结果表明，氧化还原降解反应的专一性差，在得到寡糖的同时，在光谱图中出现一些比较复杂的副产物的结构信息。最后，根据MTT法的原理，以有体皮肤成纤维细胞为材料，通过紫外线辐射产生自由基造成氧化损伤，研究降解产物对成纤维细胞的保护作用。当无紫外线辐射时，降解产物对成纤维细胞具有显著的促进生长增殖作用：当经UVa、UBb辐射时则可以显著地表现出对损伤的保护作用，并且这种促进生长和保护作用呈显著的量效关系，表明降解产物具有清除基自由基的作用。但是，因为氧化还原降解以应的机理尚不十分明的以及琼羟胶的特殊结构，使得反应的副产物很难预测，也就使得分离工作难以进行，所以，根据目前所得的信息，尚不能确定是降解产物的什么级分产生的以上两种生物活性。

拟人机器人交互作用力的研究

拟人机器人的动力学具有高度非线性、高度耦合的特点,分析清楚各组成部分之间的交互作用力是实施高级控制方法的基础。文章在以往分析移动机械手的基础上,从整体建模的角度入手,对拟人机器人的交互作用力提出了一个新的模型,即神经网络模型。利用该模型对一个特殊的单一手臂运动的例子进行了拟合,其结果是收敛的。这说明提出的模型是有效的,此后,我们将陆续给出研究成果。

华北克拉通孔玆岩系中石榴子石花岗岩的成因：熔体与残留相混合机制的研究

The most widespread rock associations in the Western Block of North China Craton are khondalites distributed mainly in Jining, Liangcheng and Datong. A large quantitiy of garnet-bearing granites are contained in the khondalites. A great deal of research has been carried out on them by previous researchers. Studies of these garnet-bearing granites consist essentially of structural characteristics, petrography and geochemistry, and finally geochronological determinations. Summing up these researches, it will not be difficult to see that all of these authors have regarded these large numbers of garnets (up to 20%) contained in granites as crystallized products from magmas, but they have not proved this from petrological perspective. Theoretically, there are possibly three kinds of petrogenesis as to these garnets. The first one is that they have been transferred to the granites from khondalites by melt when anatexis happened to khondalites, and they, in essence, are residual metamorphic garnets; The second one is that when the khondalites were being melted, these garnets were produced from biotite dehydration melting, and the newly formed garnets intruded together with the melt and eventually molded the garnet-bearing granites. Garnets of this possible kind either showed independent crystals, or garnets from khondalites took place secondary growth under favorable temperature and pressure conditions for their crystallization; The last possibility is that these garnets were crystallized from magmas in which suitable pressure, temperature and composition were available. These garnets, generally, should be fine-grained. The aim of this study is, through examining the mineral chemistry of the garnets and the whole rock chemistry, to ascertain under which kind of mechanism, in the world, did these garnets form? Besides, we try to calculate the temperatures under which khondalites began melting and reactions of the garnets and the cooled melts happened by garnet-biotite thermometry. The whole rock chemistry analyses of the garnet-bearing granites tell us that all the samples are strongly peraluminous (A/CNK greater than 1.1) on the A/NK vs. A/CNK plot. On the SiO2-K2O plot, the granites are mainly constrained to be high-K calc-alkaline and calc-alkaline series, consistent with previous researches. On the ACF（（Al2O3-Na2O-K2O）-FeO(T)-CaO） discrimination plot, all the six garnet-bearing granite samples drop into the area of S-type granites. The relationship between CaO/Na2O and SiO2 shows that the overwhelming majority of garnet-bearing granites have a CaO/Na2O value over 0.3, revealing that they probably come from metagreywacke precursors or mediate-felsic orthogeneisses compositionally similar to them. Detailed EPMA analyses conducted on the garnets contained in the garnet-bearing granites show that all the garnets are dominated by almandine and pyrope, which occupy 92-96% (Weight Percentage) of each garnet analyzed, typical of granulite facies. Their chemical composition is entirely different from those crystallized in magmas, but extremely similar to those of typical granulite facies metapelites in khondalites and typical granulites, indicating all the garnets to be metamorphogenic. In addition, REEs distribution patterns of the garnets are totally different from typical biotite granites and peraluminous granites. In other words, both LREE and HREE of our garnets are evidently lower than those from these two kinds of rocks. Moreover, compared to the REE pattern of the garnets from typical amphibolites, LREE content of our garnets is obviously higher and HREE content is a little lower. However, REE patterns of our garnets are completely in harmony with those of garnets from typical granulites. So, the REE patterns of garnets, again, prove that all the garnets we studied are metamorphogenic. Biotites appear in two forms, being as inclusions in the garnet and as selvages immediately adjacent to the garnet, respectively. Two reactions and their corresponding temperatures, with the help of petrography and Garnet-Biotite geothermometers, could be obtained, which are Bt+ Pl+ Qtz→Kfs+ Opx+ Grt+ melt as positive reaction and Kfs+ Grt+ melt→Bt+ Pl+ Qtz as reverse reaction, respectively. Summing up the discussion above, we declare that the garnet-bearing granites distributed in the Western Block of North China Craton are the mixture of melts and restites resulted from biotite dehydration melting. The garnets contained in the restites are the products from biotite dehydration melting and restites from the khondalites, respectively.

反射地震转换波成像方法研究

Multi-waves and multi-component get more and more attentions from oil industry. On the basis of existent research results, My research focuses on some key steps of OBC 4C datum processing. OBC datum must be preprocessed quite well for getting a good image. We show a flow chart of preprocess including attenuation of noise on multi-component datum、elimination ghost by summing P and Z and rotation of horizontal components. This is a good foundation for the coming steps about OBC processing. How to get exact converted point location and to analyze velocity are key points in processing reflection seismic converted wave data. This paper includes computing converted point location, analyzing velocity and nonhyperbolic moveout about converted waves. Anisotropic affects deeply the location of converted wave and the nonhyperbolic moveout. Supposed VTI, we research anisotropic effect on converted wave location and the moveout. Since Vp/Vs is important, we research the compute method of Vp/Vs from post-stack data and pre-stack data. It is a part of the paper that inversing anisotropic parameter by traveltime. Pre-stack time migration of converted wave is an focus, using common-offset Kirchhoff migration, we research the velocity model updating in anisotropic media. I have achieved the following results: 1) using continued Fractions, we proposed a new converted point approximate equation, when the offset is long enough ,the thomsen’s 2 order equation can’t approximate to the exact location of converted point, our equation is a good approximate for the exact location. 2) our new methods about scanning nonhyperbolic velocity and Vp/Vs can get a high quality energy spectrum. And the new moveout can fit the middle and long offset events. Processing the field data get a good result. 3) a new moveout equation, which have the same form as Alkhalifah’s long offset P wave moveout equation, have the same degree preciseness as thomsen’s moveout equation by testing model data. 4) using c as a function of the ratio offset to depth, we can uniform the Li’s and thomsen’s moveout equation in a same equation, the model test tell us choice the reasonable function C can improve the exact degree of Li’s and thomsen’s equation. 5) using traveltime inversion ,we can get anisotropic parameter, which can help to flat the large offset event and propose a model of anisotropic parameter which will useful for converted wave pre-stack time migration in anisotropic media. 6）using our pre-stack time migration method and flow, we can update the velocity model and anisotropic parameter model then get good image. Key words: OBC, Common converted Point (CCP), Nonhyperbolic moveout equation, Normal moveout correction, Velocity analysis, Anisotropic parameters inversion, Kirchhoff anisotropic pre-stack time migration, migration velocity model updating

人工源宽角地震资料中的广角效应及其对地壳速度成像的影响

The seismic wide-angle reflection/refraction method is the one of the most effective method for probing the crustal and upper mantle structure. It mainly uses the wide-angle reflection information from the boundary in the crust and the top boundary of the upper mantle to rebuild the crust and upper mantle structure. Through analyzing the reflection and transmission coefficients of various incident waves on the interface, we think relative to the pre-critical angle reflection information the post critical angle reflection information that received by wide-angle seismic data exists a time-shift effect with the offset variation, and then it must cause the error for velocity analysis and structure image. The feature of the wide-angle seismic wave field of the fourteen representative crust columns tell us that the wide-angle effects in the different representative tectonic units for the interface depth and the interval velocity in crust. We studied the features of the wide-angle seismic wave field through building the crust model and inverse its travel time by GA method to know the wide-angle influence on crustal velocity image. At last we finished the data processing of the Tunxi-Wenzhou wide-angle seismic profile. The results are as following: (1) Through building crust model, we labeled the travel time for all the phases by ray tracing method and remove wide-angle effects method, it revealed the wide-angle effect exists in the seismic data. (2) The travel time inversion by GA method can tell us that the depth by traditional ray tracing method is shallower than the result by remove wide-angle effects method, the latter can recover the crust structure model in effect. (3) We applied the two method mentioned before to the fourteen representative crust columns in China. It indicates that the removed wide-angle effect method in travel time inversion is reasonable and effective. (4) The real data processing from Tunxi-Wenzhou wide-angle seismic profile give us the basic structure through the two ways. The main influence exhibits in the difference of the interval velocity of the curst, and the wide-angle effects in shallow interface are stronger than the deep interface.

边坡动力响应分析及应用研究

As a marginal subject, dynamic responses of slopes is not only an important problem of engineering geology (Geotechnical problem), but also of other subjects such as seismology, geophysics, seismic engineering and engineering seismic and so on. Owning to the gulf between different subjects, it is arduous to study dynamic responses of slopes and the study is far from ripeness. Studying on the dynamic responses of slopes is very important in theories as well as practices. Supported by hundreds of bibliographies, this paper systemically details the development process of this subject, introduces main means to analyze this subject, and then gives brief remarks to each means respectively. Engineering geology qualitative analysis is the base of slopes dynamic responses study. Because of complexity of geological conditions, engineering geology qualitative analysis is very important in slopes stability study, especially to rock slopes with complex engineering geology conditions. Based on research fruits of forerunners, this paper summarizes factors influencing slopes dynamic stability into five aspects as geology background, stratums, rock mass structure, and topography as well as hydrogeology condition. Based on rock mass structure controlling theory, engineering geology model of the slope is grouped into two typical classes, one is model with obvious controlling discontinuities, which includes horizontal bedded slope, bedding slope, anti-dip slope, slide as well as slope with base rock and weathered crust; the other is model without obvious controlling discontinuities, which includes homogeneous soil slope, joint rock mass slope. Study on slope failure mechanism under dynamic force, the paper concludes that there are two effects will appear in slope during strong earthquake, one is earthquake inertia force, the other is ultra pore pressure buildup. The two effects lead to failure of the slope. To different types of slope failure, the intensity of two effects acting on the slope is different too. To plastic flow failure, pore pressure buildup is dominant; to falling rock failure and toppling failure, earthquake inertia force is dominant in general. This paper briefly introduces the principle of Lagrangian element method. Through a lot of numerical simulations with FLAC3D, the paper comprehensively studies dynamic responses of slopes, and finds that: if the slope is low, displacement, velocity and acceleration are linear enlarging with elevation increasing in vertical direction; if the slope is high enough, displacement, velocity and acceleration are not linear with elevation any more, on the other hand, they fluctuate with certain rhythm. At the same time, the rhythm appears in the horizontal direction in the certain area near surface of the slope. The distribution form of isoline of displacement, velocity and acceleration in the section of the slope is remarkably affected by the slope angle. In the certain area near the slope surface, isoline of displacement, velocity and acceleration is parallel to the surface of the slope, in the mean time, the strike direction of the extreraum area is parallel to the surface of the slope too. Beyond this area, the isoline direction and the strike direction of the extremum area turn to horizontal with invariable distance. But the rhythm appearing or not has nothing to with the slope angle. The paper defines the high slope effect and the low slope effect of slopes dynamic responses, discusses the threshold height H^t of the dynamic high slope effect, and finds that AW is proportional to square root of the dynamic elastic moduli El P , at the same time, it is proportional to period Tof the dynamic input. Thus, the discriminant of H^t is achieved. The discriminant can tell us that to a slope, if its height is larger than one fifth of the wavelength, its response regular will be the dynamic high slope effect; on the other hand, its response regular will be the dynamic low slope effect. Based on these, the discriminant of different slopes taking on same response under the same dynamic input is put forward in this paper. At the same time, the paper studies distribution law of the rhythm extremum point of displacement, velocity and acceleration, and finds that there exists relationship of N = int among the slope height H, the number of the rhythm extremum VHlhro) point N and ffthre- Furthermore, the paper points out that if N^l, the response of the slope will be dynamic high slope effect; \fN

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

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.