裂缝储层预测技术应用研究

In practice,many fracture reservoir was found,which has giant potential for exporation. For example,in limestone fracture reservoir,igneous rock fracture reservoir and shale fracture reservoir ,there are high yield oil wells found. The fracture reservoir has strong anisotropy and is very difficult to explore and produce.Since 1990’s,the techniques that use structure information and P-WAVE seismic attributes to detect fracture developed very rapidly,include stress and strain analysis,using amplitude,interval velocity,time-difference,azimuthal AVO analysis etc. Based on research and develop these advanced techniques of fracture detect,this paper selected two typical fracture reservoir as target area,according to the characters of research area,selected different techniques to pridect the fracture azimuth and density of target,and at last ,confirmed the favored area. This paper includes six parts:the first chapter mainly addresses the domestic and international research actuality about the fracture prediction and the evolement in ShengLi oil field,then according to the temporal exploration requirement,a research route was established; Based on the close relationship between structural fracture and the geotectonic movement and the procedure of rock distortion，the second chapter research the structural fracture predicting technique which is realized by computing the strain in every geotectonic movement ,which is by use of the forward and inversion of the growing history of structure; The third chapter discussed many kind of traditional techniques for fracture reservoir prediction,and point out their disadvantages.then research and develop the coherence volume computing technique which can distinguish from faults,the seismic wave absorbing technique,and other fracture predicting technique which is by use of seismic attributes ,such as azimuthal AVO FVO etc; The fourth chapter first establish the geological and petrophysical model by use of the existed log and drill well information, then research the variation of amplitude and seismic wave which is caused by fractures.based on it , the fracture predicting technique which is by use of variation of azimuthal impedance is researched;The fifth chapter is a case study,it selects shale fracture reservoir in LuoJia area as target,selects several kind of techniques to apply ,at last ,the fracture distribution of target reservoir and favored area were gotten;the sixth chapter is another case study,it selects limestone fracture reservoir in BoShen6 buried hill as target,selects several kind of techniques to apply,similarly favored area were gotten. Based on deeply research and development of the new techniques for fracture reservoir exploration, This paper selects two fracture reservoirs the most typical in ShengLi as targets to be applied ,good results show up a good application way ,which can be used for reference for future fracture exploration,and it can bring materially economic and social benefit.

膨胀土化学改性实验研究

Expansive soil is a kind of typical unsaturated soil with characteristics of high swelling-shrinking deformation, cracks and over consolidation. It is very harmful to civil engineering, As a new processing method deal with expansive soil, Chemistry treatment has widespread applied in developed countries such as Europe and America, and also gained remarkable result. Based on the embankment filling soil improving testing projects in Meng-Xin freeway, this paper proposed a new processing method to expansive soil embankment wrapped with PAS-treated soil, experimental study of expansive soil chemical improved by PAS is been carried out. The water content change is the external factor which causes expansive soil to have swelling-shrinkage deformation. this reflected that the soil body swelling-shrinkage characteristic mainly depends on its mineral ingredient and the soil-water mutual function. This paper takes expansive soil as one kind of ordinary high plastic clay from angle of clay-water mutual function explained the expansive soil swelling-shrinkage deformation mechanism on microscopic. And take this swelling-shrinkage mechanism as the master line, Cooperates with the China Academy of Chemistry, we developed the new method PAS treatment, trough ionic exchange, joint, package and flocculation, the stronger static electricity function weakened the level through adsorption and the stronger static electricity function, PAS can weakened the negative charge repulsion between levels, causes the electric potential to reduce, diffusion layer thickness to be thinner, and improves the water affinity performance of expansive soil effectively. Moreover the space network architecture compromised with PAS and soil enhanced the joint strength between the clay particles , enable the soil body to have comparatively high strength and the distortion rate. pointed proposed the PAS modified principle. Combine with the construction of experimented road, this paper sums up and presents the construction craft and technology requirement of PAS treatment to expansive soil embankment. Through many experimental studied the basic physical property, the intensity characteristic and water stability changes of expansive soil and PAS-treated soil. The results of study indicate that adding lime into the expansive soil can reduce the content of clay gain obviously, reduce the plasticity notably, increase the strength greatly, control the property of swelling and shrinking effectively, and can meliorate the stability of sucking water clearly. Simultaneity PAS don’t change the cultivate capacity of the soil, the modified slope of the embankment can adopt plant fixed slope method as ecology protection. Finally the processing effect of use different treatment has analyzed through numerical simulation, summarized the PAS chemical wrapping treatment process in the actual project application, and appraised its processing effect and the project efficiency. The research indicated that PAS chemical treatment is one effective method to improve expansive soil. Compare with long-distance replacement, especially in the high plastic expansive soil massive distribution area, PAS treatment has the very greatly economical superiority to be promoted. The study in the paper not only afforded technique method to Meng-Xin expressway construction but also important for improvement of the expressway construction theory in swelling soil areas. Key words: PAS; expansive soil; swelling-shrinkage deformation mechanism; wrapping embankment; chemical modified treatment.

