In situ electrochemical scanning tunnelling microscopy investigation of structure for horseradish peroxidase and its electrocatalytic property

Immobilization of protein molecules is a fundamental problem for scanning tunnelling microscopy (STM) measurements with high resolution. In this paper, an electrochemical method has been proved to be an effective way to fix native horseradish peroxidase (HRP) as well as inactivated HRP from electrolyte onto a highly oriented pyrolytic graphite (HOPG) surface. This preparation is suitable for both ex situ and in situ electrochemical STM (ECSTM) measurements. In situ STM has been successfully employed to observe totally different structures of HRP in three typical cases: (1) in situ ECSTM reveals an oval-shaped pattern for a single molecule in neutral buffer solution, which is in good agreement with the dimension determined as 6.2 x 4.3 x 1.2. nm(3) by ex situ STM for native HRP; (2) in situ ECSTM shows that the adsorbed HRP molecules on HOPG in a denatured environment exhibit swelling globes at the beginning and then change into a V-shaped pattern after 30 min; (3) in situ ECSTM reveals a black hole in every ellipsoidal sphere for inactivated HRP in strong alkali solution. The cyclic voltammetry results indicate that the adsorbed native HRP can directly catalyse the reduction of hydrogen peroxide, demonstrating that a direct electron transfer reduction occurred between the enzyme and HOPG electrode, whereas the corresponding cyclic voltammograms for denatured HRP and inactivated HRP adsorbed on HOPG electrodes indicate a lack of ability to catalyse H2O2 reduction, which confirms that the HRP molecules lost their biological activity. Obviously, electrochemical results powerfully support in situ STM observations.

A STUDY OF CATALYTIC PROPERTIES OF PEROVSKITE-TYPE OXIDES LAMNYCO1-YO3 .1. THE RELATIONSHIP BETWEEN THE TYPE OF OXYGEN SPECIES AND THE CATALYTIC PROPERTY IN AMMONIA OXIDATION

The type of oxygen species in perovskite-type oxides LaMnyCo1-yO3 (y = 0.0, 0.1, 0.3, 0.5, 0.7, 0.9, 1.0) has been studied by means of XRD, XPS and TPD. The catalytic activity in ammonia oxidation was also investigated. It was found that there were three desorption peaks in TPD curve corresponding to three types of oxygen species (alpha, beta, beta'). The desorption temperatures were 293 K less-than-or-equal-to T(alpha) less-than-or-equal-to 773 K, 773 K less-than-or-equal-to T(beta) less-than-or-equal-to K and T(beta') greater-than-or-equal-to 1073 K respectively. The relationship among the composition, structure and the catalytic property of.the catalyst was correlated and could be explainned with a model based on solid defect reaction and the interaction between Co and Mn ions. The adsorption strength and quantity of a oxygen are proportional to the catalytic activity. The, result indicates that the synergetic effect between B-site ions seems to the benefit of the ammonis oxidation reaction.

STATISTICAL THERMODYNAMICS OF POLYDISPERSE POLYMER BLENDS WITH STRONG INTERACTION

A statistical thermodynamics theory of polydisperse polymer mixtures with strong interaction between dissimilar components based on a lattice fluid model is formulated. Expressions for the free energy, equation of state, phase stability and spinodal for a polydisperse, binary polymer mixture with strong interaction are derived.

STUDY ON STRUCTURE-PROPERTY RELATIONSHIP OF POLYIMIDE BLENDS .1. A DYNAMIC MECHANICAL STUDY OF THE MISCIBILITY OF POLYIMIDE POLYIMIDE BLENDS

Three pairs of polyimide/polyimide blends (50/50 wt%) with different molecular structures were prepared by two ways, i.e. mixing of the polyamic acid precursors with subsequent imidization, and direct solution mixing of the polyimides. The blends were studied with DMA technique. The results obtained show that all the blends prepared with these two different ways are miscible, as there existed only one glass transition temperature(Tg) for all the blends. It is suggested that the miscibility of these polyimide/polyimide blends is a result of the strong inter-molecular charge-transfer interaction between the chains of their components.

Evaluation of antioxidant property of extract and fractions obtained from a red alga, Polysiphonia urceolata

