Relationship between space distribution of the benthic macroinvertebrates community and trophic state in a Neotropical reservoir (Itupararanga, Brazil)

The purpose of this work was to verify the benthic macroinvertebrates community responses through environmental factors along a headwater tropical reservoir. Samplings were taken with a Van-Veen grab along the reservoir in littoral and profundal regions and in the headwater, next to the dam and the middle of the reservoir. Samples were taken during both wet and dry seasons. Dissolved oxygen concentrations, electric conductivity, temperature and pH near the sediment have been performed in situ, at every sampling station by using a multiprobe and Secchi disc. Total water phosphorus and chlorophyll a concentrations were analyzed to determine the trophic state index. Sediment's organic matter, total phosphorus, nitrogen concentrations and granulometric composition were measured. In order to verify which environmental variables would have more influence over the benthic macroinvertebrates community, a canonical correspondence analysis (CCA) was performed. The total number of recorded taxa was 28. Among them, the family Chironomidae (Diptera) was the richest group (19 taxa). It can be proposed that the benthic macroinvertebrates community may be influenced by environmental conditions such as nutrient and organic matter availability, as well as dissolved oxygen concentration. Macroinvertebrates are adequate bioindicators of water quality due to their sensibility to environmental changes mentioned before. Chironomus sp, Limnodrilus hoffmeisteri and Branchiura sowerbyi comprises a group that can be considered bio-indicators of eutrophic conditions. A second group can be considered as indicator of mesotrophic conditions. The presence of two or more members from that group which comprises Tanytarsini spp, Fissimentum sp, Pelomus sp and Goeldichironomus sp, like predominant taxa, may indicates mesotrophic conditions.

The incidents of housing policies in the process of nesting the urban space (Buenos Aires, 1870-1999)

This article describes the housing policies displayed historically in B.A. that affected the processes of configuration and of giving hierarchy to the urban space. Since the process of modernization of B.A., at the end of the XIXth century, housing measures, urbanistic projects and political decisions have influenced the building of the city and the space distribution of their inhabitants: in this way, they have integrated some people and excluded others. No wonder, that, historically, popular sectors have established themselves in the South of the city (and in outskirt villages). These zones have been disregarded by the state, which has invested less in these areas. We will see how the connection between state housing policies and the population redistribution in the city, confirms the persistence of a strong process of urban and residential segregation that tends to expel the inhabitants of popular sectors and to attract the ones of the middle-class and high class sectors.

CLASSICAL DISTRIBUTION-FUNCTIONS DERIVED FROM WIGNER DISTRIBUTION-FUNCTIONS

A mapping which relates the Wigner phase-space distribution function associated with a given stationary quantum-mechanical wavefunction to a specific solution of the time-independent Liouville transport equation is obtained. Two examples are studied.

Space analysis of water erosion in the red-yellow latosol under conilon coffee cultivation

This work aimed to study the space behavior of the water erosion in a red-yellow latosol. Then a study was developed in an area with colinon coffee cultivation in an Experimental Farm of Bananal do Norte of INCAPER in Cachoeiro de Itapemirim - ES. Soil samples were obtained from 0,0 to 0,20 m depth in an irregular grid with 109 samples. The analyzed variables were granulometric fractions, erodibility (K), natural erosion potential (PNE), soil loss (A) and erosion risk (RE). All the variables showed space dependency with moderate index of space dependency and similar standard of space distribution. The soil loss is related with the space distribution of the granulometric fractions.

Mapping Wigner distribution functions into semiclassical distribution functions

A mapping that relates the Wigner phase-space distribution function of a given stationary quantum mechani-cal wave function, a solution of the Schrödinger equation, to a specific solution of the Liouville equation, both subject to the same potential, is studied. By making this mapping, bound states are described by semiclassical distribution functions still depending on Planck's constant, whereas for elastic scattering of a particle by a potential they do not depend on it, the classical limit being reached in this case. Following this method, the mapped distributions of a particle bound in the Pöschl-Teller potential and also in a modified oscillator potential are obtained.

The stochastic exit problem for dynamical systems

