The second anti-lipopolysaccharide factor (EsALF-2) with antimicrobial activity from Eriocheir sinensis

The anti-lipopolysaccharide factor CALF) is a small basic protein that can bind and neutralize lipopolysaccharide (LPS), mediating degranulation and activation of an intracellular coagulation cascade. In the present study, cDNA of the second Eriocheir sinensis ALF (designated as EsALF-2) was cloned and the full-length cDNA of EsALF-2 was of 724 bp, consisting of an open reading frame (ORF) of 363 bp encoding a polypeptide of 120 amino acids. The deduced amino acid of EsALF-2 shared 82% similarity with EsALF-1 from E. sinensis and about 53-65% similarity with ALFs from other crustaceans. The potential tertiary structures of EsALF-1 and EsALF-2 contained two highly conserved-cysteine residues to define the LPS binding site, but the N-terminal of EsALF-1 formed a single additional alpha-helix compared to EsALF-2, implying that EsALF-1 and EsALF-2 might represent different biological functions in E. sinensis. The mRNA transcript of EsALF-2 was detected in all examined tissues of healthy crabs, including haemocytes, hepatopancreas, gill, muscle, heart and gonad, which suggested that EsALF-2 could be a multifunctional molecule for the host immune defense responses and thereby provided systemic protection against pathogens. The mRNA expression of EsALF-2 was up-regulated after Listonelln anguillarum and Pichia pastoris challenge and the recombinant protein of EsALF-2 showed antimicrobial activity against L. anguillarum and P. pastoris. indicating that EsALF-2 was involved in the immune defense responses in Chinese mitten crab against L. anguillarum and P. pastoris. These results together indicated that there were abundant and diverse ALFs in E. sinensis with various biological functions and these ALFs would provide candidate promising therapeutic or prophylactic agents in health management and diseases control of crab aquaculture. (C) 2010 Elsevier Ltd. All rights reserved.

On SAR wind speed ambiguities and related geophysical model functions

Recent investigations show that normalized radar cross sections for C-band microwave sensors decrease under high wind conditions with certain incident angles instead of increase, as is the case for low to moderate wind speeds. This creates the problem of ambiguities in high wind speed retrievals from synthetic aperture radar (SAR). In the present work, four geophysical model functions (GMFs) are studied, namely the high wind C-band model 4 (CMOD4HW), C-band model 5 (CMOD5), the high wind vertical polarized GMF (HWGMF_VV), and the high wind horizontal polarized GMF (HWGMF_HH). Our focus is on model behaviours relative to wind speed ambiguities. We show that, except for CMOD4HW, the other GMFs exhibit the wind speed ambiguity problem. To consider this problem in high wind speed retrievals from SAR, we focus on hurricanes and propose a method to remove the speed ambiguity using the dominant hurricane wind structure.

Strongly nonlinear long interfacial waves in uniform basic flows and their dispersion relationship

In this paper, long interfacial waves of finite amplitude in uniform basic flows are considered with the assumption that the aspect ratio between wavelength and water depth is small. A new model is derived using the velocities at arbitrary distances from the still water level as the velocity variables instead of the commonly used depth-averaged velocities. This significantly improves the dispersion properties and makes them applicable to a wider range of water depths. Since its derivation requires no assumption on wave amplitude, the model thus can be used to describe waves with arbitrary amplitude.

Forms and functions of inorganic carbon in the Jiaozhou Bay sediments

Inorganic carbon forms and their influencing factors, mutual transformation and contribution to carbon cycling in the Jiaozhou Bay sediments were discussed. The results show that inorganic carbon in sediments could be divided into five forms: NaCl form, NH3 center dot H2O form, NaOH form, NH2OH center dot HCl form and HCl form. Thereinto, NH2OH center dot HCl form and HCl form account for more than 70% of total inorganic carbon. There was close relationship among every form of inorganic carbon and their correlativity was clearly different with different sedimentary environment except the similar strong positive correlation among NH2OH center dot HCl form, HCl form and total inorganic carbon in all regions of the Jiaozhou Bay. All forms of inorganic carbon were influenced by organic carbon, pH, Eh, Es, nitrogen and phosphorus in sediments, but their influence had different characteristics in different regions. Every farm of inorganic carbon transformed into each other continuously during early diagenesis of sediments and the common phenomenon was that NaCl form, NH3 center dot H2O form, NaOH form and NH2OH center dot HCl form might transform into steady HCl form. NaCl form, NH3 center dot H2O form, NaOH form and NH2OH center dot HCl form could participate in carbon recycle and they are potential carbon source; HCl form may be buried for a long time in sediments, and it may be one of the final resting places of atmospheric CO2. Inorganic carbon which entered into sediments was about 4.98 x 10(10) g in the Jiaozhou Bay every year, in which about 1.47x10(10) g of inorganic carbon might be buried for a long time and about 3.51. x 10(10) g of inorganic carbon might return into seawater and take part in carbon recycling.