青藏高原东部电性结构特征及其构造意义

Eastern Himalayan Syntaxis (EHS) and its surroundings (eastern margin of Tibet) is one of the most complicated tectonic areas in the world. As the exhaust opening of the balanced materials of the Tibetan Plateau during the collision of Indan and Eurasian plates, the deep structure beneath EHS surrounding region is referred to as the key to the study of the dynamics of the plateau. EHS3D project, sponsored by NSFC, has been proposed to explore the deep electric features of the area. During the first stage of EHS3D(2006-2008), MT+LMT measurements have been conducted along two lines from Chayu to Qingshuihe (EHS3D-3) and Chayu to Ruoergai (EHS3D-2). This paper will discuss the MT models of EHS3D-3 line. By the data procrssing, including distortion analysis, Robust estimation and strike decomposition, rotated apparent resitivities and phases have been obtained for each station. Then conventional 2-D inversion algorithms (NLCG and RRI) were employed to produce 2-D models. The final preferred 2-D model suggests that the upper crust consists of resistive blocks while in mid-lower crust there are two extensive conductive bodies beneath Lhasa block and Qiangtang terrain respectively. Jinshajiang suture is a gradient belt and Bangong-Nujiang suture appear a conductive belt dipping to the north. . We concluded that the formation of the two conductive bodies attributed to the partial melt and fluids in the lower crust. The regional electric strike derived from decomposition analysis indicates that the crust and upper mantle move in different manners. The upper crust moves like slips of rigid blocks along major slip faults while the lower crust creeps as a flow in the conductive channels.

粘弹性介质中的地震波传播与波形反演研究

The real earth is far away from an ideal elastic ball. The movement of structures or fluid and scattering of thin-layer would inevitably affect seismic wave propagation, which is demonstrated mainly as energy nongeometrical attenuation. Today, most of theoretical researches and applications take the assumption that all media studied are fully elastic. Ignoring the viscoelastic property would, in some circumstances, lead to amplitude and phase distortion, which will indirectly affect extraction of traveltime and waveform we use in imaging and inversion. In order to investigate the response of seismic wave propagation and improve the imaging and inversion quality in complex media, we need not only consider into attenuation of the real media but also implement it by means of efficient numerical methods and imaging techniques. As for numerical modeling, most widely used methods, such as finite difference, finite element and pseudospectral algorithms, have difficulty in dealing with problem of simultaneously improving accuracy and efficiency in computation. To partially overcome this difficulty, this paper devises a matrix differentiator method and an optimal convolutional differentiator method based on staggered-grid Fourier pseudospectral differentiation, and a staggered-grid optimal Shannon singular kernel convolutional differentiator by function distribution theory, which then are used to study seismic wave propagation in viscoelastic media. Results through comparisons and accuracy analysis demonstrate that optimal convolutional differentiator methods can solve well the incompatibility between accuracy and efficiency, and are almost twice more accurate than the same-length finite difference. They can efficiently reduce dispersion and provide high-precision waveform data. On the basis of frequency-domain wavefield modeling, we discuss how to directly solve linear equations and point out that when compared to the time-domain methods, frequency-domain methods would be more convenient to handle the multi-source problem and be much easier to incorporate medium attenuation. We also prove the equivalence of the time- and frequency-domain methods by using numerical tests when assumptions with non-relaxation modulus and quality factor are made, and analyze the reason that causes waveform difference. In frequency-domain waveform inversion, experiments have been conducted with transmission, crosshole and reflection data. By using the relation between media scales and characteristic frequencies, we analyze the capacity of the frequency-domain sequential inversion method in anti-noising and dealing with non-uniqueness of nonlinear optimization. In crosshole experiments, we find the main sources of inversion error and figure out how incorrect quality factor would affect inverted results. When dealing with surface reflection data, several frequencies have been chosen with optimal frequency selection strategy, with which we use to carry out sequential and simultaneous inversions to verify how important low frequency data are to the inverted results and the functionality of simultaneous inversion in anti-noising. Finally, I come with some conclusions about the whole work I have done in this dissertation and discuss detailly the existing and would-be problems in it. I also point out the possible directions and theories we should go and deepen, which, to some extent, would provide a helpful reference to researchers who are interested in seismic wave propagation and imaging in complex media.

