南海北部大陆边缘盆地热演化史及其与油气关系

This thesis focuses on the present-day thermal field features, evolution and their connections to hydrocarbon generation of the three continental margin basins-the Yinggehai (Yingge Sea), Qiongdongnan(southeast Qiong), and Pear River Mouth basins-in northern South China Sea, based on available data from drillings, loggings, seismic cross-sections, BHTs, thermal indicators (Ro%, inclusion, etc) and geopressure measurements. After studying of present-day distribution of geothermal field and thermal disturbance of fluid in the sedimentary strata, the author discovered that the distribution of gas fields in Yinggehai Basin are closely related to the distribution of anomalously high thermal gradient area, whereas it is not the case for the Pear River Mouse Basin. And detailed processing of the fluid inclusion data indicates that geothermal fluids activated frequently in this area, and they may mainly be derived upward from the overpressure and hydrocarbon-generating beds, 3000-4500 m in depth. Therefore, the abnormal gradients in sedimentary beds were mainly caused by the active geothermal fluids related to hydrocarbon migrating and accumulating in this area. Because of the effect of overpressure retarding on vitrinite reflectance, the thermal indicators for thermal history reconstruction should be assessed before put into use. Although some factors, such as different types of kerogen, heating ratio, activities of thermal fluids and overpressure, may have effects on the vitrinite reflectance, under the circumstance that thermal fluids and overpressure co-exist, overpressure retarding is dominant. And the depth and correction method of overpressure retarding were also determined in this paper. On the basis of reviewing the methods of thermal history studies as well as existing problems, the author believes that the combination of thermal-indicator-inversion and tectono-thermal modeling is an effective method of the thermal history reconstruction for sedimentary basins. Also, a software BaTherMod for modeling thermal history of basins was successfully developed in this work. The Yinggehai Basin has been active since Tertiary, and this was obviously due to its tectonic position-the plate transition zone. Under the background of high thermal flow, long-term quick subsidence and fluid activities were the main reasons that lead to high temperature and overpressure in this basin. The Zhujiangkou Basin, a Tertiary fault-basin within the circum-Pacific tectonic realm, was tectonically controlled by the motion of the Pacific Plate and resembles the other petroliferous basins in eastern China. This basin developed early, and characterized intensive extension in the early stage and weak activity in the later stage of its development. Whereas the Qiongdongnan Basin was in a weak extension early and intensity of extension increased gradually. The relative geographical locations and the extensional histories of three basins ilustrate that the northern continental margin of South China Sea spread from south to north. On the other hand, the Qiongdongnan and Yinggehai Basins may have been controlled by the same tectonic regime since later Tertiary, whereas the Zhujiangkou Basin was not meaningfully influenced. So, the tectono-thermal evolution character of the Qiongdonnan basin should be closely to the other two. It may be concluded that the three basins have been developed within the active continental margin since Tertiary, and the local lithosphere might undergo intensive extension-perhaps two or three times of episodic extension occurred. Extension lead to large tectonoc subsidence and extreme thick Tertiary sediments for hydrocarbon generation in the basins. In response to the periodic extension of the basins, the palaeothermal flow were also periodical. The three basins all have the characteristics of multi-phase thermal evolutions that is good for oil-gas generation. And the overpressure expands the depth range of oil-gas habitat, which is meaningful to petroleum exploration in this region.

Dynamic multi-partitioning for parallel finite element applications

The central product of the DRAMA (Dynamic Re-Allocation of Meshes for parallel Finite Element Applications) project is a library comprising a variety of tools for dynamic re-partitioning of unstructured Finite Element (FE) applications. The input to the DRAMA library is the computational mesh, and corresponding costs, partitioned into sub-domains. The core library functions then perform a parallel computation of a mesh re-allocation that will re-balance the costs based on the DRAMA cost model. We discuss the basic features of this cost model, which allows a general approach to load identification, modelling and imbalance minimisation. Results from crash simulations are presented which show the necessity for multi-phase/multi-constraint partitioning components

Heat and mass transfer in two-phase porous materials under intensive microwave heating

