土壤侵蚀高维指数非线性模型及其参数辨识

模型验证和数据库组建基础上,用WinEPIC模型定量模拟研究了黄土高原半湿润区长武、半干旱区固原和半干旱偏旱区海原20～30年内苜蓿草地水分生产潜力、10m土层土壤有效含水量和土壤湿度剖面分布特征的动态变化.结果表明:长武、固原和海原苜蓿草地水分生产潜力模拟值随降水量变化而呈现波动性降低趋势,其平均值分别为8.81、3.83和2.48t.hm-2;长武、固原和海原苜蓿草地10m土层逐月土壤有效含水量模拟值均呈现明显的波动性降低趋势,模拟初期,4～8年生苜蓿草地土壤干燥化趋势十分强烈,此后,随降水量变化长期在较低水平上波动;随着苜蓿生长年限的延长,苜蓿草地土壤干层逐年加深、加厚,长武、固原和海原土壤干层分布深度达到10m所需时间依次为6、6和4年,此后苜蓿草地降水渗深以下土层长期维持较为稳定的干燥化状态;苜蓿草地水分持续利用的合理年限为半湿润区8～10年,半干旱区6～8年,半干旱偏旱区4～6年.

油藏动态模型和剩余油分布模式

This paper belong to national "973" technological project undertaken by Shengli Oilfield. Work area is composed of turbidite reservoir of S212 and delta reservoir of S283 of Sheng2 unit in Shengtuo Oilfield that has a 36 years water injection development history. Change of the macroscopic, microscopic and filterational parameters and its mechanism have been studied in the 4 water-cut stages i.e. the primary, moderate, high and supper-high stage by using multi-disciplinary theories and approaches, computer techniques and physical simulation comprehensively. Dynamic reservoir models to different water-cut stages have been established. The study of forming mechanism and distribution of residual oil revealed the main types and spatial distribution of residual oil in different water-cut stages and the distribution mode has also been built up. Macroscopic, microscopic and filterational parameters selecting principle, optimizing and selecting standard, matching standard and laws and related database of various dynamic parameters in different water-cut stages have been established, which laid good basis for revealing reservoir macroscopic, microscopic and filterational parameters' dynamic change and residual oil distribution. The study indicated that in general, the macroscopic, microscopic and filterational parameters will slowly increase and become better in both shallow turbidite and delta reservoirs with the increasing of water cut, but different reservoirs have their own characteristics and change laws. Parameters of I~2 unit, whose petrophysical properties are better, increase more quickly than 8~3, whose petrophysical properties are more unfavorable. The changes was relatively quickly in high water-cut stage, while relatively slowly from primary to moderate and from high to supper-high water-cut stage. This paper firstly put forward that reservoir macroscopic, microscopic and filterational parameters are controlled by dynamic geological function of reservoir fluid, which is considered the major reason of reservoir parameters' dynamic changes and residual oil formation and distribution during reservoir development. Physical simulation of filterational parameters verified that forming mechanism and distribution of residual oil in different water-cut stages are also controlled by dynamic geological function of reservoir fluid. The idea of fluid geological function during reservoir development developed the theory of development geology, and has important practical values. This paper firstly constructed dynamic geological and mathematical models and five modes of residual oil distribution in Shengtuo Oilfield, and achieved four-dimensional forecast of residual oil distribution in different watercut stages. Dynamic changes and mechanism of macroscopic, microscopic and fliterational parameters of reservoir and their change process have been revealed. Forecast of residual oil distribution has been achieved by computers. This paper established the related theories, approaches and techniques for residual oil study, characterization and in different water-cut stages, and realized dynamic forecast of residual oil. It gained remarkable economic benefit and social effect in guiding field development. These theories and techniques had important meaningfulness for residual oil prediction in the terrestrial faulted basins not only in Shengli Oilfield but also in the east of China. Furthermore, this study has developed the theory of development geology.

未来气候变化对黄土高原黑河流域水资源的影响

气候变化对黄土高原的水资源有重要影响,对其影响进行评估可以为区域发展提供重要的决策依据。基于分布式水文模型SWAT和4种全球环流模式的各3种排放情景,评估了2010～2039年黄土高塬沟壑区黑河流域水资源对气候变化的潜在响应。结果表明,黑河流域2010～2039年的年均降水变化-2.3%～7.8%,年均最高和最低温度分别升高0.7～2.2℃和1.2～2.8℃,年均径流量变化-19.8%～37.0%,1.2m剖面年均土壤水分含量变化-5.5%～17.2%,年均蒸散量普遍增长0.1%～5.9%;水文气象变量变化趋势复杂,但T检验表明年降水、径流、土壤水分和蒸散增长的概率较大。对于季节变化,降水可能在12～7月份和9月份增长,8月份和10～11月份减少;径流在4～7月份和9～10月份增加,11～3月份和8月份减少;土壤水分在各月都增长;蒸散11～6月份普遍增长,7～10月份减少的可能性较大。未来气候将发生显著变化并对水资源有重要影响,需采取必要的措施来减缓其不利影响。