Antioxidant activity (AA), total phenolic content, and reducing power of the crude extract, fractions, and subfractions derived from a red alga, Polysiphonia urceolata, were evaluated and determined. The antioxidative activity was measured using the alpha,alpha-diphenyl-beta-pierylhydrazyl (DPPH) radical scavenging assay and the P-carotene-linoleate assay systems, and compared with that of butylated hydroxytoluene (BHT), gallic acid (GA), and ascorbic acid (AscA). The results showed that the crude extract and the ethyl acetate-soluble fraction exhibited higher AA than BHT in the DPPH assay model, at all of four concentration levels tested (from 0.4 to 50 mu g/ml), while, in the beta-carotene-linoleate assay system, the crude extract and the ethyl acetate-soluble fraction exhibited similar or, in most cases, higher AA than GA and AscA at the same concentrations (from 10 to 200 mu g/ml). The ethyl acetate-soluble fraction was further fractionated into seven subfractions F1-F7 by silica gel vacuum liquid chromatography. F1 was found to be the most effective subfraction in both assay systems. The total phenolic content and reducing power were determined using the Folin-Ciocalteu and the potassium ferricyanide reduction methods, respectively. Statistical analysis indicated a significant association between the antioxidant potency and total phenolic content as well as between the antioxidant potency and reducing power. (c) 2005 Elsevier Ltd. All rights reserved.

Origin, development and evolutionary model of shelf desertization environment in late stage of Upper Pleistocene

Based on the study of palaeo-environmental evolution in the shelves of the Eastern China Seas, the concept of ''shelf desertization'' in the late stage of Upper Pleistocene is defined; the environmental background and evolutionary process of shelf desertization are analysed. Study on the records of subbottom profiling and the data of core samples from shelf areas revealed that during low sea-level stages, the sedimentary environment in the exposed shelf plains was dominated by aeolian depositional process under cold and dry climatic conditions, i.e. under the action of strong winter-monsoon winds. Parts of the exposed marine strata were disintegrated, and aeolian sand dunes were formed on the disintegrated marine deposits, from which the finer sediment grains were blown away by wind and deposited in the downwind areas to form the derivative loess deposits. Thus a desertization environmental system was formed in the exposed shelf plains of the Eastern China Seas.

Impact of difference accuracy on computational properties of vertical grids for a nonhydrostatic model

Starting from nonhydrostatic Boussinesq approximation equations, a general method is introduced to deduce the dispersion relationships. A comparative investigation is performed on inertia-gravity wave with horizontal lengths of 100, 10 and 1 km. These are examined using the second-order central difference scheme and the fourth-order compact difference scheme on vertical grids that are currently available from the perspectives of frequency, horizontal and vertical component of group velocity. These findings are compared to analytical solutions. The obtained results suggest that whether for the second-order central difference scheme or for the fourth-order compact difference scheme, Charny-Phillips and Lorenz ( L) grids are suitable for studying waves at the above-mentioned horizontal scales; the Lorenz time-staggered and Charny-Phillips time staggered (CPTS) grids are applicable only to the horizontal scales of less than 10 km, and N grid ( unstaggered grid) is unsuitable for simulating waves at any horizontal scale. Furthermore, by using fourth-order compact difference scheme with higher difference precision, the errors of frequency and group velocity in horizontal and vertical directions produced on all vertical grids in describing the waves with horizontal lengths of 1, 10 and 100 km cannot inevitably be decreased. So in developing a numerical model, the higher-order finite difference scheme, like fourth-order compact difference scheme, should be avoided as much as possible, typically on L and CPTS grids, since it will not only take many efforts to design program but also make the calculated group velocity in horizontal and vertical directions even worse in accuracy.

A microwave emissivity model of sea surface under wave breaking

With the effective medium approximation theory of composites, a remedial model is proposed for estimating the microwave emissivity of sea surface under wave breaking driven by strong wind on the basis of an empirical model given by Pandey and Kakar. In our model, the effects of the shapes of seawater droplets and the thickness of whitecap layer (i.e. a composite layer of air and sea water droplets) over the sea surface on the microwave emissivity are investigated by calculating the effective dielectric constant of whitecaps layer. The wind speed is included in our model, and the responses of water droplets shapes, such as sphere and ellipsoid, to the emissivity are also discussed at different microwave frequencies. The model is in good agreement with the experimental data of microwave emissivity of sea surface at microwave frequencies of 6.6, 10.7 and 37GHz.