Computational results for the intensive microwave heating of porous materials are presented in this work. A multi-phase porous media model has been developed to predict the heating mechanism. Combined finite difference time-domain and finite volume methods were used to solve equations that describe the electromagnetic field and heat and mass transfer in porous media. The coupling between the two schemes is through a change in dielectric properties which were assumed to be dependent both on temperature and moisture content. The model was able to reflect the evolution of both temperature and moisture fields as well as energy penetration as the moisture in the porous medium evaporates. Moisture movement results from internal pressure gradients produced by the internal heating and phase change.

Hydraulics of stepped chutes: experimental-numerical-theoretical study

O estudo de escoamentos turbulentos em descarregadores em degraus tem sido um desafio para os investigadores. A macro-rugosidade do leito, a ondulação da superfície livre, a intermitência da localização da secção inicial de entrada de ar e o escoamento bi-fásico a jusante da secção inicial de entrada de ar fazem com que a caracterização do escoamento deslizante sobre turbilhões em descarregadores em degraus não seja simples. Actualmente, é possível combinar técnicas de medição fiáveis com simulações numéricas e análise teórica. Nesta dissertação, o estudo experimental baseia-se em resultados experimentais obtidos em duas instalações experimentais: a instalação A, do Laboratório Nacional de Engenharia Civil (LNEC), dotada de um descarregador em degraus com declive de 1V:0.75H e a instalação B, do Instituto Superior Técnico (IST), dotada de um descarregador em degraus com declive de 1V:2H. Uma sonda de ar, um tubo de Pitot modificado e vários hidrómetros permitiram o estudo do escoamento deslizante sobre turbilhões mono- e bifásico ao longo do descarregador em degraus da instalação A. Na bacia de dissipação de energia, a observação visual do escoamento foi auxiliada pelas leituras efectuadas em tomadas de pressão localizadas na soleira da bacia e pelas leituras da altura do escoamento efectuadas em réguas graduadas localizadas nas paredes da bacia. Na instalação B foram utilizados um tubo de Pitot e vários hidrómetros para estudar a região não arejada do escoamento deslizante sobre turbilhões no descarregador em degraus. O estudo numérico é baseado em simulações numéricas da região não arejada do escoamento deslizante sobre turbilhões sobre descarregadores em degraus com o código comercial de CFD FLOW-3D® de modo a reproduzir as condições ensaiadas experimentalmente. As simulações beneficiaram da técnica de blocos múltiplos (multi-block) num sistema de coordenadas cartesianas, da determinação da superfície livre pelo método TruVOF e da utilização de dois modelos de turbulência: os modelos k- e RNG k-. Por último, o estudo teórico consistiu em desenvolver um modelo simplificado 1D para determinar as características hidráulicas principais do trecho não arejado do escoamento deslizante sobre turbilhões em descarregadores em degraus. O modelo foi desenvolvido a partir das equações de Navier-Stokes, conjuntamente com resultados experimentais e numéricos. Os resultados apresentados nesta dissertação contribuem para o conhecimento do escoamento deslizante sobre turbilhões em descarregadores em degraus, nomeadamente na região não arejada, na secção inicial de entrada de ar e na região arejada. A hidráulica dos dissipadores de energia a jusante de descarregadores em degraus, em particular das bacias tipo III do USBR, é também objecto de estudo. Em relação à região não arejada do escoamento em descarregadores com declive acentuado, e com base em resultados experimentais e numéricos, são propostas expressões para estimar o desenvolvimento da altura equivalente de água, da espessura da camada limite, da concentração média de ar, do coeficiente de energia cinética, da dissipação de energia, do factor de resistência e do coeficiente n da fórmula de Manning. São ainda propostas expressões adimensionais para a energia cinética turbulenta e sua dissipação. Para declives moderados, são propostas expressões para estimar o desenvolvimento da altura equivalente de água, do coeficiente de energia cinética e da energia específica residual. São ainda apresentados valores do expoente 1/N da expressão adimensional da distribuição de velocidades, quer para descarregadores com declive acentuado quer com declive moderado. Em conformidade com outros estudos centrados em escoamentos de parede e com derivações teóricas, para a região não arejada do escoamento em descarregadores em degraus de acentuado declive, observa-se que o factor de resistência depende da macro-rugosidade criada pelos degraus e da geometria da secção transversal e que o coeficiente n da fórmula de Manning aumenta com a rugosidade. A descrição estatística da turbulência do escoamento é igualmente explorada, contribuindo para o conhecimento da estrutura do escoamento. Observou-se que para números de Reynolds rugoso não superiores a 6.8x104 a energia cinética turbulência e a sua dissipação cumprem leis de semelhança. Estas expressões adimensionais estão de acordo com os resultados obtidos por outros autores para escoamentos completamente desenvolvidos em canais abertos e no escoamento em rios com leito de gravilha. Em acréscimo, a taxa de dissipação de energia, quer para descarregadores de declive acentuado quer de moderado declive, é baixa. Por último, observa-se que os valores da média temporal da concentração de ar entre 0 e 1 medidos na região não arejada do escoamento dizem respeito não só ao ar capturado entre ondas de água, na zona de ondulação da superfície livre, mas também ao ar emulsionado no escoamento, i.e., sob a forma de bolhas de ar, quando perto da secção média inicial de entrada de ar, devido à diferença entre localizações instantânea e média temporal. Foram revistas metodologias e fórmulas para estimar a localização da secção inicial de entrada de ar e apresentadas expressões para estimar a concentração média de ar e a altura equivalente de água nessa secção. Relativamente à região de escoamento arejado em descarregadores em degraus com declive acentuado, os resultados experimentais apresentados nesta dissertação permitiram estimar a influência da definição da superfície livre nos parâmetros hidráulicos da região do escoamento arejado e estimar a máxima elevação do escoamento nesta região do escoamento. Com base nos resultados experimentais obtidos na bacia de dissipação de energia do tipo III do USBR localizada a jusante do descarregador em degraus da instalação A, observou-se que os perfis da altura piezométrica e da altura do escoamento tendem a seguir o perfil recomendado pelo USBR para bacias tipo III. A excepção ocorre à entrada da bacia, onde as alturas piezométricas apresentadas nesta dissertação excedem largamente as apresentadas pelo USBR. É ainda observado que, tal como entre as bacias tipo I e tipo III do USBR, o ressalto hidráulico estabiliza muito mais rapidamente numa bacia tipo III a jusante de um descarregador em degraus do que uma bacia tipo I a jusante do mesmo descarregador em degraus. Finalmente, observa-se que os blocos de amortecimento a colocar no descarregador não têm influência visível nos resultados da altura piezométrica nem da altura do escoamento ao longo da bacia. Relativamente às simulações numéricas do escoamento não arejado, a proximidade entre resultados experimentais e numéricos permite validar o modelo teórico e a integração numérica usados no FLOW-3D®. As simulações desenvolvidas também mostraram que o modelo de turbulência k- permite representar as características do escoamento não arejado em descarregadores em degraus, uma vez que não foram observadas diferenças significativas entre as simulações com este modelo e com o modelo RNG k-. Finalmente, observou-se que o modelo de entrada de ar usado no FLOW-3D® é válido para estimar a localização da secção inicial de entrada de ar. Por último, a proximidade entre os resultados obtidos da aplicação do modelo teórico desenvolvido no âmbito desta dissertação e os resultados experimentais indica que as hipóteses e simplificações consideradas no desenvolvimento do modelo são adequadas.

