基于灰色-马尔科夫模型的我国能源消费结构优化调整

本文以我国1997—2010年的消费结构数据为基础,运用灰色-马尔科夫链模型对数据进行了分析和预测,然后通过MATLAB预测了我国未来十年的能源消费结构,结合我国政府承诺的非化石能源目标对模型预测结果做了进一步的修正,最后根据能源发展趋势给出了能源消费结构的优化策略。     

Designing a mixed performance measurement system for environmental supply chain management using evolutionary game theory and balanced scorecard: A case study of an auto industry supply chain

This paper proposes a mixed performance measurement system using a combination of evolutionary game theory and the balanced scorecard (BSC) in environmental supply chain management (ESCM) that measures and evaluates business operations using the four different perspectives of finance, customer, internal business process, and learning and growth. ESCM plays an important role in the supply chain which leads to the reduction, reuse and recycling of resources involved in both upstream and downstream activities. This paper presents guidance for practical managers in evaluating and measuring ESCM by developing a knowledge-based BSC and evolutionary game theory. The primary purpose of this paper is to apply the proposed method in a case study to one of Iran's biggest auto industry supply chain SAIPA Company. The results of this study indicate that the adoption of ESCM, in the absence of regulatory pressures and cost-saving measures is triggered by public pressures and its implementation is limited by organizational factors and strategic myopia.     

Long-range atmospheric transport of three toxaphene congeners across Europe. Modeling by chained single-box FATEMOD program

Background, aims, and scope  Since toxaphene (polychlorocamphene, polychloropinene, or strobane) mixtures were applied for massive insecticide use in the 1960s to replace the use of DDT, some of their congeners have been found at high latitudes far away from the usage areas. Especially polychlorinated bornanes have demonstrated dominating congeners transported by air up to the Arctic areas. Environmental fate modeling has been applied to monitor this phenomenon using parallel zones of atmosphere around the globe as interconnected environments. These zones, shown in many meteorological maps, however, may not be the best way to configure atmospheric transport in air trajectories. The latter could also be covered by connecting a chain of simple model boxes. We aim to study this alternative approach by modeling the trajectory chain using catchment boxes of our FATEMOD model. Polychlorobornanes analyzed in biota of the Barents Sea offered one case to study this modeling alternative, while toxaphene has been and partly still is used massively at southern East Europe and around rivers flowing to the Aral Sea. Materials and methods  Pure model substances of three polychlorobornanes (toxaphene congeners P26, P50, and P62) were synthesized, their environmentally important thermal properties measured by differential scanning calorimetry, as evaluated from literature data, and their temperature dependences estimated by the QSPR programs VPLEST, WATSOLU, and TDLKOW. The evaluated property parameters were used to model their atmospheric long-range transport from toxaphene heavy usage areas in Ukraine and Aral/SyrDarja/AmuDarja region areas, through East Europe and Northern Norway (Finnmarken) to the Barents Sea. The time period used for the emission model was June 1997. Usual weather conditions in June were applied in the model, which was constructed by chaining FATEMOD model boxes of the catchment’s areas along assumed maximal air flow trajectories. Analysis of the three chlorobornanes in toxaphene mixtures function as a basis for the estimates of emission levels caused by its usage. High estimate (A) was taken from contents in a Western product chlorocamphene and low estimate (B) from mean contents in Russian polychloroterpene products to achieve modeled water concentrations. Bioaccumulation to analyzed lipid of aquatic biota at the target region was estimated by using statistical calculation for persistent organic pollutants in literature. Results  The results from model runs A and B (high and low emission estimate) for levels in sea biota were compared to analysis results of samples taken in August 1997 at Barents Sea. The model results (ng g⁻¹ lw): 4–95 in lipid of planktovores and 7–150 in lipid of piscivores, were in fair agreement with the analysis results from August 1997: 21–31 in Themisto libellula (chatka), 26–42 in Boreocadus saida (Polar cod), and 5–27 in Gadus morhua (cod) liver. Discussion  The modeling results indicate that the application of chained simple multimedia catchment boxes on predicted trajectory is a useful method for estimation of volatile airborne persistent chemical exposures to biota in remote areas. For hazard assessment of these pollutants, their properties, especially temperature dependences, must be estimated by a reasonable accuracy. That can be achieved by using measurements in laboratory with pure model compounds and estimation of properties by thermodynamic QSPR methods. The property parameters can be validated by comparing their values at an environmental temperature range with measured or QSPR-estimated values derived by independent methods. The chained box method used for long-range air transport modeling can be more suitable than global parallel zones modeling used earlier, provided that the main airflow trajectories and properties of transported pollutants are predictable enough. Conclusions  Long-range air transport modeling of persistent, especially photo-resistant organic compounds using a chain of joint simple boxes of catchment’s environments is a feasible method to predict concentrations of pollutants at the target area. This is justified from model results compared with analytical measurements in Barents Sea biota in August 1997: three of six modeled values were high and the other three low compared to the analysis results. The order of magnitude level was similar in both modeled (planktovore and piscivore) and observed (chatka and polar cod) values of lipid samples. The obtained results were too limited to firm validation but are sufficient to justify feasibility of the method, which prompts one to perform more studies on this modeling system. Recommendations and perspectives  For assessment of the risk of environmental damages, chemical fate determination is an essential tool for chemical control, e.g., for EU following the REACH rules. The present conclusion of applicability of the chained single-box multimedia modeling can be validated by further studies using analyses of emissions and target biota in various other cases. To achieve useful results, fate models built with databases having automatic steps for most calculations and outputs accessible to all chemical control professionals are essential. Our FATEMOD program catchments at environments and compound properties listed in the database represent a feasible tool for local, regional, and, according our present test results, for global exposure predictions. As an extended use of model, emission estimates can be achieved by reversed modeling from analysis results of samples corresponding to the target area. This article is dedicated to the memory of Professor Alexander B Terentiev (who passed away in November 2006), our true friend. With his Institute of Organo-Element Compounds, Russian Academy of Science, Moscow, he was an important main organizer of the six joint Finnish–Russian seminars (every third year since 1989) on the field (‘Chemistry and Ecology of Organo-Element Compounds’). He prompted us especially to search properties and environmental fates for various polyhalogen compounds. We remember him for his friendly character and great sense of humor.     

