准噶尔盆地低阻油藏实验研究及储层评价方法

Low resistivity reservoir is a special reservoir which is different from normal reservoir in identification and evaluation.Through core experiment and analysis, the achievement of which resistivity is resulted from clay additive electric conductivities and high bound water saturation in Junggar basin is gained. For accurately evaluating low resistivity, a good many of experiment have been completed, such as resistivity index and formation factor in hi^jher temperature and higher pressure, semi-permeability board, cation exchange, bound water, NMR (nucleus magnetism response), non-Nad water in different temperature and salinity, the experiments result show that lower resistivity has complex relation with these electric-parameters and chloric ion content in non-NaCl water.Based on comprehensive interpretation of NMR and normal resistivity data, the volume of moved water, bound water, moved oil and residual oil in the strata can be determined quantitatively and which have significant influence on reservoir recognition and perforation optimized.Experiment data (SEM mold, thin section, X ray diffraction, mercury penetration) can be used to analysis low resistivity forming and the relation between low resistivity and pore texture, to set up relation between porosity, permeability and petrophysical property. The reservoir was sorted, evaluated and described. The oil bedding in southern margin of Junggar basin is low porosity, low resistivity reservoir.Based on invasion theory of electric well-logging, modelling and inversion of resistivity well-logging are accomplished. For enhancing low resistivity resulted from higher bound water saturation and cation exchange, invasion period, invasion radius, the relation between fluid distribution in pore and response of laterolog logging have been studied. Virgin zone resistivity, invasion zone resistivity and invasion radius were inversed and which enhanced evaluation accuracy of reservoir. The method was used to process well-logging data in Luliang oilfield and southern margin in Junggar basin, and reservoir resistivity was enhanced effectively, appropriate oil saturation gained and it has better effect on oil exploration.

黄骅坳陷高尘头地区储层非均质性及成藏研究

Gaochentou region is located in the southwest direction of Gaochentou village in Huanghua city of Hebei province. In regionally structural position, It lies in Qikou sag In the middle part of Huanghua depression, which belongs to the east part of the south Dagang structure zone in the middle part of Huanghua depression. Its' very beneficial at regional structure in Gaochentou , and It becomes the advantage area for oil and gas gathered and preserved, Sandstone reservoir of Dongying Formation is main bearing bed .Dongying Formation in Gaochentou region of Huanghua depression is consisted of set of mudstone and sandstone interbeds by deposited delta fades . Dongying Formation can be divided into 3 members from above to below: the first member of Dongying Formation (FMDF), the second member of Dongying Formation (SMDF), and third member of Dongying Formation (TMDF). The lithology of the upper part of FMDF was consisted of mostly middle-grained and fine-grained sandstone, and it is small for the oil-bearing area of the sand bodies .The lithology of the lower part is coarse-grained sandstone bodies which are well connected between sandstone bodies of wells, and the lower part was main bed of oil production in Dongying Formation; SMDF and TMDF are consisted of larger scale set of mudstone, in which the sandbodies are lenticular and pinch out quickly, and the lithology was mostly fine sandstone and silt stone, in which there are little oil and gas .Because the reservoirs in this area are largely influenced by the factors such as lithology, fault and others, and the reservoirs have the strong,heterogeneity , there exists the problem of oil-down and water-up for vertical distribution of oil and gas bearing. It is not very clearly for the three dimension distribution of sandstone , and the geology researchs is not enough. So, it can't satisfy the need of further development and production for Gaochentou oilfield.Having the key problem of oil-down and water-up and the mechanism of the reservoir for Gaochentou area, There are as follow study works, the first, is study of the high-resolution correlation of sequence stratigraphy and sedimentary microfacies. Dongying Formation was divided into three parasequence sets and each parasequence set was divided into different amount of parasequences. FMDF, as the main oil and gas producing bed, can be divided into seven parasequences. Oil and gas are discovered in six parasequences except the seventh. On the basis of study of sedimentary microfacies, the sediments of Dongying Formation are considered deposited mainly in delta front subfacies. The microfacies types of Dongying Formation are sub-water distirbutary channel, sub-water natural bank, inter distributary channel bay, distributary channel mouth dam, and delta front mat sand.Seismic facies analysis and logging-constrained inversion technique were applied by Author for transverse prediction of sandstone reservoir. Having 4 modes of interwell single sandbodies correlation technique, Author have described distribution characteristics of sandbodies, and established geological reservoir model of Gaochentou reservoir.Author presented that the reservoirs characteristic have very strong heterogeneity ,and In the section of sandstone interlayed with mudstone,the folium sandstone interlayed with each other, and the wedge shaped sandbodies pinched out in the mudstone. So the pinch-out up sandstone trap and lenticular sandstone trap are easily formed. They are most small scale overlying pinches out in the place of slope. This article applies the concept of deep basin oil to resolve reasonably the problem of which the oil is below the water in Gaochentou area. Combined with the study of sedimentary facies, reservoir and other aspects, the mechanism and patterns of deep basin oil are studied on the basis of characteristics in Gaochentou area.On the basis of the above study, the mechanism of the oil and gas' migration and accumulation in isotropic sandstone and heterogeneous sandstone are thoroughly analyzed through experiments on physical modeling. Experiments on physical modeling show that the discrepancy between sand layers with different permeability and thickness has important influence on the direction, path, and injection layer of oil's migration. At the beginning of the injection of oil and gas in high permeability sand layer, the pressure is low, the migration resistance is small, and the oil and gas are more easily displacing the water in sand. So it can act as good transformation layer or reservoir. But at the beginning of the injection of oil and gas in sand layer with low permeability, the pressure is high, the migration resistance is big, and the oil and gas are more difficultly displacing the water in sand. So it can only act as bad or worse transformation layer or reservoir. Even if it cannot act as transformation layer or reservoir, it can act as water layer or dry layer. The discrepancy between sand layers on permeability and thickness can make discrepancy in injection of oil and gas between different layers. Consequently it leads to small amount of oil and gas injection in sand layers with low permeability. Ultimately it affects the oil's accumulation and distribution in different sand layers.At Last, combining analysis of the structure and pool forming condition, The thesis has established models of reservoir formation to predict the advantage distribution of oil and gas bearing , and put forward the prospective target It is not only of theoretical signification for explosion and importance, but also has realistic value in guiding the progressive petroleum exploration and exploitation.