The role of Qiongzhou Strait in the seasonal variation of the South China Sea circulation

An analysis of the water level and current data taken in Qiongzhou Strait in the South China Sea (SCS) over the last 37 years (1963 to 1999) was made to examine the characteristics of tidal waves and residual flow through the strait and their roles in the seasonal variation of the SCS circulation. The observations reveal that Qiongzhou Strait is an area where opposing tidal waves interact and a source of water transport to the Gulf of Beibu (Gulf of Tonkin), SCS. A year-round westward mean flow with a maximum speed of 10-40 cm s(-1) is found in Qiongzhou Strait. This accounts for water transport of 0.2-0.4 Sv and 0.1-0.2 Sv into the Gulf of Beibu in winter-spring and summer-autumn, respectively. The outflow from Qiongzhou Strait may cause up to 44% of the gulf water to be refreshed each season, suggesting that it has a significant impact on the seasonal circulation in the Gulf of Beibu. This finding is in contrast to our current understanding that the seasonal circulation patterns in the South China Sea are primarily driven by seasonal winds. Several numerical experiments were conducted to examine the physical mechanisms responsible for the formation of the westward mean flow in Qiongzhou Strait. The model provides a reasonable simulation of semidiurnal and diurnal tidal waves in the strait and the predicted residual flow generally agrees with the observed mean flow. An analysis of the momentum equations indicates that the strong westward flow is driven mainly by tidal rectification over variable bottom topography. Both observations and modeling suggest that the coastal physical processes associated with tidal rectification and buoyancy input must be taken into account when the mass balance of the SCS circulation is investigated, especially for the regional circulation in the Gulf of Beibu.

Fourier有限差分法深度偏移与正演模拟

With the development of seismic exploration, the target becomes more and more complex, which leads to a higher demand for the accuracy and efficiency in 3D exploration. Fourier finite-difference (FFD) method is one of the most valuable methods in complex structure exploration, which keeps the ability of finite-differenc method in dealing with laterally varing media and inherits the predominance of the phase-screen method in stablility and efficiency. In this thesis, the accuracy of the FFD operator is highly improved by using simulated annealing algorithm. This method takes the extrapolation step and band width into account, which is more suitable to various band width and discrete scale than the commonely-used optimized method based on velocity contrast alone. In this thesis, the FFD method is extended to viscoacoustic modeling. Based on one-way wave equation, the presented method is implemented in frequency domain; thus, it is more efficient than two-way methods, and is more convenient than time domain methods in handling attenuation and dispersion effects. The proposed method can handle large velocity contrast and has a high efficiency, which is helpful to further research on earth absorption and seismic resolution. Starting from the frequency dispersion of the acoustic VTI wave equation, this thesis extends the FFD migration method to the acoustic VTI media. Compared with the convetional FFD method, the presented method has a similar computational efficiency, and keeps the abilities of dealing with large velocity contrasts and steep dips. The numerical experiments based on the SEG salt model show that the presented method is a practical migration method for complex acoustical VTI media, because it can handle both large velocity contrasts and large anisotropy variations, and its accuracy is relatively high even in strong anisotropic media. In 3D case, the two-way splitting technique of FFD operator causes artificial azimuthal anisotropy. These artifacts become apparent with increasing dip angles and velocity contrasts, which prevent the application of the FFD method in 3D complex media. The current methods proposed to reduce the azimuthal anisotropy significantly increase the computational cost. In this thesis, the alternating-direction-implicit plus interpolation scheme is incorporated into the 3D FFD method to reduce the azimuthal anisotropy. By subtly utilizing the Fourier based scheme of the FFD method, the improved fast algorithm takes approximately no extra computation time. The resulting operator keeps both the accuracy and the efficiency of the FFD method, which is helpful to the inhancements of both the accuracy and the efficiency for prestack depth migration. The general comparison is presented between the FFD operator and the generalized-screen operator, which is valuable to choose the suitable method in practice. The percentage relative error curves and migration impulse responses show that the generalized-screen operator is much sensiutive to the velocity contrasts than the FFD operator. The FFD operator can handle various velocity contrasts, while the generalized-screen operator can only handle some range of the velocity contrasts. Both in large and weak velocity contrasts, the higher order term of the generalized-screen operator has little effect on improving accuracy. The FFD operator is more suitable to large velocity contrasts, while the generalized-screen operator is more suitable to middle velocity contrasts. Both the one-way implicit finite-difference migration and the two-way explicit finite-differenc modeling have been implemented, and then they are compared with the corresponding FFD methods respectively. This work gives a reference to the choosen of proper method. The FFD migration is illustrated to be more attractive in accuracy, efficiency and frequency dispertion than the widely-used implicit finite-difference migration. The FFD modeling can handle relatively coarse grids than the commonly-used explicit finite-differenc modeling, thus it is much faster in 3D modeling, especially for large-scale complex media.