青藏高原东部电性结构特征及地球动力学意义

As the largest and highest plateau on the Earth, the Tibetan Plateau has been a key location for understanding the processes of mountain building and plateau formation during India-Asia continent-continent collision. As the front-end of the collision, the geological structure of eastern Tibetan Plateau is very complex. It is ideal as a natural laboratory for investigating the formation and evolution of the Tibetan Plateau. Institute of Geophysics, Chinese Academy of Sciences (CAS) carried out MT survey from XiaZayii to Qingshuihe in the east part of the plateau in 1998. After error analysis and distortion analysis, the Non-linear Conjugate Gradient inversion(NLCG), Rapid Relaxation Inversin (RRI) and 2D OCCAM Inversion algorithms were used to invert the data. The three models obtained from 3 algorithms provided similar electrical structure and the NLCG model fit the observed data better than the other two models. According to the analysis of skin depth, the exploration depth of MT in Tibet is much more shallow than in stable continent. For example, the Schmucker depth at period 100s is less than 50km in Tibet, but more than 100km in Canadian Shield. There is a high conductivity layer at the depth of several kilometers beneath middle Qiangtang terrane, and almost 30 kilometers beneath northern Qiangtang terrane. The sensitivity analysis of the data predicates that the depth and resistivity of the crustal high conductivity layer are reliable. The MT results provide a high conductivity layer at 20~40km depth, where the seismic data show a low velocity zone. The experiments show that the rock will dehydrate and partially melt in the relative temperature and pressure. Fluids originated from dehydration and partial melting will seriously change rheological characteristics of rock. Therefore, This layer with low velocity and high conductivity layer in the crust is a weak layer. There is a low velocity path at the depth of 90-110 km beneath southeastern Tibetan Plateau and adjacent areas from seismology results. The analysis on the temperature and rheological property of the lithosphere show that the low velocity path is also weak. GPS measurements and the numerical simulation of the crust-mantle deformation show that the movement rate is different for different terranes. The regional strike derived from decomposition analysis for different frequency band and seismic anisotropy indicate that the crust and upper mantle move separately instead of as a whole. There are material flow in the eastern and southeastern Tibetan Plateau. Therefore, the faults, the crustal and upper mantle weak layers are three different boundaries for relatively movement. Those results support the "two layer wedge plates" geodynamic model on Tibetan formation and evolution.

Fast Contour Matching Using Approximate Earth Mover's Distance

Weighted graph matching is a good way to align a pair of shapes represented by a set of descriptive local features; the set of correspondences produced by the minimum cost of matching features from one shape to the features of the other often reveals how similar the two shapes are. However, due to the complexity of computing the exact minimum cost matching, previous algorithms could only run efficiently when using a limited number of features per shape, and could not scale to perform retrievals from large databases. We present a contour matching algorithm that quickly computes the minimum weight matching between sets of descriptive local features using a recently introduced low-distortion embedding of the Earth Mover's Distance (EMD) into a normed space. Given a novel embedded contour, the nearest neighbors in a database of embedded contours are retrieved in sublinear time via approximate nearest neighbors search. We demonstrate our shape matching method on databases of 10,000 images of human figures and 60,000 images of handwritten digits.

