温度、压力耦合背景下多相流体运移的数学模型研究

The Mathematical modeling of multiphase fluid flow is an important aspect of basin simulation, and also is a topic of geological frontier. Based on coupling relation of temperature, pressure and fluid flow, this dissertation discusses the modeling which conform to geological regularities of fluid migration. The modeling that is multi-field and multiphase includes heat transport equation, pressure evolvement equation, solution transport equation and fluid transport equation. The finite element method is effective numerical calculation methods. Author applies it to solve modeling and implements the finite element program, and the modeling is applied to Ying-Qiong Basin. The channels of fluid vertical migration are fault, fracture and other high penetrability area. In this thesis, parallel fracture model and columnar channel model have been discussed, and a characteristic time content and a characteristic space content been obtained to illustrate the influences of stratigraphic and hydrodynamic factors on the process. The elliptoid fracture model is established and its approximately solution in theory is gotten. Three kinds of modeling are applied to analyze the transient variation process of fluid pressure in the connected permeable formations. The elliptoid fracture model is the most similar geology model comparing with the other fracture models so the research on this fracture model can enhance the understanding to fluid pressure. In the non-hydrodynamic condition, because of the difference between water density and nature gas density, nature gas can migrate upon by float force. A one-dimension mathematical model of nature gas migration by float force is established and also applied to analyze the change in the saturation of gas. In the process of gas migration its saturation is non-continuous. Fluid flow is an important factor which influences the distribution of the temperature-field, the change of temperature can influence fluid property (including density, viscidity, and solubility),a nd the temperature field has coupling relations to the fluid pressure field. In this dissertation one-dimension and two-dimension thermal convection modeling is developed and also applied to analyze convective and conductive heat transfer. Author has established one-dimension and two-dimension mathematical modeling in which fluid is a mixture of water and nature gas based on the coupling relation between temperature and pressure, discussed mixture fluid convection heat transfer in different gas saturation, and analyzed overpressure form mechanism. Based on geothermal abnormity and pore pressure distribution in Dongfong 1-1, Yinggehai Basin, South China Sea, one-dimension mathematical modeling of coupling temperature and pressure is established. The modeling simulates the process that fluid migrates from deep to shallow and overpressure forms in shallow. When overpressure is so large that fractures appear and overpressure is released. As deep fluid flow to shallow, the high geothermal then forms in shallow. Based on the geological characteristics in Ya13-1, two-dimension mathematical modeling of coupling temperature and pressure is established. Fluid vertically flows in fault and then laterally migrates in reservoir. The modeling simulates the geothermal abnormity and pore pressure distribution in reservoir.

地震波场模拟及叠前深度偏移的李群方法