地下水位上升造成的黄土远程滑坡形成机理与数值模拟研究

Based on field survey, laboratory testing and numerical modeling, engineering characteristics of undisturbed loess and the mechanism of long-runout loess landslides caused by underground water level rise, as well as the formation conditions and spatial distribution of landslides, are systematically studied and analyzed. Loess landslides at south Plateau of Jingyang County are mainly classified as flowslide, slide and fall. Flowslide is the main type characteristic of high velocity, long runout and multi-stages. The steep relief composed of loose structured loess-old aged soil serials and the rise of groundwater table are the predominant conditions for landslides in the study area. To study loess mechanic poperties and loess landslides mechanisims, isotropically and anisotropically consolidated undrained compression(ICU and ACU) tests and constant-deviator-drained compression (CQD) tests were carried out on undisturbed samples. The results of undrained compression tests performed at the in-situ stress level show that the soils are of consistently strain-softening in the stress-strain relations and cause high excess pore pressure. The steady-state line and the potential region of instability are obtained from ICU and ACU test results. A necessary condition for liquefaction is that the soil state initially lies in or is brought into the potential instability region. In addition, a strong strain-softening model is also formed. CQD tests demonstrate that the mobilized friction angle is far less than the steady-state angle and that the soil experiences undrained contractive failure suddenly at very small strains when its stress path during drained loading tries to cross the potential instability region，thus validates the proposed instability region. Based on the location of the region of potential instability and the stress state of slope soil, a method of static liquefaction analysis is proposed for loess landslides caused by rise in groundwater table. Compared with other liquefaction analysis methods, this method overcomes the limitations inherent in conventional slope stability method and undrained brittleness index method. Triaxial tests composed of constant water content (CW) and wetting tests at constant deviator stress are performed on undisturbed unsaturated samples. The stress-strain relation of CW tests takes on strain-hardening behavior; The results of wetting tests at constant deviator stress designed to study the mechanics of failure of unsaturated loess caused by an increase in the degree of saturation (wetting) shows that a contractive failure occurs in the undisturbed samples. On the basis of the above triaxial test results, the initiation of static liquefaction is presented for long-runout loess landslides caused by rise in groundwater table, that is, the loess slope soil gradually transfer from unsaturated to saturated state under the infiltration of irrigation. A contractive failure occurs in the local region at very small strain by increasing the pore-water pressure at constant deviator stresses under drained conditons. It is the contractive failrue resulting from rise of pore pressure that leads to high excess pore pressure in the neighbour soil which reduces shear resistance of soil. The neighbour soils also fail due to the rapid increase in pore-water pressure. Thus a connected failure surface is developed quickly and a flowslide occurs. Based on the saturated-unsaturated seepage theory, transient seepage is computed using the finite element method on loess slope under groundwater table rise. Pore-water pressure distribution for every time step after irrigation are obtained. The phreatic surface in the slope increases with the groundwater table. Pore-water pressure distribution within 8m above the phreatic surface changes very quickly，but the water content and pore water pressure in the region ranging from 8m above the phreatic surface up to ground surface is almost not affected and the matric suction usually is kept at 100～120 kPa. Based on the results of laboratory tests and seepage flow analysis, the development process of loess landslide is modeled considering groundwater table rise. The shearing plastic zone first occurs at the slope toe where the soil is soaked for long term during rise in groundwater table. As irrigation continues, the shearing plastic zone gradually extends to the interior soils, with the results that the tensile plastic zone occurs at the slope crown. As time goes on, both the shearing plastic zone and tensile plastic zone continue to extend. Then a connected plastic zone is formed and fowslide occurs. In comparision to laboratory test results, the results of numerical simulation quite well verify the presented mechanism of static liquefaction of long-runout loess landslides caused by rise in groundwater table.