The problem of "exit against a flow" for dynamical systems subject to small Gaussian white noise excitation is studied. Here the word "flow" refers to the behavior in phase space of the unperturbed system's state variables. "Exit against a flow" occurs if a perturbation causes the phase point to leave a phase space region within which it would normally be confined. In particular, there are two components of the problem of exit against a flow:

i) the mean exit time

ii) the phase-space distribution of exit locations.

When the noise perturbing the dynamical systems is small, the solution of each component of the problem of exit against a flow is, in general, the solution of a singularly perturbed, degenerate elliptic-parabolic boundary value problem.

Singular perturbation techniques are used to express the asymptotic solution in terms of an unknown parameter. The unknown parameter is determined using the solution of the adjoint boundary value problem.

The problem of exit against a flow for several dynamical systems of physical interest is considered, and the mean exit times and distributions of exit positions are calculated. The systems are then simulated numerically, using Monte Carlo techniques, in order to determine the validity of the asymptotic solutions.

中子相衬层析成像方法

介绍了中子相衬和中子层析的工作原理，提出了中子相衬层析成像方法．这种方法兼备了相衬和层析的各自优点，不采用干涉术，能高分辨率地重构待测样品的三维空间位相分布，能测试和分辨吸收系数非常相近样品的空间结构．

Mapeamento dos remanescentes de manguezal, a partir da interpretação de ortofotos coloridas, e análise dos vetores de pressão antrópica, como subsídio ao planejamento para conservação do ecossistema no Estado do Rio de Janeiro

Manguezal é um ecossistema costeiro que ocorre nas regiões tropicais e subtropicais do planeta, ocupando a zona entremarés dos oceanos, e sendo caracterizado pela presença de vegetação arbórea adaptada à condições adversas de salinidade, substrato, baixa oxigenação e submersão periódica. A pressão sobre os manguezais do Estado do Rio de Janeiro vem se intensificando nas últimas décadas, e estão associadas a vetores de pressão como os aterros, desmatamentos, queimadas, corte seletivo de madeira, captura predatória de moluscos e crustáceos, lançamento de efluentes de origens diversas, a superexplotação dos recursos pesqueiros e a utilização de técnicas e apetrechos inadequados. Considerando a inexistência de mapeamento integrado e atualizado dos remanescentes de manguezal, indicando sua localização e dimensionamento, o presente estudo veio suprir essa demanda, construindo uma ferramenta consistente para a análise dos principais vetores a que estão expostos, subsidiando a proposição de ações para a conservação e monitoramento desse ecossistema. Essas ações consideram a necessidade de preservação da biodiversidade, da manutenção da atividade pesqueira, da estabilidade da linha de costa, e da subsistência de diversas populações que habitam a região costeira. O mapeamento dos manguezais do Estado do Rio de Janeiro foi elaborado a partir da interpretação visual de ortofotografias coloridas do ano de 2005, na escala 1:10.000, tendo sido realizadas checagens de campo para identificação da verdade terrestre. Os remanescentes mapeados totalizam uma área de aproximadamente 17.720 ha, estando distribuídos por sete regiões hidrográficas localizadas na zona costeira fluminense. Esses ocorrem com mais freqüência, e com maiores dimensões, nas regiões da baía da Ilha Grande, Guandu(Sepetiba) e baía de Guanabara. O estudo contemplou ainda o levantamento e sistematização de dados cartográficos e de sensoriamento remoto, e a identificação e análise dos principais vetores de pressão que atuam sobre esses, a partir da adaptação da metodologia da Análise de Cadeia Causal. Nessa análise foram identificados como principais problemas ambientais dos manguezais fluminenses, a Modificação de habitats e comunidades, a Poluição, e a Exploração não sustentável dos recursos pesqueiros, todos associados aos diferentes vetores de pressão já relacionados. Por fim, foram apresentadas propostas de ações para subsidiar a implementação da Política Estadual para a Conservação dos Manguezais do Estado do Rio de Janeiro, contemplando os níveis operacional, de planejamento, e político. A reativação do Grupo Técnico Permanente sobre Manguezais é de vital importância para a retomada dessas discussões e para a implementação de ações, apoiado na ampliação dos conhecimentos sobre esse rico ecossistema e, integrando e fortalecendo a atuação dos diversos atores envolvidos, buscando assim garantir a integridade dos manguezais fluminenses