The processes of seismic wave propagation in phase space and one way wave extrapolation in frequency-space domain, if without dissipation, are essentially transformation under the action of one parameter Lie groups. Consequently, the numerical calculation methods of the propagation ought to be Lie group transformation too, which is known as Lie group method. After a fruitful study on the fast methods in matrix inversion, some of the Lie group methods in seismic numerical modeling and depth migration are presented here. Firstly the Lie group description and method of seismic wave propagation in phase space is proposed, which is, in other words, symplectic group description and method for seismic wave propagation, since symplectic group is a Lie subgroup and symplectic method is a special Lie group method. Under the frame of Hamiltonian, the propagation of seismic wave is a symplectic group transformation with one parameter and consequently, the numerical calculation methods of the propagation ought to be symplectic method. After discrete the wave field in time and phase space, many explicit, implicit and leap-frog symplectic schemes are deduced for numerical modeling. Compared to symplectic schemes, Finite difference (FD) method is an approximate of symplectic method. Consequently, explicit, implicit and leap-frog symplectic schemes and FD method are applied in the same conditions to get a wave field in constant velocity model, a synthetic model and Marmousi model. The result illustrates the potential power of the symplectic methods. As an application, symplectic method is employed to give synthetic seismic record of Qinghai foothills model. Another application is the development of Ray+symplectic reverse-time migration method. To make a reasonable balance between the computational efficiency and accuracy, we combine the multi-valued wave field & Green function algorithm with symplectic reverse time migration and thus develop a new ray+wave equation prestack depth migration method. Marmousi model data and Qinghai foothills model data are processed here. The result shows that our method is a better alternative to ray migration for complex structure imaging. Similarly, the extrapolation of one way wave in frequency-space domain is a Lie group transformation with one parameter Z and consequently, the numerical calculation methods of the extrapolation ought to be Lie group methods. After discrete the wave field in depth and space, the Lie group transformation has the form of matrix exponential and each approximation of it gives a Lie group algorithm. Though Pade symmetrical series approximation of matrix exponential gives a extrapolation method which is traditionally regarded as implicit FD migration, it benefits the theoretic and applying study of seismic imaging for it represent the depth extrapolation and migration method in a entirely different way. While, the technique of coordinates of second kind for the approximation of the matrix exponential begins a new way to develop migration operator. The inversion of matrix plays a vital role in the numerical migration method given by Pade symmetrical series approximation. The matrix has a Toepelitz structure with a helical boundary condition and is easy to inverse with LU decomposition. A efficient LU decomposition method is spectral factorization. That is, after the minimum phase correlative function of each array of matrix had be given by a spectral factorization method, all of the functions are arranged in a position according to its former location to get a lower triangular matrix. The major merit of LU decomposition with spectral factorization (SF Decomposition) is its efficiency in dealing with a large number of matrixes. After the setup of a table of the spectral factorization results of each array of matrix, the SF decomposition can give the lower triangular matrix by reading the table. However, the relationship among arrays is ignored in this method, which brings errors in decomposition method. Especially for numerical calculation in complex model, the errors is fatal. Direct elimination method can give the exact LU decomposition But even it is simplified in our case, the large number of decomposition cost unendurable computer time. A hybrid method is proposed here, which combines spectral factorization with direct elimination. Its decomposition errors is 10 times little than that of spectral factorization, and its decomposition speed is quite faster than that of direct elimination, especially in dealing with a large number of matrix. With the hybrid method, the 3D implicit migration can be expected to apply on real seismic data. Finally, the impulse response of 3D implicit migration operator is presented.

鄂尔多斯盆地苏里格气田储层次生孔隙性质与分布预测

Permian reservoir in Sulige area of Ordos Basin, on which this paper focused, belongs to fluvial-delta lithofacies. The majority formations in this area are complicated channel sand deposit with serious inhomogeneity which makes natural gas exploration be very tough in this area. This inhomogeneity can be found everywhere both in large horizontal area and vertical profile of inner and interbedded formations.This paper studied the inhomogeneity characteristic of Permian formation in sulige area of Ordos Basin according to the logging data.Correlating with core data, a criterion to distinguish different type of reservoirs by using logging data is determined after the study of logging response is done considering the diverse conditions of deposit environments, lithology and reservoir space. The characteristic relationships between the various type formations and logging responses fully and systemically are established.It investigated reservoir parameter calculation methods amply. Combining the conventional and special logging data, basing on the feature of low porosity -permeability formation of sulige area, a set of methods to calculate reservoir parameters was formed including primary porosity, secondary porosity, fracture porosity, permeability and water saturation under the conditions of both low porosity-permeability and inhomogenous reservoirs. One thing should be pay close attention is the parameter M for calculating saturation. It is found that the M in low porosity -permeability formation decreases as the porosity decrease, which is opposite to the law that M increases as the porosity decrease in the formation with intermediate to high porosity and permeability. This view has innovated the traditional theory and offered theory basis for the logging interpretation of low porosity - permeability reservoir. Meanwhile it also improved the Arqi formula theoretically and enhanced the logging interpretation accuracy and rescued a number of formations which has been thought to be hopeless according to the old theory.By using advantage logging interpretation procedure, a territorial synthetic geology evaluation to the inhomogeneous reservoir was completed basing on the single well interpretation. All the reservoir evaluation parameters including sand formation thickness, primary porosity, secondary porosity were calculated and evaluated. The rules of changing and development for sand formation thickness, sand physical properties and secondary porous were found at different formations of upper part of the Member 8 of Shihezi, lower part of the Member 8 of Shihezi, the Member 1 of shanxi and the Member 2 of shanxi individually. Evaluation and Correlation of these five formations were also completed and one conclusion was arrived: upper part of the Member 8 of Shihezi formation has the best performance followed by the lower part of the Member 8 of Shihezi, the Member 1 of shanxi and the Member 2 of shanxi formation.After studied the relationship between reservoir deposition characteristic and the natural gas richness, it is regarded that reservoir inhomogeneity is the key issue of the impaction on the natural gas. Natural gas in Sulige gas field was mainly accumulated in sands of channel bar, distributary channel and debouchure bar. Especially, the quartz sand with rich of secondary porous space has obvious better physical properties than other reservoir and usually can forms the concentration of natural gas.