基于多线程的Linux下并发服务器的实现研究

探讨了Linux平台下多线程技术和套接字网络通讯问题,在此基础上利用互斥锁和条件变量技术设计了一个面向连接的多线程并发服务器的详细算法。最后给出了基于Posix线程库的Linux系统下用C++实现多线程并发服务器的基本程序框架。实际项目应用表明这种基于多线程的并发服务器结构在完善程序功能的同时还可以有效提高其服务性能。

强化学习算法中启发式回报函数的设计及其收敛性分析

回报函数设计的好与坏对学习系统性能有着重要作用,按回报值在状态-动作空间中的分布情况,将回报函数的构建分为两种形式:密集函数和稀疏函数,分析了密集函数和稀疏函数的特点.提出启发式回报函数的基本设计思路,利用基于保守势函数差分形式的附加回报函数,给学习系统提供更多的启发式信息,并对算法的最优策略不变性和迭代收敛性进行了证明.启发式回报函数能够引导学习,加快学习进程,从而可以实现强化学习在实际大型复杂系统应用中的实时控制和调度.

公路电子收费系统中DSRC协议规范的应用设计

根据领域中专用短程通信()协议的基础规范,设计出应用于电子收费系统的物理层、数据链路层和应用层相关参数ITSDSRCDSRC选定和设置。力求在电子收费领域相关的装置与设备能遵循指定的规格标准,并增进各系统相互之间的相容性和互连性。

空间电流体系地磁反演及理论研究

The space currents definitely take effects on electromagnetic environment and also are scientific highlight in the space research. Space currents as a momentum and energy provider to Geospace Storm, disturb the varied part of geomagnetic field, distort magnetospheric configuration and furthermore take control of the coupling between magnetosphere and ionosphere. Due to both academic and commercial objectives above, we carry on geomagnetic inverse and theoretical studies about the space currents by using geomagnetic data from INTERMAGNET. At first, we apply a method of Natural Orthogonal Components (NOC) to decomposition the solar daily variation, especially for (solar quiet variation). NOC is just one of eign mode analysis, the most advantage of this method is that the basic functions (BFs) were not previously designated, but naturally came from the original data so that there are several BFs usually corresponding to the process really happened and have more physical meaning than the traditional spectrum analysis with the fixed BFs like Fourier trigonometric functions. The first two eign modes are corresponding to the and daily variation and their amplitudes both have the seasonal and day-to-day trend, that will be useful for evaluating geomagnetic activity indices. Because of the too strict constraints of orthogonality, we try to extend orthogonal contraints to the non-orthogonal ones in order to give more suitable and appropriate decomposition of the real processes when the most components did not satisfy orthogonality. We introduce a mapping matrix which can transform the real physical space to a new mathematical space, after that process, the modified components which associated with the physical processes have satisfied the orthogonality in the new mathematical space, furthermore, we can continue to use the NOC decomposition in the new mathematical space, and then all the components inversely transform back to original physical space, so that we would have finished the non-orthogonal decomposition which more generally in the real world. Secondly, geomagnetic inverse of the ring current’s topology is conducted. Configurational changes of the ring current in the magnetosphere lead to different patterns of disturbed ground field, so that the global configuration of ring current can be inferred from its geomagnetic perturbations. We took advantages of worldwide geomagnetic observatories network to investigate the disturbed geomagnetic field which produced by ring current. It was found that the ring current was not always centered at geomagnetic equator, and significantly deviated off the equator during several intense magnetic storms. The deviation owing to the tilting and latitudinal shifting of the ring current with respect to the earth’s dipole can be estimated from global geomagnetic survey. Furthermore those two configurational factors which gave a quantitative description of the ring current configuration, will be helpful to improve the Dst calibration and understand the dependence of ring current’s configuration on the plasma sheet location relative to the equator when magnetotail field warped. Thirdly, the energization and physical acceleration process of ring current during magnetic storm has been proposed. When IMF Bz component increase, the enhanced convection electric field drive the plasma injection into the inner magnetosphere. During the transport process, a dynamic heating is happened which make the particles more ‘hot’ when the injection is more deeply inward. The energy gradient along the injection path is equivalent to a kind of force, which resist the plasma more earthward injection, as a diamagnetic effect of the magnetosphere anti and repellent action to the exotically injected plasma. The acceleration efficiency has a power law form. We use analytical way to quantitatively describe the dynamical process by introducing a physical parameter: energization index, which will be useful to understand how the particle is heated. At the end, we give a scheme of how to get the from storm time geomagnetic data. During intense magnetic storms, the lognormal trend of geomagnetic Dst decreases depend on the heating dynamic of magnetosphere controlling ring current. The descending pattern of main phase is governed by the magnetospheric configuration, which can be describled by the energization index. The amplitude of Dst correlated with convection electric field or south component of the solar wind. Finally, the Dst index is predicted by upstream solar wind parameter. As we known space weather have posed many chanllenges and impacts on techinal system, the geomagnetic index for evaluating the activity space weather. We review the most popular Dst prediction method and repeat the Dst forecasting model works. A concise and convnient Key Points model of the polar region is also introduced to space weather. In summary, this paper contains some new quantitative and physical description of the space currents with special focus on the ring current. Whatever we do is just to gain a better understanding of the natural world, particularly the space environment around Earth through analytical deduction, algorithm designing and physical analysis, to quantitative interpretation. Applications of theoretical physics in conjunction with data analysis help us to understand the basic physical process govering the universe.