地震反射界面形态与介质各向异性属性联立成像：地震反射波旅行时与偏振信息综合反演

It has been a difficult problem faced by seismologists for long time that how exactly to reconstruct the earth's geometric structure and distribution of physical attributes according to seismic wave's kinematical and dynamic characteristics, obtained in seismological observation. The jointing imaging of seismic reflector and anisotropy attributes in the earth interior is becoming the research hot spot. The limitation of shoot and observation system makes that the obtained seismic data are too scarce to exactly reconstruct the geological objects. It is popular that utilizing only seismic reflection traveltimes or polarizations information make inversion of the earth's velocity distribution by fixing seismic reflector configuration (vice versa), these will lead to the serious non-uniqueness reconstruction due to short of effective data, the non-uniqueness problem of reconstructing anisotropy attributes will be more serious than in isotropy media. Obviously it is not enough to restrict the media structure only by information of seismic reflection traveltimes or polarizations, which even sometimes will lead to distorted images and misinterpretation of subsurface structure. So we try to rebuild seismic reflection structure (geometry) and media anisotropic structure (physics) in the earth interior by jointing data of seismic wave kinematics and dynamics characteristics, we carry out the new experiment step by step, and the research mainly comprises of two parts: one is the reconstruction of P-wave vertical velocity and anisotropic structure(Thomsen parameter s and 8) in the transversely isotropic media with vertical symmetrical axis(VTI) by fixing geometrical structure, and the other is the simultaneous inversion of the reflector surface conformation and seismic anisotropic structure by jointing seismic reflection traveltimes and polarizations data. Simulated annealing method is used to the first research part, linear inversion based on BG theory and Simulated annealing are applied to the second one. All the research methods are checked by model experiments, then applied to the real data of the wide-angle seismic profile from Tunxi, Anhui Province, to Wenzhou, Zhejiang Province. The results are as following The inversion results based on jointing seismic PP-wave or PSV-wavereflection traveltimes and polarizations data are more close to real model than themodels based simply on one of the two data respectively. It is shown that the methodwe present here can effectively reconstruct the anisotropy attributes in the earth'sinterior when seismic reflector structure is fixed.The layer thickness, P-wave vertical velocity and Thomsen anisotropicparameters {s and 8) could be resolved simultaneously by jointing inversion ofseismic reflection traveltimes and polarizations with the linear inversion methodbased on BG theory.The image of the reflector structure, P-wave vertical velocity and theanisotropy parameters in the crust could be obtained from the wide-angle seismicprofile from Tunxi (in Anhui Province), to Wenzhou (in Zhejiang Province). Theresults reveal the difference of the reflector geometrical structure and physicalattributes in the crust between Yangtze block and Cathaysia block, and attempt tounderstand the characteristics of the crustal stress field in the areas.