Desenvolvimento de um carro elétrico puro: instrumentação e registos energéticos

Desde meados do século XX os avanços na indústria automóvel trazem a associação da eletrónica sendo esta cada vez mais necessária. Com este trabalho pretende-se construir um computador de bordo de um veículo elétrico dando continuação a uma dissertação anterior, que consistiu em um Controlador de Motor de carro elétrico com capacidade de travagem regenerativa. Este Computador deve ser de baixo custo e é destinado à instrumentação, apresentação de balanços energéticos e funções básicas de diagnóstico do estado do veículo. Foi igualmente um objetivo desenvolver uma Interface a partir de um LCD para este Computador de Bordo. No modelo de Controlador de Motor disponível foram aplicadas alterações e adições na instrumentação de forma a conseguir um diagnóstico do veículo mais preciso e extenso. Para o computador ter conhecimento do consumo e outros parâmetros relacionados com o motor e respetivo Controlador, foi necessário estabelecer uma comunicação entre ambas unidades. Antes de se implementar um protocolo de comunicação realizou-se uma pesquisa por protocolos usados na indústria automóvel, com o intuito de saber qual o mais apropriado para o presente trabalho. A tarefa seguinte consistiu em uma pesquisa por hardware com o qual desenvolver o Computador e sua Interface. Os balanços energéticos implicaram o desenvolvimento de métodos de cálculo, efetuados com os parâmetros transmitidos pelo Controlador de Motor. No final da dissertação demonstra-se todas as funcionalidades do Computador de Bordo desenvolvido e como este é utilizado.

Procedimentos e dinâmicas inerentes ao processo dos concursos públicos e à realização de "empreitadas de obras públicas, em Portugal"

Trabalho de relatório de Estágio de natureza científica para obtenção do grau de Mestre em Engenharia Civil na Área de Especialização em Edificações

Projecto de fundações e estrutura de um edifício de escritórios utilizando os eurocódigos estruturais

Trabalho Final de Mestrado para obtenção do grau de Mestre em Engenharia Civil na Área de Especialização de Estruturas

Simulation model for driving dynamics, energy use and power supply

Dissertação para obtenção do grau de Mestre em Engenharia Electrotécnica Ramo de Energia

Estágio curricular na empresa Consulfrio

Trabalho Final de Mestrado para obtenção do grau de Mestre em Engenharia Mecânica

Controlo da fermentação do mosto e das condições ambientais de uma sala de lagares

Trabalho Final de Mestrado para obtenção do grau de Mestre em Engenharia Mecânica

Preços de transferência – o caso português

Dissertação de Mestrado apresentado ao Instituto de Contabilidade e Administração do Porto para a obtenção do grau de Mestre em Auditoria, sob orientação da Mestre Gabriela Pinheiro

Condicionamento acústico e certificação energética de edifícios