南海沉积物中的放射虫及其古环境意义

In this paper a systematic study of radiolarian from surface sediments of all the South China Sea (SCS) has initially been done for its deposit ecology, biogeography and paleoenvironment significance. The paleoenvironment information obtained by radiolarian analysis and other sedimentary method for core samples is also made use as the synthesis proxy for revealing the paleoenvironment changes in the SCS and the relations of it with the past globe change during last 200ka. Some results come out of this study as: 1) Radiolarian skeleton chemistry composition and skeleton morphological features were analyzed, chiefly dividing them as 15 types of basic morphological features: 2) Analysis of biogeographical feature demonstrates that the fauna in SCS obviously belonging to a transitional type of west Pacific - Indian Ocean and has a particularity of itself: 3) Tendency of radiolarian population distribution is lower in shallow continental shelf area, increasing gradually toward the abyssal region; 4) Nine character boundaries of radiolarian depth distribution in the sediments from whole the South China Sea may be recognized; 5)Two radiolarian transfer functions for paleotemperature and paleo-primary productivity in the South China Sea have also been established respectively. The equation of transfer function for winter paleotemperature has only an average estimated error of 0.18678 ℃ and the equation for paleo-primary productivity has the calculation accuracy of 85.31%; 6) Changes of radiolarian individuals, numbers of species and H(S) values in core NS93-5 show the completely different oceanic geographical circumstance and ecology structure in the Last Maximum Glacial with present; 7) The abundance variation of some raiolarian warm species and cold species indicate the changes of water masse features along with the paleoenvironment evolution, showing that this sea area clearly is controlled in 6 issues of oxygen isotopes by the cold water masses; 8) By comparative analysis of δ ~(18)O curves with GISP2;s ice core can core 17940 of the northern SCS, the occurring characters of D/O's events 1-21 and Heinrich's events H1-H6 have been revealed in this sea area by core NS93-5, which prove the existence of paleoclimatic tele-connections between the southern SCS and Arctic region since about 200ka BP.

中国东部上第三系硬粘土的工程地质特性及边坡开挖工程环境效应研究