Study of crystal formation in titanate glass irradiated by 800 nm femtosecond laser pulse

The structure of the titanate glass is destroyed during irradiation by the femtosecond laser pulses, and (TiO6)(8-) and (TiO4)(4-) anion units are exsolved from the network of the titanate glass. These anion units are rearranged to form some crystals such as anatase and Ba2TiO4 crystals. By Raman spectroscopy, it is found that these crystals have a strong dependence on the intensity of the femtosecond laser pulses. The relation between the generation of these crystals and space distribution of the femtosecond laser power intensity is qualitatively explained. (c) 2005 Elsevier B.V. All rights reserved.

Joint modelling of voicing label and continuous F0 for HMM based speech synthesis

Fundamental frequency, or F0 is critical for high quality speech synthesis in HMM based speech synthesis. Traditionally, F0 values are considered to depend on a binary voicing decision such that they are continuous in voiced regions and undefined in unvoiced regions. Multi-space distribution HMM (MSDHMM) has been used for modelling the discontinuous F0. Recently, a continuous F0 modelling framework has been proposed and shown to be effective, where continuous F0 observations are assumed to always exist and voicing labels are explicitly modelled by an independent stream. In this paper, a refined continuous F0 modelling approach is proposed. Here, F0 values are assumed to be dependent on voicing labels and both are jointly modelled in a single stream. Due to the enforced dependency, the new method can effectively reduce the voicing classification error. Subjective listening tests also demonstrate that the new approach can yield significant improvements on the naturalness of the synthesised speech. A dynamic random unvoiced F0 generation method is also investigated. Experiments show that it has significant effect on the quality of synthesised speech. © 2011 IEEE.

西藏齿突蟾蝌蚪冷适应策略研究