O presente relatório desenvolvido no âmbito do estágio curricular realizado na empresa – ASL & Associados, com duração de seis meses, tem como objetivo a elaboração e desenvolvimento de estudos acústicos de edifícios e elaboração de certificados energéticos de edifícios existentes. Inicialmente, será feita uma descrição da empresa acolhedora expondo a sua missão, valores e o portfólio abarcando projetos e revisão de projetos, certificação energética, ensaios acústicos, gestão e fiscalização de empreendimentos. As atividades desenvolvidas no estágio englobam a certificação energética de edifícios existentes aplicando a nova regulamentação térmica (REH), apresentando as simplificações de cálculo aplicadas a edifícios existentes e dois exemplos específicos de uma fração autónoma e um edifício unifamiliar. No campo da acústica foram estudados três edifícios, um habitacional, um misto e um de serviços. Para o condicionamento acústico dos casos de estudo serão aplicados métodos de cálculo utilizados na empresa e um método de cálculo para avaliação e classificação da qualidade acústica, o qual foi desenvolvido pelo Laboratório Nacional de Engenharia Civil – LNEC. Por fim, foram elaborados as ponderações sobre a experiência vivida, como considerações sobre os resultados obtidos tendo em conta os objetivos propostos para o desenvolvimento do estágio.

Modelação e dimensionamento de estruturas

O presente relatório refere-se ao trabalho desenvolvido durante o período de estágio curricular enquadrado no ciclo de estudos do Mestrado em Engenharia Civil do ISEP, Instituto Superior de Engenharia do Porto. O estágio desenvolveu-se ao longo de seis meses, desde Fevereiro até Julho de 2015, na empresa FASE S.A., em ambiente de gabinete, para a obtenção do grau de mestre em Engenharia Civil. Durante o período de estágio foram desenvolvidos projetos de estabilidade de estruturas metálicas e de betão armado, desde a fase de conceção, pré-dimensionamento, modelação numérica e análise de resultados, dimensionamento final, até à produção das peças desenhas e escritas constituintes de um projeto de estruturas. Foi possível analisar e dimensionar estruturas através de programas de cálculo automático e ferramentas de cálculo que serão referidas no presente relatório. Neste relatório será descrita pormenorizadamente a elaboração de um projeto em betão armado em todas as suas vertentes. Apresentam-se também as metodologias de cálculo empregues. No culminar, enunciam-se algumas conclusões de carácter geral decorrentes do trabalho desenvolvido.

El control geopolítico de los recursos naturales en Colombia: El caso del carbón en el Departamento del Cesar

El interés de este trabajo de investigación es explicar el tema del control geopolítico de los recursos naturales en Colombia en el periodo 2002-2011. Se analiza y explica cómo a partir del control territorial y el régimen de extracción desarrollados por la Drummond Company Inc. en el Departamento del Cesar se puede establecer su índice de poder geopolítico. Siguiendo los postulados de Michael Klare, respecto al desarrollo de una carrera internacional por los recursos energéticos mundiales, y los postulados de Kepa Sodupe, respecto a la posibilidad de reconfigurar los métodos de cálculo del poder, se avanza hacia los resultados de la investigación. Estos resultados permiten identificar de forma documentada el índice de poder geopolítico de una multinacional extranjera en el territorio colombiano.

Fine-root turnover rates of European forests revisited: an analysis of data from sequential coring and ingrowth cores

Background and Aims Forest trees directly contribute to carbon cycling in forest soils through the turnover of their fine roots. In this study we aimed to calculate root turnover rates of common European forest tree species and to compare them with most frequently published values. Methods We compiled available European data and applied various turnover rate calculation methods to the resulting database. We used Decision Matrix and Maximum-Minimum formula as suggested in the literature. Results Mean turnover rates obtained by the combination of sequential coring and Decision Matrix were 0.86 yr−1 for Fagus sylvatica and 0.88 yr−1 for Picea abies when maximum biomass data were used for the calculation, and 1.11 yr−1 for both species when mean biomass data were used. Using mean biomass rather than maximum resulted in about 30 % higher values of root turnover. Using the Decision Matrix to calculate turnover rate doubled the rates when compared to the Maximum-Minimum formula. The Decision Matrix, however, makes use of more input information than the Maximum-Minimum formula. Conclusions We propose that calculations using the Decision Matrix with mean biomass give the most reliable estimates of root turnover rates in European forests and should preferentially be used in models and C reporting.