The engineering geological properties of Neogene hard clays and related engineering problems are frontiers in the fields of Engineering Geology, Soil Mechanics and Rock Mechanics. Recently, it has been recognized that Neogene hard clay is the intermediate type of material between the soil and the rock. Many aspects of them, such as sampling, testing, calculating and engineering process, are special, which could not be researched by means of theories and methods of traditional Soil Mechanics of Rock Mechanics. In order to get real knowledge and instruct the engineering practice, intersect studying of multiple disciplines, including Engineering Geology, Soil Mechanics and Rock Mechanics, etc., is necessary. Neogene hard clay is one of the important study objects of regional problem rocks & soils in our country, which extensively distributed in China, especially in Eastern China. Taking the related areas along the middle line of the Project of Transferring Water from the South to the North (e.g. Nanyang basin, Fangcheng-Baofeng area and Handan-Yongnian area), South-west of Shandong, Xu-Huai area and Beijing area, etc. as main study areas, the paper divided Neogene hard clays into reduction environment dominated origin and oxidation environment dominated origin, which distributed on areas western and eastern to Mount Taihangshan respectively. Intermediate types are also existed in some areas, which mainly distribute near the edges of depositional basins; they are usually of transitions between diluvial and lacustrine deposits. As to Neogene hard clays from Eastern China, the clay particle content is high, and montmorillonite or illite/montmorillonite turbostratic mineral is the dominating clay mineral. The content of effective montmorillonite is very high in each area, which is the basis for the undesirable engineering properties of Neogene hard clays. For hard clays from the same area, the content of effective montmorillonite in gray-greenish hard clay is much higher than that in purple-brownish or brown-yellowish hard clay, which is the reason why the gray-greenish hard clay usually has outstanding expansive property. On the other hand, purple-brownish or brown-yellowish hard clay has relatively less montmorillonite, so its property is better. All of these prove that the composition (clay mineral) of Neogene hard clay is the control factor for the engineering properties. Neogene hard clays have obvious properties such as fissured, overconsolidated and expansive, which are the main reasons that many engineering problems and geological harzards usually occur in Neogene hard clays. The paper systematically elaborates the engineering properties of Neogene hard clays from Eastern China, analyses the relationships between engineering properties and basic indexes. The author introduces the ANN method into the prediction of engineering property indexes of hard clays, which provides a new way for quantitatively assessment and prediction of engineering property indexes. During investigation in the field, the author found that there exists obvious seam-sheared zone between different hard clays in Miocene Xiacaowan formation in Xu-Huai area. Similar phenomenon also exists near the borderline between Neogene hard clays and underlying coal measures in the Southwest of Shandong province, which could be observed in the cores. The discovery of seam-sheard zone has important theoretical and practical significance for engineering stability analysis and revealing the origin of fissures in Neogene hard clays. The macrostructure, medium structure and microstructure together control the engineering properties of hard clays. The author analyses and summarizes the structural effects on hard clays in detail. The complex of the strength property of hard clays is mostly related to the characteristics of fissures, which is one of the main factors that affect the choice of shear strength parameters. So structure-control theory must be inseparably combined with composition-control theory during the engineering geological and rock/soil mechanics research of hard clays. The engineering properties, such as fissured, overconsolidated and expansive, control the instability of engineering behaviors of Neogene hard clays under the condition of excavation, i.e. very sensitive to the change of existence environment. Based on test data analysis, the author elaborates the effects of engineering environment change on the engineering properties. Taking Nanyang basin as example, the author utilizes FEM to study the effects of various factors on stability of cutting canal slopes, than sets forth the characteristics, development laws and formation mechanism of the deformation and failure of hard clay canal slopes, summarizes the protection and reinforcement principles, as well as the protection and remedy steps. On the basis of comparison of engineering properties of domestic and foreign Neogene muddy deposits, in the view of whole globe and associated with the geological characteristics of China, the paper demonstrates that the intermediate type of the material between the soil and the rock, named "hard clay/soft rock", which can not be separated abruptly, really exists in China. The author has given a preliminary classification based on its geological origin and distribution law, which is very significant for promoting the mixture of Engineering Geology, Soil Mechanics and Rock Mechanics. In the course of large scales engineering construction in China, many engineering experiences and testing data are gained, summarizing these testing results and automatically managing them with computer technology are very necessary. The author develops a software named "Hard Clay-Soft Rock Engineering Geological Information Management and Analysis System (HRGIMS)", realizes the automatic and visual management of geo-engineering information, on the basis of information management, the functions of test data analysis and engineering property prediction are strengthened. This system has well merits for practice and popularization.