A Constant-Factor Approximation Algorithm for Embedding Unweighted Graphs into Trees

We present a constant-factor approximation algorithm for computing an embedding of the shortest path metric of an unweighted graph into a tree, that minimizes the multiplicative distortion.

Shape from Contour

The problem of using image contours to infer the shapes and orientations of surfaces is treated as a problem of statistical estimation. The basis for solving this problem lies in an understanding of the geometry of contour formation, coupled with simple statistical models of the contour generating process. This approach is first applied to the special case of surfaces known to be planar. The distortion of contour shape imposed by projection is treated as a signal to be estimated, and variations of non-projective origin are treated as noise. The resulting method is then extended to the estimation of curved surfaces, and applied successfully to natural images. Next, the geometric treatment is further extended by relating countour curvature to surface curvature, using cast shadows as a model for contour generation. This geometric relation, combined with a statistical model, provides a measure of goodness-of-fit between a surface and an image contour. The goodness-of-fit measure is applied to the problem of establishing registration between an image and a surface model. Finally, the statistical estimation strategy is experimentally compared to human perception of orientation: human observers' judgements of tilt correspond closely to the estimates produced by the planar strategy.

Low-temperature oxidation of CO over Pd/CeO2-TiO2 catalysts with different pretreatments

A comprehensive study of the low-temperature oxidation of CO was conducted over Pd/TiO2, Pd/CeO2, and Pd/CeO2-TiO2 pretreated by a series of calcination and reduction processes. The catalysts were characterized by N-2 adsorption, XRD, H-2 chemisorption, and diffuse-reflectance infrared Fourier transform spectroscopy. The results indicated that Pd/CeO2-TiO2 has the highest activity among these catalysts, whether in the calcined state or in the reduced state. The activity of all of the catalysts can be improved significantly by the pre-reduction, and it seems that the reduction at low temperature (LTR. 150 degrees C) is more effective than that at high temperature (HTR, 500 degrees C), especially for Pd/CeO2 and Pd/TiO2. The catalysts with various supports and pretreatments are also different in the reaction mechanisms for CO oxidation at low temperature. Over Pd/TiO2, the reaction may proceed through a surface reaction between the weakly adsorbed CO and oxygen (Langmuir-Hinshelwood). For Ce-containing catalysts, however, an alteration of reaction mechanism with temperature and the involvement of the oxygen activation at different sites were observed, and the light-off profiles of the calcined Pd/CeO2 and Pd/CeOi-TiO2 show a distortion before CO conversion achieves 100%. At low temperature, CO oxidation proceeds mainly via the reaction between the adsorbed CO on Pd-0 sites and the lattice oxygen of surface CeO2 at the Pd-Ce interface, whereas at high temperature it proceeds via the reaction between the adsorbed CO and oxygen. The high activity of Pd/CeO2-TiO2 for the low-temperature CO oxidation was probably due to the enhancements of both CO activation, caused by the facilitated reduction of Pd2+ to Pd-0, and oxygen activation, through the improvement of the surface oxygen supply and the oxygen vacancies formation. The reduction pretreatment enhances metal-support interactions and oxygen vacancy formation and hence improves the activity of CO oxidation. (c) 2005 Elsevier Inc. All rights reserved.

PROPAGATION OF SONIC BOOMS THROUGH A REAL, STRATIFIED ATMOSPHERE