地幔对流理论及其在解释地球物理观测的应用

With the improving of mantle convection theory, the developing of computing method and increasing of the measurement data, we can numerically simulate more clearly about the effects on some geophysical observed phenomenons such as the global heat flow and global lithospheric stress field in the Earth's surface caused by mantle convection, which is the primary mechanism for the transport of heat from the Earth's deep interior to its surface and the underlying force mechanism of dynamics in the Earth.Chapter 1 reviews the historical background and present research state of mantle convection theory.In Chapter 2, the basic conception of thermal convection and the basic theory about mantle flow.The effects on generation and distribution of global lithospheric stres s field induced by mantle flow are the subject of Chapter 3. Mantle convection causes normal stress and tangential stresses at the bottom of the lithosphere, and then the sublithospheric stress field induces the lithospheric deformation as sixrface force and results in the stress field within the lithosphere. The simulation shows that the agreement between predictions and observations is good in most regions. Most of subduction zones and continental collisions are under compressive. While ocean ridges, such as the east Pacific ridge, the Atlantic ridge and the east African rift valley, are under tensile. And most of the hotspots preferentially occur in regions where calculated stress is tensile. The calculated directions of the most compressive principal horizontal stress are largely in accord with that of the observation except for some regions such as the NW-Pacifie subduction zone and Qinghai-Tibet Plateau, in which the directions of the most compressive principal horizontal stress are different. It shows that the mantel flow plays an important role in causing or affecting the large-scale stress field within the lithosphere.The global heat flow simulation based on a kinematic model of mantle convection is given in Chapter 4. Mantle convection velocities are calculated based on the internal loading theory at first, the velocity field is used as the input to solve the thermal problem. Results show that calculated depth derivatives of the near surface temperature are closely correlated to the observed surface heat flow pattern. Higher heat flow values around midocean ridge systems can be reproduced very well. The predicted average temperature as a function of function of depth reveals that there are two thermal boundary layers, one is close to the surface and another is close to the core-mantle boundary, the rest of the mantle is nearly isothermal. Although, in most of the mantle, advection dominates the heat transfer, the conductive heat transfer is still locally important in the boundary layers and plays an important role for the surface heat flow pattern. The existence of surface plates is responsible for the long wavelength surface heat flow pattern.In Chapter 5, the effects on present-day crustal movement in the China Mainland resulted from the mantle convection are introduced. Using a dynamic method, we present a quantitative model for the present-day crustal movement in China. We consider not only the effect of the India-Eurasia collision, the gravitational potential energy difference of the Tibet Plateau, but also the contribution of the shear traction on the bottom of the lithosphere induced by the global mantle convection. The comparison between our results and the velocity field obtained from the GPS observation shows that our model satisfactorily reproduces the general picture of crustal deformation in China. Numerical modeling results reveal that the stress field on the base of the lithosphere induced by the mantle flow is probably a considerable factor that causes the movement and deformation of the lithosphere in continental China with its eflfcet focuing on the Eastern China A numerical research on the small-scale convection with variable viscosity in the upper mantle is introduced in Chapter 6. Based on a two-dimensional model, small-scale convection in the mantle-lithosphere system with variable viscosity is researched by using of finite element method. Variation of viscosity in exponential form with temperature is considered in this paper The results show that if viscosity is strongly temperature-dependent, the upper part of the system does not take a share in the convection and a stagnant lid, which is identified as lithosphere, is formed on the top of system because of low temperature and high viscosity. The calculated surface heat flow, topography and gravity anomaly are associated well with the convection pattern, namely, the regions with high heat flow and uplift correspond to the upwelling flow, and vice versa.In Chapter 7, we give a brief of future research subject: The inversion of lateral density heterogeneity in the mantle by minimizing the viscous dissipation.

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