孤东油田馆陶组河流相储层流动单元模型与剩余油分布规律研究

Since C.L. Hearn presented the concept of flow unit in 1984, its generation mechanisms and controlling factors have been studied in many aspects using different methods by researchers. There are some basic methods to do the research, and there are several concepts and classification standards about flow unit. Based on previous achievements and using methodologies from sedimentary geology, geophysics, seismic stratigraphy, and reservoir engineering, the author systemically studies the factors controlling flow unit, puts forward a series of methods for recognition, classification and evaluation of flow unit. The results obtained in this paper have important significance not only for understanding the flow unit, but also for revealing the distribution of remaining oil. As a case, this paper deals with the reservoir rocks in Guantao Group of Gudong Oilfield. Zhanhua Sag, Jiyang Depression in Bohaiwan Basin. Based on the study of stratigraphic, depositional and structural characteristics, the author establishes reservoir geological models, reveals the geological characteristics of oil-bearing reservoir of fluvial facies, points out the factors controlling flow unit and geological parameters for classification of flow unit. and summarizes methods and technologies for flow unit study when geological, well-logging and mathematical methods are used. It is the first attempt in literatures to evaluate reservoir by well-logging data constrained by geological conditions, then a well-logging evaluation model can be built. This kind of model is more precise than ever for calculating physical parameters in flow unit. In a well bore, there are six methods to recognize a flow unit. Among them, the activity function and intra-layer difference methods are the most effective. Along a section, the composition type of flow unit can be located according amplitude and impedance on seismic section. Slice method and other methods are used to distinguish flow unit. In order to reveal the distribution laws of flow unit in space, the author create a new method, named combination and composition of flow unit. Based on microscopic pore structure research, the classification methods of flow unit are developed. There are three types of flow unit in the reservoir of fluvial facies. They have their own lithology, petrophysics and pore structure character. Using judgement method, standard functions are built to determine the class of flow unit of fluvial facies. Combining reservoir engineering methods, the distribution laws of remaining oil in different types, or in different part of a flow unit are studied. It is evident that the remaining oil is controlled by the type of flow unit. The author reveals the relationship between flow unit and remaining oil distribution, builds the flowing models, predicts the variation of reservoir parameters in space, put forward different methods developing remaining oil in different flow unit. Especially, based on the results obtained in this paper, some suggestions for the adjustment of the developing flow units have been applied in Districts No.4 and No.7, and good results have been yielded. So, the results of this paper can guide oil field development. They are useful and significant for developing the remaining oil and enhancing the oil recovery efficiency.

小波域地震资料的处理方法（WSSP）研究