本文通过对高海拔两栖类西藏齿突蟾(Scutiger boulengeri)蝌蚪在实验室特定低温条件下的冷适应微空间行为分布的动态变化分析、温度耐受性实验及在不同适应温度的乳酸脱氢酶(LDH)同工酶的酶量与活性比较分析, 探讨了高海拔两栖类蝌蚪的部分冷适应策略。 西藏齿突蟾蝌蚪在不同温度的行为分布是一连续、动态过程，需用多种检验方法综合利用才能进行判断；在15℃, 除低海拔分布的西藏齿突蟾种群外所有实验物种蝌蚪均符合负二项分布、NeymanⅡ型分布；在10℃, 高海拔两栖类蝌蚪均符合负二项分布、NeymanⅡ型分布；在5℃、0℃低温时，高海拔两栖类不同分组的西藏齿突蟾蝌蚪的负二项分布、NeymanⅡ型分布均呈现明显差异, 这可能与高海拔两栖类蝌蚪在低温条件下通过不断地改变其行为分布方式来避免自身被冻伤有关。野外观察表明：高海拔两栖类蝌蚪常选择与流动河水相连的静水水体这种微生境中生存, 蝌蚪应对环境温度极端变化会不断改变其行为分布方式来选择最佳生存温度以避免极端高、低温对自身身体的伤害, 这种对微生境的利用能力对高海拔两栖类蝌蚪耐受极端环境温度的变化极其重要。 两栖类蝌蚪的温度耐受性实验表明不同的驯化温度可以改变西藏齿突蟾蝌蚪、两栖类仙琴水蛙蝌蚪的最适温度、逃避温度，并具有显著影响。 随着驯化温度5℃、10℃逐渐升高, 其最适温度、逃避温度也在一定范围内升高，但驯化温度对低海拔的仙琴水蛙蝌蚪的最适温度、逃避温度的改变效应大于高海拔的西藏齿突蟾蝌蚪的改变效应, 仙琴水蛙蝌蚪对温度的耐受范围、最适温度和逃避温度的ARRS值都大于西藏齿突蟾蝌蚪, 这说明仙琴水蛙蝌蚪对环境温度变化的适应能力大于西藏齿突蟾蝌蚪。 高海拔地区不同分组的两栖类蝌蚪, 在0℃适应温度时, LDH5条带的酶相对含量最高，而在5℃、10℃、15℃适应温度时，LDH5条带的酶相对含量明显都降低, 这表明酵解作用是高海拔两栖类蝌蚪的一些组织在低温﹑缺氧环境中的重要供能方式。高海拔两栖类蝌蚪同一分组的LDH总酶活性总是表现为10℃适应温度的总酶活性最高，而对低海拔的两栖类蝌蚪则是0℃适应温度的总酶活性最高, 这说明高海拔两栖类蝌蚪的LDH同工酶A、B两亚基基因活性在10℃时最高, 而低海拔两栖类蝌蚪的LDH同工酶A、B两亚基基因活性在0℃时最高。同时发现在15℃适应温度组的高海拔两栖类蝌蚪的LDH电泳图谱都有第6条带,有可能由LDH - C亚基组成, 对高海拔两栖类蝌蚪的LDH - C亚基只在15℃适应温度下才表达的机理还有待进一步的研究。 高海拔两栖类西藏齿突蟾蝌蚪通过行为分布方式的改变来选择最佳的生存温度, 这种温度选择过程与野外特定的微生境的存在密切相关, 现在由于人类对河道的不合理利用正在导致高海拔两栖类蝌蚪赖以生存的这种微生境逐渐消失, 这种微生境的消失将加速高海拔的两栖类种群数量衰退的进程。高海拔两栖类物种蝌蚪在低温（0℃）上表现出的同工酶多谱带说明，其A、B两亚基都有所表达，及其参与代谢的方式也是正常的，而低海拔两栖类物种蝌蚪只有A亚基表达的LDH5存在，因此其主要参与酵解过程，这种通过动物自身生理代谢方式的改变来适应极端环境温度条件的变化是高海拔两栖类蝌蚪能适应低温环境的重要策略。但高海拔物种的适应温度变化范围显著小于低海拔物种，对环境温度的变化适应能力有限，特别是对高温区域，因此全球气候变化可能对高海拔物种影响更为显著。 The partly cold-adaptation stratagem of the high altitude amphibian tadpole were researched in the laboratory by analyzing the high altitude amphibian tadpole of Scutiger boulengeri mainly on endpoints related to the dynamic variation of the micro-spatial behavior distribution patterns, the experiment of the temperature tolerance, and the enzyme content and activity of the lactic acid dehydrogenase(LDH) isozyme in special temperature condition. The behavior distribution of the Scutiger boulengeri tadpole is continuous and variable, but it can be figured out by multple testing ways. At 15℃, all of the experiment amphibian tadpoles behavior distribution fit both for the negative binomial distribution and NeymanⅡtype distribution except for the low altitude Scutiger boulengeri tadpoles. At 10℃, all of the high altitude amphibian tadpoles behavior distribution fit both for the negative binomial distribution and NeymanⅡtype distribution. At lower temperature, 5℃ and 0℃, the high altitude amphibian tadpoles of the Scutiger boulengeri at different groups behavior distribution fit for or don’t fit for behavior distribution respectively. It is denoted that the high altitude amphibian tadpoles probably avoid frostbiting by varying the behavior distribution patterns at low temperature condition. The high altitude amphibian tadpoles often actively select the special microhabitat which has the connected still water body and the flowing water body in the wild. It is important that tadpoles can endure the extreme temperature variety in this kind of microhabitat, because tadpoles can be better survival through select temperature condition through migrating in these kinds of microhabitats by varying their own behavior distribution patterns. Different acclimation temperature causes the significant change of preferred temperature(PT)、 avoiding temperature(AT) both in high altitude amphibian Scutiger boulengeri tadpoles and in low altitude amphibian Rana daunchina tadpoles in the temperature endurance experiment. With the acclimation temperature growing from 5℃ to 10℃. the PT and the AT of them would be uprise to some extent, but the effect of acclimation temperature on the PT and the AT of the tadpoles of Rana daunchina is more significant than the ones on the tadpoles of Scutiger boulengeri, at the same, the effects on the temperature endurance range, the ARRs of the tadpoles of Rana daunchina would be stronger than the ones on the tadpoles of Scutiger boulengeri. It is implied that the adaptation ability of tadpoles of Rana daunchina to the surroundings temperature alternation preferred to tadpoles of Scutiger boulengeri. At 0℃ acclimation temperature, the LDH5 enzyme comparative content of the high altitude amphibian tadpoles at different groups was highest, but it becomes lower at 5℃、10℃、15℃ acclimation temperature. It indicated that the alcoholysis role was the important ways of applying energy for special tissue of the high altitude amphibian tadpoles in low-temperature and low-oxygen condition. The total enzyme activity of the LDH of the high altitude amphibian tadpoles in the same group always keeps the highest at 10℃ acclimation temperature, but the low altitude amphibian tadpoles’ was maximum at 0℃. It was denoted that the gene activity of LDH -A and LDH – B submit was highest at 10℃ acclimation temperature for the high altitude amphibian tadpoles, but the low altitude amphibian tadpoles’ was maximum at 0℃. Meanwhile, the LDH electrophoretogram of the high altitude amphibian tadpoles always composed of 6 stripes at 15℃ acclimation temperature，the extra stripe probably was composed by LDH-C submit。It is unknown why LDH-C expresses only under high temperature。. The high altitude amphibian tadpoles can select the most optimal temperature by changing their behavior distribution patterns ceaselessly, but this course of selecting the most suitable temperature correlated with the special microhabitat in the wild closely. Nowadays, this kind of microhabitat which the high altitude amphibian tadpoles rely on are lossing gradually for human being exploit the riverway unreasonably. The disappearing of the microhabitat would accelerate the decline of the high altitude amphibian population. Compare to one band of LDH5, which only composed by the LDH-A submit, presents in the low altitude amphibian at 0℃, the five bands which composed by the LDH-A and LDH-B are checked out, this means the species which occurred in the highland is more adaptable to the low temperature. It is an important stratagem for the high altitude amphibian tadpoles adapt to the limited low temperature depends on the animal energy metabolism change.However, this kind of adaption is restricted, the adaption range to the temperature is much norrow in the high altitude amphibian than in the low one, especially for the high temperature side. The global climate change will be more serious for the high altitude species.