Sonic boom propagation in a quiet) stratified) lossy atmosphere is the subject of this dissertation. Two questions are considered in detail: (1) Does waveform freezing occur? (2) Are sonic booms shocks in steady state? Both assumptions have been invoked in the past to predict sonic boom waveforms at the ground. A very general form of the Burgers equation is derived and used as the model for the problem. The derivation begins with the basic conservation equations. The effects of nonlinearity) attenuation and dispersion due to multiple relaxations) viscosity) and heat conduction) geometrical spreading) and stratification of the medium are included. When the absorption and dispersion terms are neglected) an analytical solution is available. The analytical solution is used to answer the first question. Geometrical spreading and stratification of the medium are found to slow down the nonlinear distortion of finite-amplitude waves. In certain cases the distortion reaches an absolute limit) a phenomenon called waveform freezing. Judging by the maturity of the distortion mechanism, sonic booms generated by aircraft at 18 km altitude are not frozen when they reach the ground. On the other hand, judging by the approach of the waveform to its asymptotic shape, N waves generated by aircraft at 18 km altitude are frozen when they reach the ground. To answer the second question we solve the full Burgers equation and for this purpose develop a new computer code, THOR. The code is based on an algorithm by Lee and Hamilton (J. Acoust. Soc. Am. 97, 906-917, 1995) and has the novel feature that all its calculations are done in the time domain, including absorption and dispersion. Results from the code compare very well with analytical solutions. In a NASA exercise to compare sonic boom computer programs, THOR gave results that agree well with those of other participants and ran faster. We show that sonic booms are not steady state waves because they travel through a varying medium, suffer spreading, and fail to approximate step shocks closely enough. Although developed to predict sonic boom propagation, THOR can solve other problems for which the extended Burgers equation is a good propagation model.

Intersubject Regularity in the Intrinsic Shape of Human V1

Previous studies have reported considerable intersubject variability in the three-dimensional geometry of the human primary visual cortex (V1). Here we demonstrate that much of this variability is due to extrinsic geometric features of the cortical folds, and that the intrinsic shape of V1 is similar across individuals. V1 was imaged in ten ex vivo human hemispheres using high-resolution (200 μm) structural magnetic resonance imaging at high field strength (7 T). Manual tracings of the stria of Gennari were used to construct a surface representation, which was computationally flattened into the plane with minimal metric distortion. The instrinsic shape of V1 was determined from the boundary of the planar representation of the stria. An ellipse provided a simple parametric shape model that was a good approximation to the boundary of flattened V1. The aspect ration of the best-fitting ellipse was found to be consistent across subject, with a mean of 1.85 and standard deviation of 0.12. Optimal rigid alignment of size-normalized V1 produced greater overlap than that achieved by previous studies using different registration methods. A shape analysis of published macaque data indicated that the intrinsic shape of macaque V1 is also stereotyped, and similar to the human V1 shape. Previoud measurements of the functional boundary of V1 in human and macaque are in close agreement with these results.

A Neural Model of Visually Guided Steering, Obstacle Avoidance, and Route Selection

A neural model is developed to explain how humans can approach a goal object on foot while steering around obstacles to avoid collisions in a cluttered environment. The model uses optic flow from a 3D virtual reality environment to determine the position of objects based on motion discotinuities, and computes heading direction, or the direction of self-motion, from global optic flow. The cortical representation of heading interacts with the representations of a goal and obstacles such that the goal acts as an attractor of heading, while obstacles act as repellers. In addition the model maintains fixation on the goal object by generating smooth pursuit eye movements. Eye rotations can distort the optic flow field, complicating heading perception, and the model uses extraretinal signals to correct for this distortion and accurately represent heading. The model explains how motion processing mechanisms in cortical areas MT, MST, and VIP can be used to guide steering. The model quantitatively simulates human psychophysical data about visually-guided steering, obstacle avoidance, and route selection.

A Neural Model of Visually Guided Steering, Obstacle Avoidance, and Route Selection

A neural model is developed to explain how humans can approach a goal object on foot while steering around obstacles to avoid collisions in a cluttered environment. The model uses optic flow from a 3D virtual reality environment to determine the position of objects based on motion discontinuities, and computes heading direction, or the direction of self-motion, from global optic flow. The cortical representation of heading interacts with the representations of a goal and obstacles such that the goal acts as an attractor of heading, while obstacles act as repellers. In addition the model maintains fixation on the goal object by generating smooth pursuit eye movements. Eye rotations can distort the optic flow field, complicating heading perception, and the model uses extraretinal signals to correct for this distortion and accurately represent heading. The model explains how motion processing mechanisms in cortical areas MT, MST, and posterior parietal cortex can be used to guide steering. The model quantitatively simulates human psychophysical data about visually-guided steering, obstacle avoidance, and route selection.