In exploration seismology, the geologic target of oil and gas reservoir in complex medium request the high accuracy image of the structure and lithology of the medium. So the study of the prestack image and the elastic inversion of seismic wave in the complex medium come to the leading edge. The seismic response measured at the surface carries two fundamental pieces of information: the propagation effects of the medium and the reflections from the different layer boundaries in the medium. The propagation represent the low-wavenumber component of the medium, it is so-called the trend or macro layering, whereas the reflections represent the high-wavenumber component of the medium, it is called the detailed or fine layering. The result of migration velocity analysis is the resolution of the low-wavenumber component of the medium, but the prestack elastic inversion provided the resolution of the high-wavvenumber component the medium. In the dissertation, the two aspects about the migration velocity estimation and the elastic inversion have been studied.Firstly, any migration velocity analysis methods must include two basic elements: the criterion that tell us how to know whether the model parameters are correct and the updating that tell us how to update the model parameters when they are incorrect, which are effected on the properties and efficiency of the velocity estimation method. In the dissertation, a migration velocity analysis method based on the CFP technology has been presented in which the strategy of the top-down layer stripping approach are adapted to avoid the difficult of the selecting reduce .The proposed method has a advantage that the travel time errors obtained from the DTS panel are defined directly in time which is the difference with the method based on common image gather in which the residual curvature measured in depth should be converted to travel time errors.In the proposed migration velocity analysis method, the four aspects have been improved as follow:? The new parameterization of velocity model is provided in which the boundaries of layers are interpolated with the cubic spline of the control location and the velocity with a layer may change along with lateral position but the value is calculated as a segmented linear function of the velocity of the lateral control points. The proposed parameterization is suitable to updating procedure.? The analytical formulas to represent the travel time errors and the model parameters updates in the t-p domain are derived under local lateral homogeneous. The velocity estimations are iteratively computed as parametric inversion. The zero differential time shift in the DTS panel for each layer show the convergence of the velocity estimation.? The method of building initial model using the priori information is provided to improve the efficiency of velocity analysis. In the proposed method, Picking interesting events in the stacked section to define the boundaries of the layers and the results of conventional velocity analysis are used to define the velocity value of the layers? An interactive integrate software environment with the migration velocity analysis and prestack migration is built.The proposed method is firstly used to the synthetic data. The results of velocity estimation show both properties and efficiency of the velocity estimation are very good.The proposed method is also used to the field data which is the marine data set. In this example, the prestack and poststack depth migration of the data are completed using the different velocity models built with different method. The comparison between them shows that the model from the proposed method is better and improves obviously the quality of migration.In terms of the theoretical method of expressing a multi-variable function by products of single-variable functions which is suggested by Song Jian (2001), the separable expression of one-way wave operator has been studied. A optimization approximation with separable expression of the one-way wave operator is presented which easily deal with the lateral change of velocity in space and wave number domain respectively and has good approach accuracy. A new prestack depth migration algorithm based on the optimization approximation separable expression is developed and used to testing the results of velocity estimation.Secondly, according to the theory of the seismic wave reflection and transmission, the change of the amplitude via the incident angle is related to the elasticity of medium in the subsurface two-side. In the conventional inversion with poststack datum, only the information of the reflection operator at the zero incident angles can be used. If the more robust resolutions are requested, the amplitudes of all incident angles should be used.A natural separable expression of the reflection/transmission operator is represented, which is the sum of the products of two group functions. One group function vary with phase space whereas other group function is related to elastic parameters of the medium and geological structure.By employing the natural separable expression of the reflection/transmission operator, the method of seismic wave modeling with the one-way wave equation is developed to model the primary reflected waves, it is adapt to a certain extent heterogeneous media and confirms the accuracy of AVA of the reflections when the incident angle is less than 45'. The computational efficiency of the scheme is greatly high.The natural separable expression of the reflection/transmission operator is also used to construct prestack elastic inversion algorithm. Being different from the AVO analysis and inversion in which the angle gathers formed during the prstack migration are used, the proposed algorithm construct a linear equations during the prestack migration by the separable expression of the reflection/transmission operator. The unknowns of the linear equations are related to the elasticity of the medium, so the resolutions of them provided the elastic information of the medium.The proposed method of inversion is the same as AVO inversion in , the difference between them is only the method processing the amplitude via the incident angle and computational domain.

地震偏移宽角吸收边界条件和Toeplitz矩阵反演中的边界效应研究