Automated two-dimensional separation flow system with electrochemical preconcentration, stripping, capillary electrophoresis and contactless conductivity detection for trace metal ion analysis

This paper describes the automation of a fully electrochemical system for preconcentration, cleanup, separation and detection, comprising the hyphenation of a thin layer electrochemical flow cell with CE coupled with contactless conductivity detection (CE-C(4)D). Traces of heavy metal ions were extracted from the pulsed-flowing sample and accumulated on a glassy carbon working electrode by electroreduction for some minutes. Anodic stripping of the accumulated metals was synchronized with hydrodynamic injection into the capillary. The effect of the angle of the slant polished tip of the CE capillary and its orientation against the working electrode in the electrochemical preconcentration (EPC) flow cell and of the accumulation time were studied, aiming at maximum CE-C(4)D signal enhancement. After 6 min of EPC, enhancement factors close to 50 times were obtained for thallium, lead, cadmium and copper ions, and about 16 for zinc ions. Limits of detection below 25 nmol/L were estimated for all target analytes but zinc. A second separation dimension was added to the CE separation capabilities by staircase scanning of the potentiostatic deposition and/or stripping potentials of metal ions, as implemented with the EPC-CE-C(4)D flow system. A matrix exchange between the deposition and stripping steps, highly valuable for sample cleanup, can be straightforwardly programmed with the multi-pumping flow management system. The automated simultaneous determination of the traces of five accumulable heavy metals together with four non-accumulated alkaline and alkaline earth metals in a single run was demonstrated, to highlight the potentiality of the system.