The modeling formula based on seismic wavelet can well simulate zero - phase wavelet and hybrid-phase wavelet, and approximate maximal - phase and minimal - phase wavelet in a certain sense. The modeling wavelet can be used as wavelet function after suitable modification item added to meet some conditions. On the basis of the modified Morlet wavelet, the derivative wavelet function has been derived. As a basic wavelet, it can be sued for high resolution frequency - division processing and instantaneous feature extraction, in acoordance with the signal expanding characters in time and scale domains by each wavelet structured. Finally, an application example proves the effectiveness and reasonability of the method. Based on the analysis of SVD (Singular Value Decomposition) filter, by taking wavelet as basic wavelet and combining SVD filter and wavelet transform, a new de - noising method, which is Based on multi - dimension and multi-space de - noising method, is proposed. The implementation of this method is discussed the detail. Theoretical analysis and modeling show that the method has strong capacity of de - noising and keeping attributes of effective wave. It is a good tool for de - noising when the S/N ratio is poor. To give prominence to high frequency information of reflection event of important layer and to take account of other frequency information under processing seismic data, it is difficult for deconvolution filter to realize this goal. A filter from Fourier Transform has some problems for realizing the goal. In this paper, a new method is put forward, that is a method of processing seismic data in frequency division from wavelet transform and reconstruction. In ordinary seismic processing methods for resolution improvement, deconvolution operator has poor part characteristics, thus influencing the operator frequency. In wavelet transform, wavelet function has very good part characteristics. Frequency - division data processing in wavelet transform also brings quite good high resolution data, but it needs more time than deconvolution method does. On the basis of frequency - division processing method in wavelet domain, a new technique is put forward, which involves 1) designing filter operators equivalent to deconvolution operator in time and frequency domains in wavelet transform, 2) obtaining derivative wavelet function that is suitable to high - resolution seismic data processing, and 3) processing high resolution seismic data by deconvolution method in time domain. In the method of producing some instantaneous characteristic signals by using Hilbert transform, Hilbert transform is very sensitive to high - frequency random noise. As a result, even though there exist weak high - frequency noises in seismic signals, the obtained instantaneous characteristics of seismic signals may be still submerged by the noises. One method for having instantaneous characteristics of seismic signals in wavelet domain is put forward, which obtains directly the instantaneous characteristics of seismic signals by taking the characteristics of both the real part (real signals, namely seismic signals) and the imaginary part (the Hilbert transfom of real signals) of wavelet transform. The method has the functions of frequency division and noise removal. What is more, the weak wave whose frequency is lower than that of high - frequency random noise is retained in the obtained instantaneous characteristics of seismic signals, and the weak wave may be seen in instantaneous characteristic sections (such as instantaneous frequency, instantaneous phase and instantaneous amplitude). Impedance inversion is one of tools in the description of oil reservoir. one of methods in impedance inversion is Generalized Linear Inversion. This method has higher precision of inversion. But, this method is sensitive to noise of seismic data, so that error results are got. The description of oil reservoir in researching important geological layer, in order to give prominence to geological characteristics of the important layer, not only high frequency impedance to research thin sand layer, but other frequency impedance are needed. It is difficult for some impedance inversion method to realize the goal. Wavelet transform is very good in denoising and processing in frequency division. Therefore, in the paper, a method of impedance inversion is put forward based on wavelet transform, that is impedance inversion in frequency division from wavelet transform and reconstruction. in this paper, based on wavelet transform, methods of time - frequency analysis is given. Fanally, methods above are in application on real oil field - Sansan oil field.

工程地质三维建模与可视化技术应用开发

Synthetic Geology Information System (SGIS) is an important constituent part of the theory of Engineering Geomechanics Mate-Synthetic (EGMS), and is the information system more suited for the collection, storage, management, analysis and processing to the information coming from engineering geology,' geological engineering and geotechnical engineering. Its contents involve various works and methods of the investigation, design, and construction in different stages of the geological engineering. Engineering geological and three-dimensional modeling and visualization is the fundamental part of the SGIS, and is a theory, method and technique by which, adopting the computer graphics and image processing techniques, the data derived from engineering geological survey and the calculated results obtained from the geomechanical numerical simulation and analysis are converted to the graphics and images displayed on the computer screen and can be processed interactively. In this paper, the significance and realizing approaches of the three-dimensional modeling and visualization for the complex geological mass in the engineering geology are discussed and the methods of taking advantage of the interpolation and fitting for the scattered and field-surveyed data to simulate the geological layers, such as the topography and earth surface, the groundwater table and the stratum boundary, are researched into. At the mean time, in mind the characteristics of the structure of the basic data for three-dimensional modeling, its visual management can be resolved into the engineering surveyed database management module, plot parameter management module and data output module and the requirement for basic data management can be fulfilled. In the paper, the establishment and development of the three-dimensional geological information system are probed tentatively, and an instance of three-dimensional visual Engineering Distribution Information System (EDIS), theConstruction Management Information System for an airport, in which the functions, such as the real-time browse among the three-dimensional virtual-reality landscapes of the airport construction from start to finish, the information query to the airport facility and the building in the housing district and the recording and playback of the animation sets for the browse and the takeoff and landing of the planes, is developed by applying the component-mode three-dimensional virtual-reality geological information system (GIS) software development kits (SDK), so the three-dimensional visual management platform is provided for the airport construction. Moreover, in the gaper, integrated with the three-dimensional topography visualization and its application in the Sichuan-Tibet Highways, the method of the digital elevation model (DEM) data collection from the topographic maps is described, and the three-dimensional visualization and the roaming about the terrain along the highway are achieved through computer language programming. Understanding to the important role played by the varied and unique topographical condition in the gestation and germination of the highly-dense, frequently-arising and severely-endangered geological hazards can be deepened.