A global shallow water model using high order multi-moment constrained finite volume method and icosahedral grid

A novel accurate numerical model for shallow water equations on sphere have been developed by implementing the high order multi-moment constrained finite volume (MCV) method on the icosahedral geodesic grid. High order reconstructions are conducted cell-wisely by making use of the point values as the unknowns distributed within each triangular cell element. The time evolution equations to update the unknowns are derived from a set of constrained conditions for two types of moments, i.e. the point values on the cell boundary edges and the cell-integrated average. The numerical conservation is rigorously guaranteed. in the present model, all unknowns or computational variables are point values and no numerical quadrature is involved, which particularly benefits the computational accuracy and efficiency in handling the spherical geometry, such as coordinate transformation and curved surface. Numerical formulations of third and fourth order accuracy are presented in detail. The proposed numerical model has been validated by widely used benchmark tests and competitive results are obtained. The present numerical framework provides a promising and practical base for further development of atmospheric and oceanic general circulation models. (C) 2009 Elsevier Inc. All rights reserved.

Intraseasonal variability of Indian Ocean sea surface temperature during boreal winter: Madden-Julian Oscillation versus submonthly forcing and processes

[ 1] Intraseasonal variability of Indian Ocean sea surface temperature (SST) during boreal winter is investigated by analyzing available data and a suite of solutions to an ocean general circulation model for 1998 - 2004. This period covers the QuikSCAT and Tropical Rainfall Measuring Mission (TRMM) observations. Impacts of the 30 - 90 day and 10 - 30 day atmospheric intraseasonal oscillations (ISOs) are examined separately, with the former dominated by the Madden-Julian Oscillation (MJO) and the latter dominated by convectively coupled Rossby and Kelvin waves. The maximum variation of intraseasonal SST occurs at 10 degrees S - 2 degrees S in the wintertime Intertropical Convergence Zone (ITCZ), where the mixed layer is thin and intraseasonal wind speed reaches its maximum. The observed maximum warming ( cooling) averaged over ( 60 degrees E - 85 degrees E, 10 degrees S - 3 degrees S) is 1.13 degrees C ( - 0.97 degrees C) for the period of interest, with a standard deviation of 0.39 degrees C in winter. This SST change is forced predominantly by the MJO. While the MJO causes a basin-wide cooling ( warming) in the ITCZ region, submonthly ISOs cause a more complex SST structure that propagates southwestward in the western-central basin and southeastward in the eastern ocean. On both the MJO and submonthly timescales, winds are the deterministic factor for the SST variability. Short-wave radiation generally plays a secondary role, and effects of precipitation are negligible. The dominant role of winds results roughly equally from wind speed and stress forcing. Wind speed affects SST by altering turbulent heat fluxes and entrainment cooling. Wind stress affects SST via several local and remote oceanic processes.

太阳辐射对热层和电离层变化性的影响