叠前叠后联合反演技术及其在QHD32-6油田的应用研究

What geophysical inversion studied includes the common mathematics physical property of inversion and the constitution and appraisal method of solution in geophysics domain, i.e. using observed physical phenomenon from the earth surface to infer space changing and physical property structure of medium within the earth. Seismic inversion is a branch of geophysical inversion. The basic purpose of seismic inversion is to utilizing seismic wave propagating law in the medium underground to infer stratum structure and space distribution of physical property according to data acquisition, processing and interpretation, and then offer the vital foundation for exploratory development. Poststack inversion is convenient and swift, its acoustic impedance inversion product can reflect reservoir interior changing rule to a certain degree, but poststack data lack abundant amplitude and travel time information included in prestack data because of multiple superimpose and weaken the sensitiveness which reflecting reservoir property. Compared with poststack seismic inversion, prestack seismic inversion has better fidelity and more adequate information. Prestack seismic inversion, including waveform inversion, not only suitable for thin strata physical property inversion, it can also inverse reservoir oil-bearing ability. Prestack seismic inversion and prestack elastic impedance inversion maintain avo information, sufficiently applying seismic gathering data with different incident angle, partial angle stack, gradient and intercept seismic data cube. Prestack inversion and poststack inversion technology were studied in this dissertation. A joint inversion method which synthesize prestack elastic wave waveform inversion, prestack elastic impedance inversion and poststack inversion was proposed by making fully use of prestack inversion multiple information and relatively fast and steady characteristic of poststack inversion. Using the proposed method to extract rock physics attribute cube with clear physical significance and reflecting reservoir characterization, such as P-wave and S-wave impedance, P-wave and S-wave velocity, velocity ratio, density, Poisson ratio and Lame’s constant. Regarding loose sand reservoir in lower member of Minghuazhen formation, 32-6 south districts in Qinhuangdao,as the research object, be aimed at the different between shallow layer loose sand and deep layer tight sand, first of all, acquire physical property parameters suitable for this kind of heavy oil pool according to experimental study, establishing initial pressure and shear wave relational model; Afterwards, performing prestack elastic wave forward and inversion research, summarizing rules under the guidance of theoretical research and numerical simulation, performing elastic impedance inversion, calculating rock physics attributes; Finally, predicting sand body distribution according to rock physics parameters, and predicting favorable oil area combine well-logging materials and made good results.

孤南-三合村洼陷含油气系统研究