超参数控制下的波阻抗约束反演

The seismic survey is the most effective geophysical method during exploration and development of oil/gas. As a main means in processing and interpreting seismic data, impedance inversion takes up a special position in seismic survey. This is because the impedance parameter is a ligament which connects seismic data with well-logging and geological information, while it is also essential in predicting reservoir properties and sand-body. In fact, the result of traditional impedance inversion is not ideal. This is because the mathematical inverse problem of impedance is poor-pose so that the inverse result has instability and multi-result, so it is necessary to introduce regularization. Most simple regularizations are presented in existent literature, there is a premise that the image(or model) is globally smooth. In fact, as an actual geological model, it not only has made of smooth region but also be separated by the obvious edge, the edge is very important attribute of geological model. It's difficult to preserve these characteristics of the model and to avoid an edge too smooth to clear. Thereby, in this paper, we propose a impedance inverse method controlled by hyperparameters with edge-preserving regularization, the inverse convergence speed and result would be improved. In order to preserve the edge, the potential function of regularization should satisfy nine conditions such as basic assumptions edge preservation and convergence assumptions etc. Eventually, a model with clear background and edge-abnormity can be acquired. The several potential functions and the corresponding weight functions are presented in this paper. The potential functionφLφHL andφGM can meet the need of inverse precision by calculating the models. For the local constant planar and quadric models, we respectively present the neighborhood system of Markov random field corresponding to the regularization term. We linearity nonlinear regularization by using half-quadratic regularization, it not only preserve the edge, and but also simplify the inversion, and can use some linear methods. We introduced two regularization parameters (or hyperparameters) λ2 and δ in the regularization term. λ2 is used to balance the influence between the data term and the transcendental term; δ is a calibrating parameter used to adjust the gradient value at the discontinuous position(or formation interface). Meanwhile, in the inverse procedure, it is important to select the initial value of hyperparameters and to change hyperparameters, these will then have influence on convergence speed and inverse effect. In this paper, we roughly give the initial value of hyperparameters by using a trend- curve of φ-(λ2, δ) and by a method of calculating the upper limit value of hyperparameters. At one time, we change hyperparameters by using a certain coefficient or Maximum Likelihood method, this can be simultaneously fulfilled with the inverse procedure. Actually, we used the Fast Simulated Annealing algorithm in the inverse procedure. This method overcame restrictions from the local extremum without depending on the initial value, and got a global optimal result. Meanwhile, we expound in detail the convergence condition of FSA, the metropolis receiving probability form Metropolis-Hasting, the thermal procession based on the Gibbs sample and other methods integrated with FSA. These content can help us to understand and improve FSA. Through calculating in the theoretic model and applying it to the field data, it is proved that the impedance inverse method in this paper has the advantage of high precision practicability and obvious effect.

用宽频带远震体波波形研究中国及邻近地区地壳构造及地质意义初探