细粒沉积物多组分粒度特征及其成因机制研究

By using high-resolution laser grain size instrument Mastersizer 2000, the grain size distribution of windblown depositions (loess and sandy dunes), aqueous sediments (lake, river, riverside and foreshore sand), weathering crust, sloping materials and other fine-grain sediments are systemically measured. The multimodal characteristics of grain size distribution of these sediments are carefully studied. The standard patterns and their grain size characteristics of various sediments are systemically summarized. The discrepancies of multimodal distribution among windblown depositions, aqueous sediments and other sediments are concluded and the physical mechanisms of grain size multimodal distribution of various sediments are also discussed in this paper. The major conclusions are followed: 1. The multimodal characteristic of grain size distribution is a common feature in all sediments and results from properties of transportation medium, dynamic intensity, transportation manner and other factors. 2. The windblown depositions are controlled by aerodynamic forcing, resulting in that the median size of the predominant mode gradually decreases form sandy dunes to loess. Similarly, the aqueous sediments are impacted by dynamic forces of water currents and the median grain size of the predominant mode decreases gradually from river to lake sediments. Because the kinetic viscidity of air is lower than of water, the grain size of modes of windblown depositions is usually finer than that of corresponding modes of aqueous sediments. Typical characteristics of sediments grain size distribution of various sediments have been summarized in the paper: (1) Suspended particles which diameters are less than 75μm are dominant in loess and dust. There are three modes in loess’ grain size distribution: fine, median and coarse (the median size is <1μm、1-10μm、10-75μm, respectively). The coarse mode which percentage is larger than that of others is controlled by source distance and aerodynamic intensity of dust source areas. Some samples also have a saltation mode which median size is about 300-500μm. Our analysis demonstrates that the interaction of wind, atmospheric turbulence, and dust grain gravity along the dust transportation path results in a multimodal grain size distribution for suspended dust. Changes in the median sizes of the coarse and medium modes are related to variation in aerodynamic forcing (lift force related to vertical wind and turbulence) during dust entrainment in the source area and turbulence intensity in the depositional area. (2) There is a predominant coarse saltation mode in grain size distribution of sandy dunes, which median size is about 100-300μm and the content is larger than that of other modes. The grain size distribution curve is near axis symmetric as a standard logarithm normal function. There are some suspended particles in some samples of sandy dunes, which distribution of the fine part is similar to that of loess. Comparing with sandy samples of river sediments, the sorting property of sandy dunes is better than of river samples although both they are the saltation mode. Thus, the sorting property is a criterion to distinguish dune sands and river sands. (3) There are 5~6 modes (median size are <1μm, 1-10μm, 10-70μm, 70-150μm, 150-400μm, >400μm respectively) in grain size distribution of lacustrine sediments. The former 4 modes are suspensive and others are saltated. Lacustrine sediments can be divided into three types: lake shore facies, transitional facies and central lake facies. The grain size distributions of the three facies are distinctly different and, at the same time, the transition among three modes is also clear. In all these modes, the third mode is a criteria to identify the windblown deposition in the watershed. In lake shore sediments, suspended particles are dominant, a saltation mode sometimes occurs and the fourth mode is the most important mode. In the transitional facies, the percentage of the fourth mode decreases and that of the second mode increases from lake shore to central lake. In the central lake facies, the second mode is dominant. A higher content of the second mode indicates its position more close to the central lake. (4) The grain size distribution of river sediments is the most complex. It consist of suspension, saltation and rolling modes. In most situations, the percentage of the saltation mode is larger than that of other modes. The percentage of suspension modes of river sediments is more than of sandy dunes. The grain size distribution of river sediments indicates dynamic strength of river currents. If the fourth mode is dominant, the dynamic forcing of river is weaker, such as in river floodplain. If the five or sixth mode is dominant, the water dynamic forcing of rivers is strong. (5) Sediments can be changed by later forcing in different degree to form some complicated deposition types. In the paper, the grain size distribution of aqueous sediments of windblown deposition, windblown sediments of aqueous deposition, weathering crust and slope materials are discussed and analyzed. 3. The grain size distribution characteristics of different sediments are concluded: (1) Modal difference: Usually there are suspended and saltation modes in the windblown deposition. The third mode is dominant in loess dust and the fifth mode is predominant in sandy dunes. There are suspended, saltation and rolling particles in aqueous sediments. In lacustrine sediments, the second and fourth mode are predominant for central lake facies or lake shore facies, respectively. In river sediments, the fourth, or fifth, or sixth mode is predominant. Suspended modes: the grain size of suspended particles of windblown depositions usually is less than 75μm. The content of suspended particles is lower or none in sandy dunes. However, suspended particles of aqueous sediments may reach 150μm. Difference in grain size of suspended modes represents difference between transitional mediums and the strength of dynamic forcing. Saltation modes: the median size of saltation mode of sandy dunes fluctuates less than that of river sediments. (2) Loess dust and lacustrine sediment: Their suspended particles are clearly different. There is an obvious pit between the second and the third modes in grain size distribution of lacustrine sediments. The phenomenon doesn’t occur in loess dust. In lacustrine sediments, the second mode can be a dominant mode, such as central lake facies, and contents of the second and the third modes change reversely. However, the percentage of the third mode is always the highest in loess dust. (3) Dune Sand and fluvial sand: In these two depositions, the saltation particles are dominant and the median sizes of their saltation modes overlay in distribution range. The fifth mode of dune sand fluctuates is sorted better than that of fluvial sand. (4) Lacustrine and fluvial sediments: In lacustrine sediments, there are 5-6 modes and suspended particles can be predominant. The second mode is dominant in central lake facies and the third mode is dominant in lake shore facies. Saltation or roll modes occurred in central lake facies may indicate strong precipitation events. In fluvial sediments, saltation particles (or rolling particles) usually dominant. 4. A estimation model of lake depth is firstly established by using contents of the second, the third and the fourth modes. 5. The paleoenvironmental history of the eastern part of SongLiao basin is also discussed by analyzing the grain size distribution of Yushu loess-like sediments in Jilin. It was found that there is a tectonic movement before 40ka B.P. in SongLiao basin. After the movement, loess dust deposited in Yushu area as keerqin desert developed. In recent 2000 years, the climate became drier and more deserts activated in the eastern part of Song-Liao basin.