Pressure drop correlation for oil-refrigerant R134a mixture flashing flow in a small diameter tube

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Numerical analysis of refrigerant flow along non-adiabatic capillary tubes using a two-fluid model

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

A transmission line model developed directly in phase domain

The phases of a transmission line are tightly coupled due to mutual impedances and admittances of the line. One way to accomplish the calculations of currents and voltages in multi-phase lines consists in representing them in modal domain, where its n coupled phases are represented by their n propagation modes. The separation line in their modes of propagation is through the use of a modal transformation matrix whose columns are eigenvectors associated with the parameters of the line. Usually, this matrix is achieved through numerical methods which do not allow the achievement of an analytical model for line developed directly in the phases domain. This work will show an analytical model for phase currents and voltages of the line and results it will be applied to a hypothetical two-phase. It will be shown results obtained with that will be compared to results obtained using a classical model. © 2012 IEEE.

A transmission line model developed directly in phase domain

The phases of a transmission line are tightly coupled due to mutual impedances and admittances of the line. One way to accomplish the calculations of currents and voltages in multi phase lines consists in representing them in modal domain, where its n coupled phases are represented by their n propagation modes. The separation line in their modes of propagation is through the use of a modal transformation matrix whose columns are eigenvectors associated with the parameters of the line. Usually, this matrix is achieved through numerical methods which do not allow the achievement of an analytical model for line developed directly in the phases domain. This work will show an analytical model for phase currents and voltages of the line and results it will be applied to a hypothetical two-phase. It will be shown results obtained with that will be compared to results obtained using a classical model © 2003-2012 IEEE.

A distributed-parameters transmission line model developed directly in the phase domain

This article shows a transmission line model developed directly in the phase domain. The proposed model is based on the relationships between the phase currents and voltages at both the sending and receiving ends of a single-phase line. These relationships, established using an ABCD matrix, were extended to multi-phase lines. The proposed model was validated by using it to represent a transmission line during short-and open-circuit tests. The results obtained with the proposed model were compared with results obtained with a classical model based on modal decomposition. These comparisons show that proposed model was correctly developed. © 2013 Taylor and Francis Group, LLC.

A two-fluid model for refrigerant flow through adiabatic capillary tubes

This work presents a numerical model to simulate refrigerant flow through capillary tubes, commonly used as expansion devices in refrigeration systems. The flow is divided in a single-phase region, where the refrigerant is in the subcooled liquid state, and a region of two-phase flow. The capillary tube is considered straight and horizontal. The flow is taken as one-dimensional and adiabatic. Steady-state condition is also assumed and the metastable flow phenomena are neglected. The two-fluid model, considering the hydrodynamic and thermal non-equilibrium between the liquid and vapor phases, is applied to the two-phase flow region. Comparisons are made with experimental measurements of the mass flow rate and pressure distribution along two capillary tubes working with refrigerant R-134a in different operating conditions. The results indicate that the present model provides a better estimation than the commonly employed homogeneous model. Some computational results referring to the quality, void fraction, velocities, and temperatures of each phase are presented and discussed.

Saturated flow boiling heat transfer and critical heat flux in small horizontal flattened tubes