The theory and approach of the broadband teleseismic body waveform inversion are expatiated in this paper, and the defining the crust structure's methods are developed. Based on the teleseismic P-wave data, the theoretic image of the P-wave radical component is calculated via the convolution of the teleseismic P-wave vertical component and the transform function, and thereby a P-wavefrom inversion method is built. The applied results show the approach effective, stable and its resolution high. The exact and reliable teleseismic P waveforms recorded by CDSN and IRIS and its geodynamics are utilized to obtain China and its vicinage lithospheric transfer functions, this region ithospheric structure is inverted through the inversion of reliable transfer functions, the new knowledge about the deep structure of China and its vicinage is obtained, and the reliable seismological evidence is provided to reveal the geodynamic evolution processes and set up the continental collisional theory. The major studies are as follows: Two important methods to study crustal and upper mantle structure -- body wave travel-time inversion and waveform modeling are reviewed systematically. Based on ray theory, travel-time inversion is characterized by simplicity, crustal and upper mantle velocity model can be obtained by using 1-D travel-time inversion preliminary, which introduces the reference model for studying focal location, focal mechanism, and fine structure of crustal and upper mantle. The large-scale lateral inhomogeneity of crustal and upper mantle can be obtained by three-dimensional t ravel-time seismic tomography. Based on elastic dynamics, through the fitting between theoretical seismogram and observed seismogram, waveform modeling can interpret the detail waveform and further uncover one-dimensional fine structure and lateral variation of crustal and upper mantle, especially the media characteristics of singular zones of ray. Whatever travel-time inversion and waveform modeling is supposed under certain approximate conditions, with respective advantages and disadvantages, and provide convincing structure information for elucidating physical and chemical features and geodynamic processes of crustal and upper mantle. Because the direct wave, surface wave, and refraction wave have lower resolution in investigating seismic velocity transitional zone, which is inadequate to study seismic discontinuities. On the contrary, both the converse and reflected wave, which sample the discontinuities directly, must be carefully picked up from seismogram to constrain the velocity transitional zones. Not only can the converse wave and reflected wave study the crustal structure, but also investigate the upper mantle discontinuities. There are a number of global and regional seismic discontinuities in the crustal and upper mantle, which plays a significant role in understanding physical and chemical properties and geodynamic processes of crustal and upper mantle. The broadband teleseismic P waveform inversion is studied particularly. The teleseismic P waveforms contain a lot of information related to source time function, near-source structure, propagation effect through the mantle, receiver structure, and instrument response, receiver function is isolated form teleseismic P waveform through the vector rotation of horizontal components into ray direction and the deconvolution of vertical component from the radial and tangential components of ground motion, the resulting time series is dominated by local receiver structure effect, and is hardly irrelevant to source and deep mantle effects. Receiver function is horizontal response, which eliminate multiple P wave reflection and retain direct wave and P-S converted waves, and is sensitive to the vertical variation of S wave velocity. Velocity structure beneath a seismic station has different response to radial and vertical component of an accident teleseismic P wave. To avoid the limits caused by a simplified assumption on the vertical response, the receiver function method is mended. In the frequency domain, the transfer function is showed by the ratio of radical response and vertical response of the media to P wave. In the time domain, the radial synthetic waveform can be obtained by the convolution of the transfer function with the vertical wave. In order to overcome the numerical instability, generalized reflection and transmission coefficient matrix method is applied to calculate the synthetic waveform so that all multi-reflection and phase conversion response can be included. A new inversion method, VFSA-LM method, is used in this study, which successfully combines very fast simulated annealing method (VFSA) with damped least square inversion method (LM). Synthetic waveform inversion test confirms its effectiveness and efficiency. Broadband teleseismic P waveform inversion is applied in lithospheric velocity study of China and its vicinage. According to the data of high quality CDSN and IRIS, we obtained an outline map showing the distribution of Asian continental crustal thickness. Based on these results gained, the features of distribution of the crustal thickness and outline of crustal structure under the Asian continent have been analyzed and studied. Finally, this paper advances the principal characteristics of the Asian continental crust. There exist four vast areas of relatively minor variations in the crustal thickness, namely, northern, eastern southern and central areas of Asian crust. As a byproduct, the earthquake location is discussed, Which is a basic issue in seismology. Because of the strong trade-off between the assumed initial time and focal depth and the nonlinear of the inversion problems, this issue is not settled at all. Aimed at the problem, a new earthquake location method named SAMS method is presented, In which, the objective function is the absolute value of the remnants of travel times together with the arrival times and use the Fast Simulated Annealing method is used to inverse. Applied in the Chi-Chi event relocation of Taiwan occurred on Sep 21, 2000, the results show that the SAMS method not only can reduce the effects of the trade-off between the initial time and focal depth, but can get better stability and resolving power. At the end of the paper, the inverse Q filtering method for compensating attenuation and frequency dispersion used in the seismic section of depth domain is discussed. According to the forward and inverse results of synthesized seismic records, our Q filtrating operator of the depth domain is consistent with the seismic laws in the absorbing media, which not only considers the effect of the media absorbing of the waves, but also fits the deformation laws, namely the frequency dispersion of the body wave. Two post stacked profiles about 60KM, a neritic area of China processed, the result shows that after the forward Q filtering of the depth domain, the wide of the wavelet of the middle and deep layers is compressed, the resolution and signal noise ratio are enhanced, and the primary sharp and energy distribution of the profile are retained.