The content of this paper is based on the research work while the author took part in the key project of NSFC and the key project of Knowledge Innovation of CAS. The whole paper is expanded by introduction of the inevitable boundary problem during seismic migration and inversion. Boundary problem is a popular issue in seismic data processing. At the presence of artificial boundary, reflected wave which does not exist in reality comes to presence when the incident seismic wave arrives at the artificial boundary. That will interfere the propagation of seismic wave and cause alias information on the processed profile. Furthermore, the quality of the whole seismic profile will decrease and the subsequent work will fail.This paper has also made a review on the development of seismic migration, expatiated temporary seismic migration status and predicted the possible break through. Aiming at the absorbing boundary problem in migration, we have deduced the wide angle absorbing boundary condition and made a compare with the boundary effect of Toepiitz matrix fast approximate computation.During the process of fast approximate inversion computation of Toepiitz system, we have introduced the pre-conditioned conjugate gradient method employing co circulant extension to construct pre-conditioned matrix. Especially, employment of combined preconditioner will reduce the boundary effect during computation.Comparing the boundary problem in seismic migration with that in Toepiitz matrix inversion we find that the change of boundary condition will lead to the change of coefficient matrix eigenvalues and the change of coefficient matrix eigenvalues will cause boundary effect. In this paper, the author has made an qualitative analysis of the relationship between the coefficient matrix eigenvalues and the boundary effect. Quantitative analysis is worthy of further research.

Study of the carbonaceous deposits formed on a Mo/HZSM-5 catalyst in methane dehydro-aromatization by using TG and temperature-programmed techniques

Carbonaceous deposits formed during the temperature-programmed surface reaction (TPSR) of methane dehydro-aromatization (MDA) over Mo/HZSM-5 catalysts have been investigated by TPH, TPCO2 and TPO, in combination with thermal gravimetric analysis (TG). The TPO profiles of the coked catalyst after TPSR of MDA show two temperature peaks: one is at about 776 K and the other at about 865 K. The succeeding TPH experiments only resulted in the diminishing of the area of the high-temperature peak, and had no effect on the area of the low-temperature peak. On the other hand, the TPO profiles of the coked catalyst after succeeding TPCO2 experiments exhibited obvious reduction in the areas of both the high-and low-temperature peaks, particularly in the area of the low-temperature peak. On the basis of TPSR, TPR and TPCO2 experiments and the corresponding TG analysis, quantitative analysis of the coke and the kinetics of its burning-off process have been studied. (C) 2001 Elsevier Science B.V. All rights reserved.

Highly ordered periodic mesoporous ethanesilica synthesized under neutral conditions

Highly ordered mesoporous ethanesilica (MES) with 2D hexagonal structure was synthesized from 1,2-bis(trimethoxysilyl) ethane under neutral conditions for the first time. Divalent salts, such as NiCl2, MgCl2, ZnCl2, ZnSO4 and Zn(NO3)(2), were used to help the formation of the ordered mesostructure. The MES samples were characterized by powder X-ray diffraction, nitrogen sorption, transmission electron microscopy, FT-IR, C-13 and Si-29 solid-state NMR and thermal gravimetric analysis. A phase transition from a disordered wormhole-like structure to an ordered P6mm structure was observed upon the addition of inorganic salts. The pore size of the MES decreases from 4.7 to 3.9 nm with increasing content of the inorganic salts. Fluoride was also found to be important for the formation of ordered MES under neutral conditions.

Una testimonianza silenziosa.Storia della Chiesa cattolica in Tunisia dal Trattato del Bardo alla ‘rivoluzione dei gelsomini’