Solar ultraviolet (UV) radiation at wavelengths less than 400 nm is an important source of energy for aeronomic processes throughout the solar system. Solar UV photons are absorbed in planetary atmospheres, as well as throughout the heliosphere, via photodissociation of molecules, photoionization of molecules and atoms, and photoexcitation toexcitation including resonance scattering. In this paper, the solar irradiances data measured by TIMED SEE, as well as the solar proxies such as F10.7 and Mg II, thermosphere neutral density of CHAMP measurements and topside ionospheric plasmas densities from DMSP, are used to analyze solar irradiance effects on the variabilities of the thermosphere and the ionosphere. First, thermosphere densities near 410 km altitude are analyzed for solar irradiance variability effects during the period 2002-2004. Correlations between the densities and the solar irradiances for different spectral lines and wavelength ranges reveal significantly different characteristics. The density correlates remarkably well with all the selected solar irradiances except the lower chromospheric O I (130.4 nm) emission. Among the chosen solar proxies, the Mg II core-to-wing ratio index, EUV (30-120 nm) and F10.7 show the highest correlations with the density for short-term (< ~27 days) variations. For both long- (> ~27 days) and short-term variations, linear correlation coefficients exhibit a decreasing trend from low latitudes towards high latitudes. The density variability can be effectively modeled (capturing 71% of the variance) using multiple solar irradiance indices, including F10.7, SEUV (the EUV 30-120 nm index), and SFUV (the FUV 120-193 nm index), in which a lag time of 1 day was used for both F10.7 and SEUV, and 5 days for SFUV. In our regression formulation SEUV has the largest contribution to the density variation (40%), with the F10.7 having the next largest contribution (32%) and SFUV accounting for the rest (28%). Furthermore, a pronounced period of about 27.2 days (mean period of the Sun's rotation) is present in both density and solar irradiance data of 2003 and 2004, and a pronounced period of about 54.4 days (doubled period of the solar rotation) is also revealed in 2004. However, soft X-ray and FUV irradiances did not present a pronounced 54.4 day period in 2004, in spite of their high correlation with the densities. The Ap index also shows 54-day periodicities in 2004, and magnetic activity, together with solar irradiance, affects the 54-day variation in density significantly. In addition, NRLMSISE00, DTM-2000 and JB2006 model predictions are compared with density measurements from CHAMP to assess their accuracy, and the results show that these models underestimate the response of the thermosphere to variations induced by solar rotation. Next, the equatorial topside ionospheric plasmas densities Ni are analyzed for solar irradiance variability effects during the period 2002-2005. Linear correlations between Ni and the solar irradiances for different wavelength ranges reveal significantly different characteristics. XUV (0-35 nm) and EUV (115-130 nm) show higher correlation with Ni for the long-term variations, whereas EUV (35-115 nm) show higher correlation for the short-term variations. Moreover, partial correlation analysis shows that the long-term variations of Ni are affected by both XUV (0-35 nm) and EUV (35-115 nm), whereas XUV (0-35 nm) play a more important role; the short-term variations of Ni are mostly affected by EUV (35-115 nm). Furthermore, a pronounced period of about 27 days is present in both Ni and solar irradiance data of 2003 and 2004, and a pronounced period of about 54 days is also revealed in 2004. Finally, prompted by previous studies that have suggested solar EUV radiation as a means of driving the semiannual variation, we investigate the intra-annual variation in thermosphere neutral density near 400 km during 2002-2005. The intra-annual variation, commonly referred to as the ‘semiannual variation’, is characterized by significant latitude structure, hemispheric asymmetries, and inter-annual variability. The magnitude of the maximum yearly difference, from the yearly minimum to the yearly maximum, varies by as much as 60% from year to year, and the phases of the minima and maxima also change by 20-40 days from year to year. Each annual harmonic of the intra-annual variation, namely, annual, semiannual, ter-annual and quatra-annual, exhibits a decreasing trend from 2002 through 2005 that is correlated with the decline in solar activity. In addition, some variations in these harmonics are correlated with geomagnetic activity, as represented by the daily mean value of Kp. Recent empirical models of the thermosphere are found to be deficient in capturing most of the latitude dependencies discovered in our data. In addition, the solar flux and geomagnetic activity proxies that we have employed do not capture some latitude and inter-annual variations detected in our data. It is possible that these variations are partly due to other effects, such as seasonal-latitudinal variations in turbopause altitude (and hence O/N2 composition) and ionosphere coupling processes that remain to be discovered in the context of influencing the intra-annual variations depicted here. Our results provide a new dataset to challenge and validate thermosphere-ionosphere general circulation models that seek to delineate the thermosphere intra-annual variation and to understand the various competing mechanisms that may contribute to its existence and variability. We furthermore suggest that the term “intra-annual” variation be adopted to describe the variability in thermosphere and ionosphere parameters that is well-captured through a superposition of annual, semiannual, ter-annual, and quatra-annual harmonic terms, and that “semiannual’ be used strictly in reference to a pure 6-monthly sinusoidal variation. Moreover, we propose the term “intra-seasonal” to refer to those shorter-term variations that arise as residuals from the above Fourier representation.

Comparison of the Dynamic Response of a Saturated Soil Between Two Models

By comparing the dynamic responses of saturated soil to Biot's and Yamamoto's models, the properties of the two models have be pointed out. First of all, an analysis has been made for energy loss of each model from the basic equations. Then the damping of elastic waves in coarse sand and fine sand with loading frequency and soil's parameters have been calculated and the representation of viscous friction and Coulomb friction in the two models has been concluded. Finally, the variations of loading wave damping and stress phase angles with water depth and soil's parameters have been obtained as loading waves range in ocean waves.

Statistical and stochastic description of microvoid evolution observed in dynamic failure of ductile materials

A brief review is presented of statistical approaches on microdamage evolution. An experimental study of statistical microdamage evolution in two ductile materials under dynamic loading is carried out. The observation indicates that there are large differences in size and distribution of microvoids between these two materials. With this phenomenon in mind, kinetic equations governing the nucleation and growth of microvoids in nonlinear rate-dependent materials are combined with the balance law of void number to establish statistical differential equations that describe the evolution of microvoids' number density. The theoretical solution provides a reasonable explanation of the experimentally observed phenomenon. The effects of stochastic fluctuation which is influenced by the inhomogeneous microscopic structure of materials are subsequently examined (i.e. stochastic growth model). Based on the stochastic differential equation, a Fokker-Planck equation which governs the evolution of the transition probability is derived. The analytical solution for the transition probability is then obtained and the effects of stochastic fluctuation is discussed. The statistical and stochastic analyses may provide effective approaches to reveal the physics of damage evolution and dynamic failure process in ductile materials.

A second-order dynamic subgrid-scale stress model

A second-order dynamic model based on the general relation between the subgrid-scale stress and the velocity gradient tensors was proposed. A priori test of the second-order model was made using moderate resolution direct numerical simulation date at high Reynolds number ( Taylor microscale Reynolds number R-lambda = 102 similar to 216) for homogeneous, isotropic forced flow, decaying flow, and homogeneous rotating flow. Numerical testing shows that the second-order dynamic model significantly improves the correlation coefficient when compared to the first-order dynamic models.

Chemical equilibrium analysis and CFD simulation of coal combustion and NOx formation

A full two-fluid model of reacting gas-particle flows with an algebraic unified second-order moment (AUSM) turbulence-chemistry model is used to simulate Beijing coal combustion and NOx formation. The sub-models are the k-epsilon-kp two-phase turbulence model, the EBU-Arrhenius volatile and CO combustion model, the six-flux radiation model, coal devolatilization model and char combustion model. The blocking effect on NOx formation is discussed. In addition, the chemical equilibrium analysis is used to predict NOx concentration at different temperature. Results of CID simulation and chemical equilibrium analysis show that, optimizing air dynamic parameters can delay the NOx formation and decrease NOx emission, but it is effective only in a restricted range. In order to decrease NOx emission near to zero, the re-burning or other chemical methods must be used.

Comparisons of static, quasi-static and dynamic 3D porous media scale network models for two-phase immiscible flow in porous media

A dynamic 3D pore-scale network model is formulated for investigating the effect of interfacial tension and oil-water viscosity during chemical flooding. The model takes into account both viscous and capillary forces in analyzing the impact of chemical properties on flow behavior or displacement configuration, while the static model with conventional invasion percolation algorithm incorporates the capillary pressure only. From comparisons of simulation results from these models. it indicates that the static pore scale network model can be used successfully when the capillary number is low. With the capillary increases due to the enhancement of water viscosity or decrease of interfacial tension, only the quasi-static and dynamic model can give insight into the displacement mechanisms.

A Dynamic Model For Ice-Induced Vibration Of Structures

A dynamic model for the ice-induced vibration (IIV) of structures is developed in the present study. Ice properties have been taken into account, such as the discrete failure, the dependence of the crushing strength on the ice velocity, and the randomness of ice failure. The most important prediction of the model is to capture the resonant frequency lock-in, which is analog to that in the vortex-induced vibration. Based on the model, the mechanism of resonant IIV is discussed. It is found that the dependence of the ice crushing strength on the ice velocity plays an important role in the resonant frequency lock-in of IIV. In addition, an intermittent stochastic resonant vibration is simulated from the model. These predictions are supported by the laboratory and field observations reported. The present model is more productive than the previous models of IIV.

Discussion on the basic mathematical models of neurons - Double synaptic weight neuron in general-purpose neurocomputer, theory and application

Based on the introduction of the traditional mathematical models of neurons in general-purpose neurocomputer, a novel all-purpose mathematical model-Double synaptic weight neuron (DSWN) is presented, which can simulate all kinds of neuron architectures, including Radial-Basis-Function (RBF) and Back-propagation (BP) models, etc. At the same time, this new model is realized using hardware and implemented in the new CASSANN-II neurocomputer that can be used to form various types of neural networks with multiple mathematical models of neurons. In this paper, the flexibility of the new model has also been described in constructing neural networks and based on the theory of Biomimetic pattern recognition (BPR) and high-dimensional space covering, a recognition system of omni directionally oriented rigid objects on the horizontal surface and a face recognition system had been implemented on CASSANN-H neurocomputer. The result showed DSWN neural network has great potential in pattern recognition.