This paper presents experimental results for flow boiling heat transfer coefficient and critical heat flux (CHF) in small flattened tubes. The tested flattened tubes have the same equivalent internal diameter of 2.2 mm, but different aspect height/width ratios (H/W) of 1/4, 1/2, 2 and 4. The experimental data were compared against results for circular tubes using R134a and R245fa as working fluids at a nominal saturation temperature of 31 degrees C. For mass velocities higher than 200 kg/m(2)s, the flattened and circular tubes presented similar heat transfer coefficients. Such a behavior is related to the fact that stratification effects are negligible under conditions of higher mass velocities. Heat transfer correlations from the literature, usually developed using only circular-channel experimental data, predicted the flattened tube results for mass velocities higher than 200 kg/m(2)s with mean absolute error lower than 20% using the equivalent diameter to account for the geometry effect. Similarly, the critical heat flux results were found to be independent of the tube aspect ratio when the same equivalent length was kept. Equivalent length is a new parameter which takes into account the channel heat transfer area. The CHF correlations for round tubes predicted the flattened tube data relatively well when using the equivalent diameter and length. Furthermore, a new proposed CHF correlation predicted the present flattened tube data with a mean absolute error of 5%. (C) 2012 Elsevier Ltd. All rights reserved.

Geometrical and kinematic properties of interfacial waves in stratified oil-water flow in inclined pipe

The stratified oil-water flow pattern is common in the petroleum industry, especially in offshore directional wells and pipelines. Previous studies have shown that the phenomenon of flow pattern transition in stratified flow can be related to the interfacial wave structure (problem of hydrodynamic instability). The study of the wavy stratified flow pattern requires the characterization of the interfacial wave properties, i.e., average shape, celerity and geometric properties (amplitude and wavelength) as a function of holdup, inclination angle and phases' relative velocity. However, the data available in the literature on wavy stratified flow is scanty, especially in inclined pipes and when oil is viscous. This paper presents new geometric and kinematic interfacial wave properties as a function of a proposed two-phase Froude number in the wavy-stratified liquid-liquid flow. The experimental work was conducted in a glass test line of 12 m and 0.026 m id., oil (density and viscosity of 828 kg/m(3) and 0.3 Pa s at 20 degrees C, respectively) and water as the working fluids at several inclinations from horizontal (-20 degrees, -10 degrees, 0 degrees, 10 degrees, 20 degrees). The results suggest a physical relation between wave shape and the hydrodynamic stability of the stratified liquid-liquid flow pattern. (C) 2011 Elsevier Inc. All rights reserved.

Scherverhalten kolloidaler Lagenphasen in Rohrströmungen

Kolloidale Suspensionen aus identischen kugelförmigen, geladenen Partikeln in wässrigen Medien stellen ein ideales Modellsystem zur Untersuchung des Gleichgewichtsverhaltens, aber auch des Nicht-Gleichgewichtsverhaltens Weicher Materie dar. So bilden derartige Systeme bei hinreichend starker und langreichweitiger elektrostatischer Repulsion fluid und kristallin geordnete Strukturen aus, die wegen der weitreichenden Analogie zu atomar kondensierter Materie als kolloidale Fluide und Kristalle bezeichnet werden. Von großem Vorteil ist dabei die Möglichkeit zur kontrollierten Einstellung der Wechselwirkung und die gute optische Zugänglichkeit für Mikroskopie und Lichtstreuung sowie die Weichheit der Materialien, aufgrund derer sich auch Zustände fernab des mechanischen Gleichgewichts gezielt präparieren lassen. Themenstellung der vorliegenden Arbeit ist die Untersuchung des Phasenverhaltens und der Fließmechanismen kolloidaler Kristalle in einer Rohrströmung. Im ersten Teil der Arbeit wird gezeigt, dass beim Fluss durch eine zylindrische Röhre Mehrphasenkoexistenz auftritt, wobei ein polykristalliner Kern von einer isotropen Scherschmelze umgeben ist. Zusätzlich treten an der Grenze zwischen diesen Phasen und an der Rohrwand Phasen hexagonal geordneter übereinander hinweggleitender Lagen auf. Der Vergleich zwischen auf der Basis der Navier-Stokes-Gleichung theoretisch berechneten und gemessenen Geschwindigkeitsprofilen zeigt, dass jede dieser Phasen für sich Newtonsches Fließverhalten aufweist. Die Gesamtviskosität ist hingegen durch die mit dem Durchsatz veränderliche Phasenzusammensetzung Nicht-Newtonsch. Damit gelang es, die erstmalig von Würth beschriebene Scherverdünnung auf eine Veränderung der Phasenzusammensetzung zurückzuführen. Im zweiten Teil der Arbeit wurde erstmals das Fließverhalten der Lagenphasen mittels Lichtstreuung und Korrelationsanalyse untersucht. Dafür wurde ein im Prinzip einfacher, aber leistungsstarker Aufbau realisiert, der es erlaubt, die zeitliche Veränderung der Bragg-Reflexe der Lagenphase in radialer und azimutaler Richtung zu verfolgen und mittels Fourieranalyse zu analysieren. In Abhängigkeit vom Durchsatz geht die zunächst rastend gleitende Lagenphase in eine frei gleitende Lagenphase über, wobei charakteristische Veränderungen der Spektren sowie der Korrelationsfunktionen auftreten, die detailliert diskutiert werden. Der Übergang im Gleitmechanismus ist mit einem Verlust der Autokorrelation der Rotationskomponente der periodischen Intra-Lagenverzerrung verbunden, während die Kompressionskomponente erhalten bleibt. Bei hohen Durchflüssen lassen die Reflexbewegungen auf das Auftreten einer Eigenschwingung der frei gleitenden Lagen schließen. Diese Schwingung lässt sich als Rotationsbewegung, gekoppelt mit einer transversalen Auslenkung in Vortexrichtung, beschreiben. Die Ergebnisse erlauben eine detaillierte Diskussion von verschiedenen Modellvorstellungen anderer Autoren.