中国旅游业绿色发展效率时空演变特征及影响机理

绿色发展是旅游业可持续发展理念的重要组成部分,是旅游业奉行以人为本、生态至上和全面发展的新价值观。在梳理旅游业绿色发展概念及内涵基础上,构建旅游业绿色发展效率评价体系,运用SBM-Undersirable模型、核密度估计、空间马尔科夫链等方法,探讨2008—2018年中国31个省（市、自治区）旅游业绿色发展效率（TGDE）时空演化特征及影响机理。研究发现：（1）时间和空间变化方面,TGDE总体处于中等偏下水平,时间上呈“W”型变化形态,“下降—上升—调整”阶段特征显著;空间呈“东—中—西”递减分布,内部差异为西部地区>东部地区>中部地区,低、中、高效率由“金字塔”向“菱形”结构转变,高效率地区集中于东部沿海,中等效率多分布于中西部地区,低效率位于胡焕庸线两侧。（2）动态演进方面,TGDE始终存在两极分化现象,但区域协调性逐步增强,具有较强平稳性,难以实现跨越式发展,空间向上转移省份比较集中,以中西部为主,向下调整省份较少,且存在明显的空间溢出效应,溢出影响具有不对称性。（3）影响机理方面,总体上,经济水平、产业结构、政府规制、教育水平和旅游资源影响因子与TGDE间存在显著的正向关系,对外开放程度的作用不显著,但各因子的影响程度、作用机理及条件具有较强地域性。     

基于贝叶斯的苯并(α)芘PBPK模型优化与健康风险评估应用

应用基于生理的药代动力学(PBPK)模型预测苯并(α)芘(BaP)暴露的人体内部剂量,基于贝叶斯的马尔科夫链蒙特卡洛模拟(MCMC)方法对模型参数进行校准和优化,最后运用已优化的模型对BaP内暴露基准值进行推导.研究发现,基于贝叶斯的MCMC方法对模型后验参数校准后,模型精度明显提高,两个数据集验证结果显示残差平方和分别降低了72%和94%.PBPK模型以BaP和子代谢物3-羟基苯并(α)芘(3-OHBaP)的体内动力学过程为结构基础,模拟BaP体内浓度分布大小为脂肪>肾脏>皮肤>缓慢灌注组织>快速灌注组织>静脉血>肝脏;3-OHBaP体内浓度分布大小为肾脏>快速灌注组织>脂肪>肺>静脉血>缓慢灌注组织>肝脏>皮肤.敏感性分析显示,快速灌注组织-血分配系数对模型输出影响最大,灵敏度系数超过了200%;排泄系数影响最小,只有肾小球过滤率K_BR的灵敏度系数超过了1%.以美国国家环境保护局推荐的参考浓度2.0×10^-6mg/m³为外暴露安全基准值,基于PBPK模型推导了职业暴露的BaP生物监测当量(BE),结果显示BE值为0.405pmol/mol肌酐(尿液3-OHBaP平均浓度),为基于人体内暴露剂量水平进行定量健康风险评估奠定了基础.     