基于量子蒙特卡罗的地震波反演方法

The primary approaches for people to understand the inner properties of the earth and the distribution of the mineral resources are mainly coming from surface geology survey and geophysical/geochemical data inversion and interpretation. The purpose of seismic inversion is to extract information of the subsurface stratum geometrical structures and the distribution of material properties from seismic wave which is used for resource prospecting, exploitation and the study for inner structure of the earth and its dynamic process. Although the study of seismic parameter inversion has achieved a lot since 1950s, some problems are still persisting when applying in real data due to their nonlinearity and ill-posedness. Most inversion methods we use to invert geophysical parameters are based on iterative inversion which depends largely on the initial model and constraint conditions. It would be difficult to obtain a believable result when taking into consideration different factors such as environmental and equipment noise that exist in seismic wave excitation, propagation and acquisition. The seismic inversion based on real data is a typical nonlinear problem, which means most of their objective functions are multi-minimum. It makes them formidable to be solved using commonly used methods such as general-linearization and quasi-linearization inversion because of local convergence. Global nonlinear search methods which do not rely heavily on the initial model seem more promising, but the amount of computation required for real data process is unacceptable. In order to solve those problems mentioned above, this paper addresses a kind of global nonlinear inversion method which brings Quantum Monte Carlo (QMC) method into geophysical inverse problems. QMC has been used as an effective numerical method to study quantum many-body system which is often governed by Schrödinger equation. This method can be categorized into zero temperature method and finite temperature method. This paper is subdivided into four parts. In the first one, we briefly review the theory of QMC method and find out the connections with geophysical nonlinear inversion, and then give the flow chart of the algorithm. In the second part, we apply four QMC inverse methods in 1D wave equation impedance inversion and generally compare their results with convergence rate and accuracy. The feasibility, stability, and anti-noise capacity of the algorithms are also discussed within this chapter. Numerical results demonstrate that it is possible to solve geophysical nonlinear inversion and other nonlinear optimization problems by means of QMC method. They are also showing that Green’s function Monte Carlo (GFMC) and diffusion Monte Carlo (DMC) are more applicable than Path Integral Monte Carlo (PIMC) and Variational Monte Carlo (VMC) in real data. The third part provides the parallel version of serial QMC algorithms which are applied in a 2D acoustic velocity inversion and real seismic data processing and further discusses these algorithms’ globality and anti-noise capacity. The inverted results show the robustness of these algorithms which make them feasible to be used in 2D inversion and real data processing. The parallel inversion algorithms in this chapter are also applicable in other optimization. Finally, some useful conclusions are obtained in the last section. The analysis and comparison of the results indicate that it is successful to bring QMC into geophysical inversion. QMC is a kind of nonlinear inversion method which guarantees stability, efficiency and anti-noise. The most appealing property is that it does not rely heavily on the initial model and can be suited to nonlinear and multi-minimum geophysical inverse problems. This method can also be used in other filed regarding nonlinear optimization.

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

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