Shear effects on the properties and separation characteristics of whey protein precipitates

Selective isoelectric whey protein precipitation and aggregation is carried out at laboratory scale in a standard configuration batch agitation vessel. Geometric scale-up of this operation is implemented on the basis of constant impeller power input per unit volume and subsequent clarification is achieved by high speed disc-stack centrifugation. Particle size and fractal geometry are important in achieving efficient separation while aggregates need to be strong enough to resist the more extreme levels of shear that are encountered during processing, for example through pumps, valves and at the centrifuge inlet zone. This study investigates how impeller agitation intensity and ageing time affect aggregate size, strength, fractal dimension and hindered settling rate at laboratory scale in order to determine conditions conducive for improved separation. Particle strength is measured by observing the effects of subjecting aggregates to moderate and high levels of process shear in a capillary rig and through a partially open ball-valve respectively. The protein precipitate yield is also investigated with respect to ageing time and impeller agitation intensity. A pilot scale study is undertaken to investigate scale-up and how agitation vessel shear affects centrifugal separation efficiency. Laboratory scale studies show that precipitates subject to higher impeller shear-rates during the addition of the precipitation agent are smaller but more compact than those subject to lower impeller agitation and are better able to resist turbulent breakage. They are thus more likely to provide a better feed for more efficient centrifugal separation. Protein precipitation yield improves significantly with ageing, and 50 minutes of ageing is required to obtain a 70 - 80% yield of α-lactalbumin. Geometric scale-up of the agitation vessel at constant power per unit volume results in aggregates of broadly similar size exhibiting similar trends but with some differences due to the absence of dynamic similarity due to longer circulation time and higher tip speed in the larger vessel. Disc stack centrifuge clarification efficiency curves show aggregates formed at higher shear-rates separate more efficiently, in accordance with laboratory scale projections. Exposure of aggregates to highly turbulent conditions, even for short exposure times, can lead to a large reduction in particle size. Thus, improving separation efficiencies can be achieved by the identification of high shear zones in a centrifugal process and the subsequent elimination or amelioration of such.

Research and design of high-speed advanced analogue front-ends for fibre-optic transmission systems

In the last decade, we have witnessed the emergence of large, warehouse-scale data centres which have enabled new internet-based software applications such as cloud computing, search engines, social media, e-government etc. Such data centres consist of large collections of servers interconnected using short-reach (reach up to a few hundred meters) optical interconnect. Today, transceivers for these applications achieve up to 100Gb/s by multiplexing 10x 10Gb/s or 4x 25Gb/s channels. In the near future however, data centre operators have expressed a need for optical links which can support 400Gb/s up to 1Tb/s. The crucial challenge is to achieve this in the same footprint (same transceiver module) and with similar power consumption as today’s technology. Straightforward scaling of the currently used space or wavelength division multiplexing may be difficult to achieve: indeed a 1Tb/s transceiver would require integration of 40 VCSELs (vertical cavity surface emitting laser diode, widely used for short‐reach optical interconnect), 40 photodiodes and the electronics operating at 25Gb/s in the same module as today’s 100Gb/s transceiver. Pushing the bit rate on such links beyond today’s commercially available 100Gb/s/fibre will require new generations of VCSELs and their driver and receiver electronics. This work looks into a number of state‐of-the-art technologies and investigates their performance restraints and recommends different set of designs, specifically targeting multilevel modulation formats. Several methods to extend the bandwidth using deep submicron (65nm and 28nm) CMOS technology are explored in this work, while also maintaining a focus upon reducing power consumption and chip area. The techniques used were pre-emphasis in rising and falling edges of the signal and bandwidth extensions by inductive peaking and different local feedback techniques. These techniques have been applied to a transmitter and receiver developed for advanced modulation formats such as PAM-4 (4 level pulse amplitude modulation). Such modulation format can increase the throughput per individual channel, which helps to overcome the challenges mentioned above to realize 400Gb/s to 1Tb/s transceivers.