Planar detonation structure for chain-branching kinetics with large activation energy and small initiation rate

The structure of the steady planar detonation wave is analyzed for three-step chain-branching kinetics consistent with hydrogen–air chemistry. The initiation and chain-branching steps are described by an Arrhenius rate. Both are thermally neutral, so that heat release is due to termination. The initiation rate is typically very low and very stiff. As a result, a small fraction of the reactant is consumed in the initiation region, which is very long but ends in an exponential chain-branching explosion. Next, the reactant is rapidly converted into chain-branching radicals, within a very thin zone, which ends when the concentration of the chain-branching radical reaches a peak, because of reactant depletion. Finally, the termination step consumes the chain-branching radicals and releases heat within a region thicker than the peak zone, but much thinner than the initiation region. The analysis is based upon two assumptions: that the chain-branching activation energy is high and that initiation is slow. The structure of the initiation and chain-branching zones is different for post-shock states within or outside the explosion region in the chain-branching diagram. In the former situation, chain-branching is already stronger than termination at the von Neumann point, and vice versa. In the no-explosion case, the initiation zone becomes very long, while the little chain-branching specie produced by initiation is directly converted into product by the termination step. Temperature increases slowly until reaching the explosion curve, when chain-branching becomes stronger than termination. The subsequent structure is similar to the explosion case.     

Spatial modelling of the probability of obtaining a detectable concentration of ammonia in Moreton Bay, Australia

Ammonia is an important water quality variable, which in excess, can be detrimental to waterways and their ecosystems. In the Ecosystem Health Monitoring Program in South-east Queensland ammonia is monitored monthly, however, often more than 50% of the ammonia observations in Moreton Bay are below detection limit, making it difficult to draw useful inferences. In this paper a clipped Gaussian random field is used to spatially model and map the probability of detectable concentration of ammonia. The methodology is applied to the Moreton Bay samples collected in February 2005. The results suggest that for this month the oceanic impacted areas have higher probability of detectable ammonia concentration than the areas closer to the main sources of anthropogenic inputs.     

A Bayesian approach to multi-source forest area estimation

In efforts such as land use change monitoring, carbon budgeting, and forecasting ecological conditions and timber supply, there is increasing demand for regional and national data layers depicting forest cover. These data layers must permit small area estimates of forest area and, most importantly, provide associated error estimates. This paper presents a model-based approach for coupling mid-resolution satellite imagery with plot-based forest inventory data to produce estimates of probability of forest and associated error at the pixel-level. The proposed Bayesian hierarchical model provides access to each pixel’s posterior predictive distribution allowing for a highly flexible analysis of pixel and multi-pixel areas of interest. The paper presents a trial using multiple dates of Landsat imagery and USDA Forest Service Forest Inventory and Analysis plot data. The results describe the spatial dependence structure within the trial site, provide pixel and multi-pixel summaries of probability of forest land use, and explore discretization schemes of the posterior predictive distributions to forest and non-forest classes. Model prediction results of a holdout set analysis suggest the proposed model provides high classification accuracy, 88%, for the trial site.

Proposing an information criterion for individual-based models developed in a pattern-oriented modelling framework

Individual-based models (IBMs) have been improved in quality and reliability in recent years with an approach called pattern-oriented modelling (POM). POM proposes guidelines to develop models reproducing multiple patterns observed on the field and to test systematically how well the IBMs reproduce them. POM studies used generally traditional methods of goodness of fit such as the sum of squares evaluation or ad hoc comparisons of fitting errors and variations. Model selection, however, can be a rigorous statistical approach based on information theory and information criteria such as the Akaike's information criterion (AIC) or the deviance information criterion (DIC). So far, it has not been tried to link POM to these rigorous techniques. The main problems to achieve that are: (a) the difficulty to have likelihood functions for IBMs’ parameters and (b) the possibility to obtain posterior distributions of IBMs’ parameters given the patterns to reproduce. In a first part, this paper answers problem (a) by proposing and explaining how to calculate a deviance measure (POMDEV) for models developed in a context of POM. And while answering the second problem, a second part of the paper proposes an information criterion for model selection in a POM context (the pattern-oriented modelling information criterion: POMIC). This criterion does not yet have the same theoretical foundation as, e.g., AIC, but uses formal analogies to the DIC. In a third part POMIC is tested with a modelling exercise. This exercise shows the potential of POMIC to use multiple patterns for selecting among multiple potential submodels and eventually select the most parsimonious and well fitting model version. We conclude that POMIC, although being a heuristically derived approach, can greatly improve the POM framework.