鄂尔多斯盆地志丹地区延长组储层综合评价

Sedimentary provenance direction，sedimentary facies，reservoir geological characteristic，pore structure； physical property characteristic，reservoir classification and evaluation ，forthermore，favorable area area are forecasted of Yanchang formation in ZhiDan region， by mainly study on the data of field section observation ，core observation, well logging explaination and routine microscope slice identification，scanning Electron Microscope and reservoir analysis of lithology and physical property , Under the guidance of such advanced theories and methods as sedimentology，reservoir sedimentology，lithological oil pool and so on，in the thesis. The following fruits are mainly achieved in this paper： Yanchang formation stratum is divided and correlated in this entire region， and the characteristic of oil layer unit is detailed discussed , respectively. According to main marker bed and supplementary ones.and research result shows that the source of provenance direction during Yanchang Formation mianly is north-east. Delta and lake are mainly developed in study area ，sub-facies and micro-facies are divided，distribution of sedimentary micro-facies in plane and palaeogeographic evolution are described，and gentle slope type- shallow water delta depositional model is established. Fine-grain arkose sandstone is the main reservoir，and which have experienced such different degree diagenesis as compaction, cementation, replacement and dissolution, and in which compaction and cementation are mainly factors to reduce sandstone physical property and dissolution effectively improved physical property during burial diagenesis procedure. All reservoirs of Yanchang Formation have entered A period of late diagenetic stage according to scheme of diagenesis period division . Intergranular porosity，dissolution porosity，fissure porosity are main pore types. And porosity structure are analyse by mercury penetration capillary pressure parameter，fine-shortness type and fine- length throat type are mainly developed. as a whole，the reservoir， with the characteristic of porosity and permeability altering apparently，strong inhomogeneity ， is a medium- porosity and medium permeability one. In plane，higher- porosity and higher-permeability are corresponded well with distributary channel area， physical property and inhomogeneity are affected by both deposition and diagenesis，and distributary channel and underwater distributary channel are favorable facies . According to such characteristic as lithology，physical property，pore structure ，diagenesis and sandstone distribution, the sandy reservoir can be classified 4 types, and the main sandy in every oil layer unit are evaluated according to the standard. The analysis result of petroleum concentration rule shows that Yanchang Formation are with not only favourable oil source rock，reservoir，covering combination ，but also good entrapment condition in study area. Lithology and structure-lithology oil pool are mainly developed ，based on condition of favorable reservoir developments，accounting for deliverability and sandstone superface elevation，zone of profitabilitis are forecasted.