The development petroleum geology has made people from studying and studying and predicting in statically and respectively the pool-forming conditions of an area such as oil source bed, reservoir, overlying formation, migration, trap and preservation, etc. to regarding these conditions as well as roles of generation, reservation and accumulation as an integrated dynamic evolution development system to do study .Meanwhile apply various simulating means to try to predict from quantitative angle. Undoubtedly, the solution of these questions will accumulate exploration process, cut down exploration cost and obtain remarkable economic and social benefits. This paper which take sedimentology ,structural geology and petroleum geology as guides and take petroleum system theory as nucleus and carry out study thinking of beginning with static factor and integration of point and face as well as regarding dynamic state factor as factor and apply study methods of integration of geology, Lab research and numerical modeling proceed integrated dissect and systematic analysis to GuNan-SanHeCun depression. Also apply methods of integration of sequence stratigraphy, biostratigraphy, petrostratigraphy and seismic data to found the time-contour stratigraphic framework and reveal time-space distribution of depositional system and meantime clarify oil-source bed, reservoir and overlying distribution regular patterns. Also use basin analysis means to study precisely the depositional history, packed sequences and evolution. Meanwhile analyze systematically and totally the fracture sequence and fault quality and fault feature, study the structural form, activity JiCi and time-space juxtaposion as well as roles of fault in migration and accumulation of oil and gas of different rank and different quality fault. Simultaneously, utilize seismic, log, analysis testing data and reservoir geology theory to do systematic study and prediction to GuNan-SanHeCun reservoir, study the reservoir types macroscopic distribution and major controlling factors, reservoir rock, filler and porosity structural features as well as distribution of reservoir physical property in 3D space and do comprehensive study and prediction to major controlling and influential factors of reservoir. Furthermore, develop deepingly organic geochemistry comprehensive study, emphasis on two overlaps of oil source rock (ESI, ES3) organic geochemistry features, including types, maturity and spatial variations of organic matter to predict their source potential .Also apply biological marks to proceed oil-to-source correlation ,thereby establish bases for distribution of petroleum system. This study recover the oil generation history of oil source rocks, evaluate source and hydrocarbon discharge potential ,infer pool-forming stages and point out the accumulation direction as well as discover the forming relations of mature oil-source rock and oil reservoir and develop research to study dynamic features of petroleum system. Meanwhile use systematic view, integrate every feature and role of pool forming and the evolution history and pool-forming history, thereby lead people from static conditions such as oil source bed, reservoir, overlying formation, migration, trap and preservation to dynamically analyzing pool-forming process. Also divide GuNan-SanHeCun depression into two second petroleum system, firstly propose to divide second petroleum system according to fluid tress, structural axis and larger faults of cutting depression, and divide lower part of petroleum system into five secondary systems. Meanwhile establish layer analysis and quantitative prediction model of petroleum model, and do quantitative prediction to secondary petroleum system.

三维油藏数值模拟正反演算法研究