The work examines the change involving the Church in Tunisia from the period of the Protectorate to the present through the fundamental moments of independence (1956) and the signing of the ‘Modus vivendi’ (1964). In the first structure of the “modern” Church, a fundamental role was played by the complex figure of the French Cardinal Charles-Allemand Lavigerie who, while giving strong impulse to setting up disinterested charitable social initiatives by the congregations (Pères Blancs, Soeurs Blanches and others), also represented the ideal of the ‘evangelizing’ (as well as colonial) Church which, despite its declared will to avoid proselytism, almost inevitably tended to slip into it. During the French Protectorate (1881-1956) the ecclesiastic institution concentrated strongly on itself, with little heed for the sensitivity of its host population, and developed its activities as if it were in a European country. From the social standpoint, the Church was mostly involved in teaching, which followed the French model, and health facilities. In the Church only the Pères Blancs missionaries were sincerely committed to promoting awareness of the local context and dialogue with the Muslims. The Catholic clergy in the country linked its religious activity close to the policy of the Protectorate, in the hope of succeeding in returning to the ancient “greatness of the African Church”, as the Eucharistic Congress in Carthage in 1930 made quite clear. The Congress itself planted the first seed in the twentyfive- year struggle that led the Tunisian population to independence in 1956 and the founding of the Republic in 1957. The conquest of independence and the ‘Modus vivendi’ marked a profound change in the situation and led to an inversion of roles: the Catholic community was given the right to exist only on the condition that it should not interfere in Tunisian society. The political project of Bourguiba, who led the Republic from 1957 to 1987, aimed to create a strongly egalitarian society, with a separation between political and religious powers. In particular, in referring to the Church, he appeared as a secularist with no hostility towards the Catholics who were, however, considered as “cooperators”, welcome so long as they were willing to place their skills at the service of the construction of the state. So, in the catholic Community was a tension between the will of being on the side of the country and that of conserving a certain distance from it and not being an integral part of it. In this process of reflection, the role of the Second Vatican Council was fundamental: it spread the idea of a Church open to the world and the other religions, in particular to Islam: the teaching of the Council led the congregations present in the country to accept the new condition. This new Church that emerged from the Council saw some important events in the process of “living together”, of “cultural mixing” and the search for a common ground between different realities. The almost contemporary arrival of Arab bishops raised awareness among the Tunisians of the existence of Christian Arabs and, at the same time, the Catholic community began considering their faith in a different way. In the last twenty years the situation has continued to change. Side by side with the priests present for decades or even those born there, some new congregations have begun to operate, albeit in small numbers: they have certainly revitalized the community of the faithful, but they sometimes appear more devoted to service “within” the Church, than to services for the population, and are thus characterized by exterior manifestations of their religion. This sort of presence has made it possible for Bourguiba's successor, Ben Ali (president from 1987 to 2011), to practice forms of tolerance even more clearly, but always limited to formal relations; the Tunisians are still far from having a real understanding of the Catholic reality, with certain exceptions connected to relations on a personal and not structured plane, as was the case in the previous period. The arrival of a good number of young people from sub-Saharan Africa, most of all students, belonging to the JCAT, and personnel of the BAD has “Africanized” the Church in Tunisia and has brought about an increase in Christians' exterior manifestations; but this is a visibility that is not blatant but discreet, with the implicit risk of the Church continuing to be perceived as a sort of exterior body, alien to the country; nor can we say, lacking proper documentation, how it will be possible to build a bridge between different cultures through the “accompaniment” of Christian wives of Tunisians. Today, the Church is living in a country that has less and less need of it; its presence, in the schools and in health facilities, is extremely reduced. And also in other sectors of social commitment, such as care for the disabled, the number of clergymen involved is quite small. The ‘revolution’ in 2011 and the later developments up to the present have brought about another socio-political change, characterized by a climate of greater freedom, but with as yet undefinable contours. This change in the political climate will inevitable have consequences in Tunisia’s approach to religious and cultural minorities, but it is far too soon to discuss this on the historical and scientific planes.

Generalized semi-geostrophic theory on a sphere

Douglas, Robert; Cullen, M.J.P.; Roulston, I.; Sewell, M.J., (2005) 'Generalized semi-geostrophic theory on a sphere', Journal of Fluid Mechanics 531 pp.123-157 RAE2008

Radio tomography: A new experimental technique

Pryse, Sian, 'Radio tomography: A new experimental technique', Surveys in Geophysics (2003) 24 pp.1-38 RAE2008

Modele fizyczne planetoid pojedynczych i podwójnych otrzymane na podstawie przeglądów fotometrycznych: przygotowanie do naukowego wykorzystania danych z kosmicznej misji Gaia

Wydział Fizyki

Tectonics Inversions & Emerald Deposits