Thermal fluctuations and boundary conditions in the lattice Boltzmann method

The lattice Boltzmann method is a popular approach for simulating hydrodynamic interactions in soft matter and complex fluids. The solvent is represented on a discrete lattice whose nodes are populated by particle distributions that propagate on the discrete links between the nodes and undergo local collisions. On large length and time scales, the microdynamics leads to a hydrodynamic flow field that satisfies the Navier-Stokes equation. In this thesis, several extensions to the lattice Boltzmann method are developed. In complex fluids, for example suspensions, Brownian motion of the solutes is of paramount importance. However, it can not be simulated with the original lattice Boltzmann method because the dynamics is completely deterministic. It is possible, though, to introduce thermal fluctuations in order to reproduce the equations of fluctuating hydrodynamics. In this work, a generalized lattice gas model is used to systematically derive the fluctuating lattice Boltzmann equation from statistical mechanics principles. The stochastic part of the dynamics is interpreted as a Monte Carlo process, which is then required to satisfy the condition of detailed balance. This leads to an expression for the thermal fluctuations which implies that it is essential to thermalize all degrees of freedom of the system, including the kinetic modes. The new formalism guarantees that the fluctuating lattice Boltzmann equation is simultaneously consistent with both fluctuating hydrodynamics and statistical mechanics. This establishes a foundation for future extensions, such as the treatment of multi-phase and thermal flows. An important range of applications for the lattice Boltzmann method is formed by microfluidics. Fostered by the "lab-on-a-chip" paradigm, there is an increasing need for computer simulations which are able to complement the achievements of theory and experiment. Microfluidic systems are characterized by a large surface-to-volume ratio and, therefore, boundary conditions are of special relevance. On the microscale, the standard no-slip boundary condition used in hydrodynamics has to be replaced by a slip boundary condition. In this work, a boundary condition for lattice Boltzmann is constructed that allows the slip length to be tuned by a single model parameter. Furthermore, a conceptually new approach for constructing boundary conditions is explored, where the reduced symmetry at the boundary is explicitly incorporated into the lattice model. The lattice Boltzmann method is systematically extended to the reduced symmetry model. In the case of a Poiseuille flow in a plane channel, it is shown that a special choice of the collision operator is required to reproduce the correct flow profile. This systematic approach sheds light on the consequences of the reduced symmetry at the boundary and leads to a deeper understanding of boundary conditions in the lattice Boltzmann method. This can help to develop improved boundary conditions that lead to more accurate simulation results.