The dynamic prediction of complex reservoir development is one of the important research contents of dynamic analysis of oil and gas development. With the increase development of time, the permeabilities and porosities of reservoirs and the permeability of block reservoir at its boundaries are dynamically changing. How to track the dynamic change of permeability and porosity and make certain the permeability of block reservoir at its boundary is an important practical problem. To study developing dynamic prediction of complex reservoir, the key problem of research of dynamic prediction of complex reservoir development is realizing inversion of permeability and porosity. To realize the inversion, first of all, the fast forward and inverse method of 3-dimension reservoir simulation must be studied. Although the inversion has been widely applied to exploration and logging, it has not been applied to3-dimension reservoir simulation. Therefore, the study of fast forward and inverse method of 3-dimension reservoir simulation is a cutting-edge problem, takes on important realistic signification and application value. In this dissertation, 2-dimension and 3-dimension fluid equations in porous media are discretized by finite difference, obtaining finite difference equations to meet the inner boundary conditions by Peaceman's equations, giving successive over relaxation iteration of 3-dimension fluid equations in porous media and the dimensional analysis. Several equation-solving methods are compared in common use, analyzing its convergence and convergence rate. The alternating direction implicit procedure of 2-dimension has been turned into successive over relaxation iteration of alternating direction implicit procedure of 3-dimension fluid equations in porous media, which possesses the virtues of fast computing speed, needing small memory of computer, good adaptability for heterogeneous media and fast convergence rate. The geological model of channel-sandy reservoir has been generated with the help of stochastic simulation technique, whose cross sections of channel-sandy reservoir are parabolic shapes. This method makes the hard data commendably meet, very suit for geological modeling of containing complex boundary surface reservoir. To verify reliability of the method, theoretical solution and numerical solution are compared by simplifying model of 3-dimension fluid equations in porous media, whose results show that the only difference of the two pressure curves is that the numerical solution is lower than theoretical at the wellbore in the same space. It proves that using finite difference to solve fluid equations in porous media is reliable. As numerical examples of 3-dimension heterogeneous reservoir of the single-well and multi-well, the pressure distributions have been computed respectively, which show the pressure distributions there are clearly difference as difference of the permeabilities is greater than one order of magnitude, otherwise there are no clearly difference. As application, the pressure distribution of the channel-sandy reservoir have been computed, which indicates that the space distribution of pressure strongly relies on the direction of permeability, and is sensitive for space distributions of permeability. In this dissertation, the Peaceman's equations have been modified into solving vertical well problem and horizontal well problem simultaneously. In porous media, a 3D layer reservoir in which contain vertical wells and horizontal wells has been calculated with iteration. For channel-sandy reservoir in which there are also vertical wells and horizontal wells, a 3D transient heterogeneous fluid equation has been discretized. As an example, the space distribution of pressure has been calculated with iteration. The results of examples are accord with the fact, which shows the modification of Peaceman's equation is correct. The problem has been solved in the space where there are vertical and horizontal wells. In the dissertation, the nonuniform grid permeability integration equation upscaling method, the nonuniform grid 2D flow rate upscaling method and the nonuniform grid 3D flow rate upscaling method have been studied respectively. In those methods, they enhance computing speed greatly, but the computing speed of 3D flow rate upscaling method is faster than that of 2D flow rate upscaling method, and the precision of 3D flow rate upscaling method is better than that of 2D flow rate upscaling method. The results also show that the solutions of upscaling method are very approximating to that of fine grid blocks. In this paper, 4 methods of fast adaptive nonuniform grid upscaling method of 3D fluid equations in porous media have been put forward, and applied to calculate 3D heterogeneous reservoir and channel-sandy reservoir, whose computing results show that the solutions of nonuniform adaptive upscaling method of 3D heterogeneous fluid equations in porous media are very approximating to that of fine grid blocks in the regions the permeability or porosity being abnormity and very approximating to that of coarsen grid blocks in the other region, however, the computing speed of adaptive upscaling method is 100 times faster than that of fine grid block method. The formula of sensitivity coefficients are derived from initial boundary value problems of fluid equations in porous media by Green's reciprocity principle. The sensitivity coefficients of wellbore pressure to permeability parameters are given by Peaceman's equation and calculated by means of numerical calculation method of 3D transient anisotropic fluid equation in porous media and verified by direct method. The computing results are in excellent agreement with those obtained by the direct method, which shows feasibility of the method. In the dissertation, the calculating examples are also given for 3D reservoir, channel-sandy reservoir and 3D multi-well reservoir, whose numerical results indicate: around the well hole, the value of the sensitivity coefficients of permeability is very large, the value of the sensitivity coefficients of porosity is very large too, but the sensitivity coefficients of porosity is much less than the sensitivity coefficients of permeability, so that the effect of the sensitivity coefficients of permeability for inversion of reservoir parameters is much greater than that of the sensitivity coefficients of porosity. Because computing the sensitivity coefficients needs to call twice the program of reservoir simulation in one iteration, realizing inversion of reservoir parameters must be sustained by the fast forward method. Using the sensitivity coefficients of permeability and porosity, conditioned on observed valley erosion thickness in wells (hard data), the inversion of the permeabilities and porosities in the homogeneous reservoir, homogeneous reservoir only along the certain direction and block reservoir are implemented by Gauss-Newton method or conjugate gradient method respectively. The results of our examples are very approximating to the real data of permeability and porosity, but the convergence rate of conjugate gradient method is much faster than that of Gauss-Newton method.