In spite of the great amount of emerald deposits throughout the world, the priorities in quality and volume of extracted rough material are the sites of Colombia (Muzo and Chivor emerald belts). This sites are know even before the Spanish conquistadores. Emeralds were extracted from Somondoco mine (today Chivor) since 1537 and from Muzo in 1567. Contrariwise to the majority of the emerald deposits of the world, which are associated with granitic rocks, the Colombian emerald deposits are associated with hydrofracturing (the main factor controlling emerald mineralization) and hydrothermal fluids, rich in beryl, chrome and vanadium, induced by a tectonic inversion of the deep Mesozoic backarc basin, which is also responsible of the majority of the petroleum systems of the foredeep and foldbelt areas (maturation of the source-rocks andcreation of structural traps). The host rocks of the emeralds are carbonaceous calsiltites (calcareous schists) rich in organic matter of Lower Cretaceous age, which are cut by calcite veins, which, often, contain emeralds, particularly when they are folded. Indeed, since long time (Cheilletz, A. and Giulliani, G., 1996) suggested a two-stage model for the formation of the Colombian emeralds : (i) Stage I is characterized by décollement planes (early compressional tectonic regime) within the carbonaceous calsiltites, hydrothermal fluid infiltration and wall-rock metasomatic alteration ; (ii) Stage II (late tectonic regime) deforms the previous veins by thrust-related folds (development of stratiform and hydraulic breccia), which are synchronous of the emerald mineralization. The resulting tectonic structures are complex fold patterns characterized by propagation anticlines with emerald veins and emerald hydraulic breccia in the apexes, as in Quipama, Tendenquema and Chivor mines. Otherwise stated, since all emerald exploitations are, presently underground, exhaustive geological and particularly structural studies are required to reduce the probability of disappointments. The color of emeralds is from light green to thick green with obvious pleochroism. They appears with different colors when observed at different angles, especially with polarized light. The emeralds from Coscuez deposits have a homogeneous intensive color and bluish tone. At Muzo deposit, the emeralds have middle or dark green color with yellowish tone. At the Chivor deposits, the emeralds have less intensive green color with slight bluish tone. The typical inclusions are albite and pyrite, as well as long bubbles with three phase-inclusions according the zones of growth and along the crystal shapes.

A study of stone fragmentation in shock wave lithotripsy by customizing the acoustic field and waveform shape

Shock wave lithotripsy is the preferred treatment modality for kidney stones in the United States. Despite clinical use for over twenty-five years, the mechanisms of stone fragmentation are still under debate. A piezoelectric array was employed to examine the effect of waveform shape and pressure distribution on stone fragmentation in lithotripsy. The array consisted of 170 elements placed on the inner surface of a 15 cm-radius spherical cap. Each element was driven independently using a 170 individual pulsers, each capable of generating 1.2 kV. The acoustic field was characterized using a fiber optic probe hydrophone with a bandwidth of 30 MHz and a spatial resolution of 100 μm. When all elements were driven simultaneously, the focal waveform was a shock wave with peak pressures p+ =65±3MPa and p−=−16±2MPa and the −6 dB focal region was 13 mm long and 2 mm wide. The delay for each element was the only control parameter for customizing the acoustic field and waveform shape, which was done with the aim of investigating the hypothesized mechanisms of stone fragmentation such as spallation, shear, squeezing, and cavitation. The acoustic field customization was achieved by employing the angular spectrum approach for modeling the forward wave propagation and regression of least square errors to determine the optimal set of delays. Results from the acoustic field customization routine and its implications on stone fragmentation will be discussed.

Articulated Pose Estimation in a Learned Smooth Space of Feasible Solutions

A learning based framework is proposed for estimating human body pose from a single image. Given a differentiable function that maps from pose space to image feature space, the goal is to invert the process: estimate the pose given only image features. The inversion is an ill-posed problem as the inverse mapping is a one to many process. Hence multiple solutions exist, and it is desirable to restrict the solution space to a smaller subset of feasible solutions. For example, not all human body poses are feasible due to anthropometric constraints. Since the space of feasible solutions may not admit a closed form description, the proposed framework seeks to exploit machine learning techniques to learn an approximation that is smoothly parameterized over such a space. One such technique is Gaussian Process Latent Variable Modelling. Scaled conjugate gradient is then used find the best matching pose in the space of feasible solutions when given an input image. The formulation allows easy incorporation of various constraints, e.g. temporal consistency and anthropometric constraints. The performance of the proposed approach is evaluated in the task of upper-body pose estimation from silhouettes and compared with the Specialized Mapping Architecture. The estimation accuracy of the Specialized Mapping Architecture is at least one standard deviation worse than the proposed approach in the experiments with synthetic data. In experiments with real video of humans performing gestures, the proposed approach produces qualitatively better estimation results.