Estimation and Spatiotemporal Analysis of Methane Emissions from Agriculture in China

National Natural Science Foundation of China [30721140307, 30590380]; Chinese Academy of Sciences (CAS) [KZCX2-YW-432]

Density-dependent effects on seston dynamics and rates of filtering and biodeposition of the suspension-cultured scallop Chlamys farreri in a eutrophic bay (northern China): An experimental study in semi-in situ flow-through systems

Effects of stocking density on seston dynamics and filtering and biodeposition by the suspension-cultured Zhikong scallop Chlamys farreri Jones et Preston in a eutrophic bay (Sishili Bay, northern China), were determined in a 3-month semi-field experiment with continuous flow-through seawater from the bay. Results showed that the presence of the scallops could strongly decrease seston and chlorophyll a concentrations in the water column. Moreover, in a limited water column, increasing scallop density could cause seston depletion due to scallop's filtering and biodeposition process, and impair scallop growth. Both filtration rate and biodeposition rate of C. farreri showed significant negative correlation with their density and positive relationship with seston concentration. Calculation predicts that the daily removal of suspended matter from water column by the scallops in Sishili Bay ecosystem can be as high as 45% of the total suspended matter; and the daily production of biodeposits by the scallops in early summer in farming zone may amount to 7.78 g m(-2), with daily C, N and P biodeposition rates of 3.06 x 10(-1), 3.86 x 10(-2) and 9.80 x 10(-3) g m(-2), respectively. The filtering and biodeposition by suspension-cultured scallops could substantially enhance the deposition of total suspended particulate material, suppress accumulation of particulate organic matter in water column, and increase the flux of C, N and P to benthos, strongly enhancing pelagic-benthic coupling. It was suggested that the filtering-biodeposition process by intensively suspension-cultured bivalve filter-feeders could exert strong top-down control on phytoplankton biomass and other suspended particulate material in coastal ecosystems. This study also indicated that commercially suspension-cultured bivalves may simultaneously and potentially aid in mitigating eutrophication pressures on coastal ecosystems subject to anthropogenic N and P loadings, serving as a eutrophic-environment bioremediator. The ecological services (e.g. filtering capacity, top-down control, and benthic-pelagic coupling) functioned by extractive bivalve aquaculture should be emphasized in coastal ecosystems. (c) 2005 Elsevier B.V. All rights reserved.

A novel three-dimensional imaging array recording and computational method by means of coded cameras reconstruction

In this paper, we propose a novel three-dimensional imaging method by which the object is captured by a coded cameras array (CCA) and computationally reconstructed as a series of longitudinal layered surface images of the object. The distribution of cameras in array, named code pattern, is crucial for reconstructed images fidelity when the correlation decoding is used. We use DIRECT global optimization algorithm to design the code patterns that possess proper imaging property. We have conducted primary experiments to verify and test the performance of the proposed method with a simple discontinuous object and a small-scale CCA including nine cameras. After certain procedures such as capturing, photograph integrating, computational reconstructing and filtering, etc., we obtain reconstructed longitudinal layered surface images of the object with higher signal-to-noise ratio. The results of experiments show that the proposed method is feasible. It is a promising method to be used in fields such as remote sensing, machine vision, etc. (c) 2006 Elsevier GmbH. All rights reserved.

A new frequency offset estimation for OFDM-based systems

A new carrier frequency offset estimation scheme in orthogonal frequency division multiplexing (OFDM) is proposed. The scheme includes coarse frequency offset estimation and fine frequency offset estimation. The coarse frequency offset estimation method we present is a improvement of Zhang's method. The estimation range of the new method is as large as the overall signal-band width. A new fine frequency offset estimation algorithm is also discussed in this paper. The new algorithm has a better performance than the Schmidl's algorithm. The system we use to calculate and simulate is based on the high rate WLAN standard adopted by the IEEE 802.11 stanidardization group. Numerical results are presented to demonstrate the performance of the proposed algorithm.

A new frequency offset estimation for OFDM-Based systems

A-new-carrier-frequency offset estimation scheme in orthogonal frequency division multiplexing (OFDM) is proposed. The scheme includes coarse frequency offset estimation and fine frequency offset estimation. The coarse frequency offset estimation method we present is a improvement of Zhang's method. The estimation range of the new method is as large as the overall signal-band width. A new fine frequency offset estimation algorithm is also discussed in this paper. The new algorithm has a better performance than the Schmidt's algorithm. The system we use to calculate and simulate is based on the high rate WLAN standard adopted by the IEEE 802.11 standardization group. Numerical results are presented to demonstrate the performance of the proposed algorithm.

Determination of sexual hormones in liquid cosmetics by polymer monolith microextraction coupled with high performance liquid chromatography

A method was presented for the determination of testosterone, methyltestosterone and progesterone in liquid cosmetics by coupling polymer monolith microextraction (PMME) to high performance liquid chromatography with UV detection. A poly (methacrylic acid-ethylene glycol dimethacrylate) monolithic capillary column was selected as the extraction medium, which showed high extraction capacity towards these compounds. To achieve optimum extraction performance, several parameters relating to PMME were investigated, including extraction flow rate and pH value, inorganic salt and organic phase concentration of the sample matrix. By simple dilution with phosphate solution and filtering, the sample solution then could be directly injected into the device for extraction. The limits of detection of testosterone, methyltestosterone and progesterone were calculated to be 2, 3, 2, 8 and 4.6 mu g/L. Good linearity was achieved in the range of 10 to 1000 mu g/L with a linear coefficient. r value above 0. 996. Excellent method reproducibility was found by intra- and inter-day precisions, yielding the relative standard deviations of < 7. 7 % and < 7. 5 %, respectively. Recovery for them in the real samples was between 83% and 119%.

吴起县不同退耕阶段农地资源生态服务价值评估分析

以陕西省吴起县为例,对三个不同退耕阶段农地资源生态服务价值进行了评估,分析了各阶段农地资源生态服务价值的变化,探讨了退耕还林(草)工程的实施对农地资源生态服务价值的影响.结果表明:退耕前、退耕初期及退耕中期吴起县农地资源的直接经济产出价值和间接生态服务价值之比分别为:1∶0.84、1∶2.63和1∶3.56.随着退耕工程的不断深入,农地资源总生态服务价值由退耕前的45 402.81万元增加到了退耕初期的81 603.81万元和退耕中期的164 925.46万元,实现了经济效益与生态效益的同步发展.

软件成本估算及其不确定性研究

软件成本估算作为软件项目可行性分析、预算、计划以及控制的基础，是软件工程中的一个重要研究领域。自上世纪60年代以来，尽管软件成本估算一直受到研究者的持续关注，但在现实环境中软件成本估算仍然是软件行业面临的一项难题，仍然有很大的深入研究与改进空间。 现实环境下，软件成本估算方法需要接受不完整且不能完全确定的信息，估算可能的软件开发工作量与开发周期，并衡量估算结果的不确定性与风险。软件成本估算方法还需要满足为使用者创造价值、实施成本低、得到人文与技术上支撑等前提条件，才能在现实环境中被接受。成本估算的应用，也还需要随着项目的演进，与涉众协商、项目计划与项目监控过程密切互动。不能处理好现实环境中的不确定性，不能解决方法实施和应用中所面临的多种关键困难，是已提出的大量成本估算方法和模型难以在现实环境得到广泛使用并发挥影响的重要原因。 本文从确定问题、方法改进、方法应用、工具支撑等多个角度，对软件成本估算进行了较为系统的研究，为解决在现实环境改进成本估算现状所面临的关键困难，以及处理软件成本估算的不确定性这个核心问题，提出了一套包含方法、过程与支撑工具的比较系统、完整的针对实际需要的解决方案。本文研究工作的主要贡献包括： 1）提出软件成本估算的问题模型。设计并实施了中国软件行业的软件成本估算现状调查，探索了我国软件成本估算的现存问题，以及改进成本估算所面临的困难。结合文献综述与产业调查结果，运用“技术接受与使用模型”以及“结果链”等方法，提出软件成本估算的问题模型。涵盖技术、人文、经济与管理多个角度，系统地总结了软件成本估算所面临的问题和潜在的改进。 2）提出集成的成本估算方法。不同于目前估算方法均依赖某种固定估算模型，而是将多种子估算模型作为有用的信息输入，在不同环境下基于历史项目数据自动生成适应该环境的集成的估算模型。 3）针对成本估算在应用环节所面临的关键困难，提出了WikiWinWin软件项目多赢协商方法并开发相应的支持工具。帮助项目涉众正确认识并使用成本估算，促进项目演进过程中涉众协商、成本估算及项目计划与执行过程的有效融合，使成本估算更有效发挥作用。 4）提出了分析框架与相应的综合性方法以处理软件成本估算的不确定性。对软件成本估算的不确定性这个核心问题，进行了系统的分析。使用贝叶斯网络以及Monte Carlo仿真对现有估算模型进行扩展，以处理估算输入的不确定性；使用集成成本估算解决估算模型本身的不确定性；并在成本估算的应用环节，以WikiWinWin方法为核心处理软件成本估算的不确定性。 5）设计开发了软件成本建模与估算支撑工具。结合本文前面提出的成本估算方法形成了工具支撑下的集成的成本建模与估算方法（InCoME方法），在处理估算的不确定性、准确性、稳定性、客观与可重复性、透明性以及建模与估算的自动化支持方面都具有较好的能力，较为全面地满足了企业在现实环境中的需要。

High sensitivity beam phase measurement system for isochronous cyclotron

The beam phase measurement system in the HIRFL is introduced. Based on the double-balanced mixer principle using rf-signal mixing and filtering techniques, a stable and sensitive phase measurement system has been developed. The phase history of the ion beam is detected by using a set of capacitive pick-up probes installed in the cyclotron. The phase information of the measurement is necessary for tuning to obtain a optimized isochronous magnetic field which induces to maximize the beam intensity and to optimize the beam quality. The result of the phase measurement is reliable and the accurancy reaches +/- 0.5 degrees.

基于主动建模的无人直升机增强LQR控制

为了解决无人直升机控制问题,通过把主动建模与LQR(Linear Quadratic Regulator)控制相结合,提出一种能补偿模型差的控制方法。该方法在悬停状态下,采用简化模型设计LQR控制器,并通过UKF(Un-scented-Kalman-Filter)在线估计简化模型与全状态模型的模型差,使用模型差作为补偿项对LQR控制增强。针对实际直升机动力学模型进行仿真,验证了基于UKF的估计和增强LQR控制的有效性。仿真实验结果证明,基于UKF的主动建模技术能够快速估计状态和参数变化,并且增强LQR控制能够使系统适应模型不确定性。

青藏高原及邻区面波频散反演--岩石圈-软流层体系结构研究

The Qinghai-Tibet Plateau lies in the place of the continent-continent collision between Indian and Eurasian plates. Because of their interaction the shallow and deep structures are very complicated. The force system forming the tectonic patterns and driving tectonic movements is effected together by the deep part of the lithosphere and the asthenosphere. It is important to study the 3-D velocity structures, the spheres and layers structures, material properties and states of the lithosphere and the asthenosphere for getting knowledge of their formation and evolution, dynamic process, layers coupling and exchange of material and energy. Based on the Rayleigh wave dispersion theory, we study the 3-D velocity structures, the depths of interfaces and thicknesses of different layers, including the crust, the lithosphere and the asthenosphere, the lithosphere-asthenosphere system in the Qinghai-Tibet Plateau and its adjacent areas. The following tasks include: (1)The digital seismic records of 221 seismic events have been collected, whose magnitudes are larger than 5.0 over the Qinghai-Tibet Plateau and its adjacent areas. These records come from 31 digital seismic stations of GSN , CDSN、NCDSN and part of Indian stations. After making instrument response calibration and filtering, group velocities of fundamental mode of Rayleigh waves are measured using the frequency-time analysis (FTAN) to get the observed dispersions. Furthermore, we strike cluster average for those similar ray paths. Finally, 819 dispersion curves (8-150s) are ready for dispersion inversion. (2)From these dispersion curves, pure dispersion data in 2°×2° cells of the areas (18°N-42°N, 70°E-106°E) are calculated by using function expansion method, proposed by Yanovskaya. The average initial model has been constructed by taking account of global AK135 model along with geodetic, geological, geophysical, receiving function and wide-angle reflection data. Then, initial S-wave velocity structures of the crust and upper mantle in the research areas have been obtained by using linear inversion (SVD) method. (3)Taking the results of the linear inversion as the initial model, we simultaneously invert the S wave velocities and thicknesses by using non-linear inversion (improved Simulated Annealing algorithm). Moreover, during the temperature dropping the variable-scale models are used. Comparing with the linear results, the spheres and layers by the non-linear inversion can be recognized better from the velocity value and offset. (4)The Moho discontinuity and top interface of the asthenosphere are recognized from the velocity value and offset of the layers. The thicknesses of the crust, lithosphere and asthenosphere are gained. These thicknesses are helpful to studying the structural differentia between the Qinghai-Tibet Plateau and its adjacent areas and among geologic units of the plateau. The results of the inversion will provide deep geophysical evidences for studying deep dynamical mechanism and exploring metal mineral resource and oil and gas resources. The following conclusions are reached by the distributions of the S wave velocities and thicknesses of the crust, lithosphere and asthenosphere, combining with previous researches. (1)The crust is very thick in the Qinghai-Tibet Plateau, varying from 60 km to 80 km. The lithospheric thickness in the Qinghai-Tibet Plateau is thinner (130-160 km) than its adjacent areas. Its asthenosphere is relatively thicker, varies from 150 km to 230 km, and the thickest area lies in the western Qiangtang. India located in south of Main Boundary thrust has a thinner crust (32-38 km), a thicker lithosphere of about 190 km and a rather thin asthenosphere of only 60 km. Sichuan and Tarim basins have the crust thickness less than 50km. Their lithospheres are thicker than the Qinghai-Tibet Plateau, and their asthenospheres are thinner. (2)The S-wave velocity variation pattern in the lithosphere-asthenosphere system has band-belted distribution along east-westward. These variations correlate with geology structures sketched by sutures and major faults. These sutures include Main Boundary thrust (MBT), Yarlung-Zangbo River suture (YZS), Bangong Lake-Nujiang suture (BNS), Jinshajiang suture (JSJS), Kunlun edge suture (KL). In the velocity maps of the upper and middle crust, these sutures can be sketched. In velocity maps of 250-300 km depth, MBT, BNS and JSJS can be sketched. In maps of the crustal thickness, the lithospheric thickness and the asthenospheric thickness, these sutures can be still sketched. In particular, MBT can be obviously resolved in these velocity maps and thickness maps. (3)Since the collision between India and Eurasian plate, the “loss” of surface material arising from crustal shortening is caused not only by crustal thickening but also by lateral extrusion material. The source of lateral extrusion lies in the Qiangtang block. These materials extrude along the JSJS and BNS with both rotation and dispersion in Daguaiwan. Finally, it extends toward southeast direction. (4)There is the crust-mantle transition zone of no distinct velocity jump in the lithosphere beneath the Qiangtang Terrane. It has thinner lithosphere and developed thicker asthenosphere. It implies that the crust-mantle transition zone of partial melting is connected with the developed asthenosphere. The underplating of asthenosphere may thin the lithosphere. This buoyancy might be the main mechanism and deep dynamics of the uplift of the Qinghai-Tibet hinterland. At the same time, the transport of hot material with low velocity intrudes into the upper mantle and the lower crust along cracks and faults forming the crust-mantle transition zone.

青藏高原东部电性结构特征及其构造意义

Eastern Himalayan Syntaxis (EHS) and its surroundings (eastern margin of Tibet) is one of the most complicated tectonic areas in the world. As the exhaust opening of the balanced materials of the Tibetan Plateau during the collision of Indan and Eurasian plates, the deep structure beneath EHS surrounding region is referred to as the key to the study of the dynamics of the plateau. EHS3D project, sponsored by NSFC, has been proposed to explore the deep electric features of the area. During the first stage of EHS3D(2006-2008), MT+LMT measurements have been conducted along two lines from Chayu to Qingshuihe (EHS3D-3) and Chayu to Ruoergai (EHS3D-2). This paper will discuss the MT models of EHS3D-3 line. By the data procrssing, including distortion analysis, Robust estimation and strike decomposition, rotated apparent resitivities and phases have been obtained for each station. Then conventional 2-D inversion algorithms (NLCG and RRI) were employed to produce 2-D models. The final preferred 2-D model suggests that the upper crust consists of resistive blocks while in mid-lower crust there are two extensive conductive bodies beneath Lhasa block and Qiangtang terrain respectively. Jinshajiang suture is a gradient belt and Bangong-Nujiang suture appear a conductive belt dipping to the north. . We concluded that the formation of the two conductive bodies attributed to the partial melt and fluids in the lower crust. The regional electric strike derived from decomposition analysis indicates that the crust and upper mantle move in different manners. The upper crust moves like slips of rigid blocks along major slip faults while the lower crust creeps as a flow in the conductive channels.

基于MRF模型的油气储层属性随机建模方法研究

Stochastic reservoir modeling is a technique used in reservoir describing. Through this technique, multiple data sources with different scales can be integrated into the reservoir model and its uncertainty can be conveyed to researchers and supervisors. Stochastic reservoir modeling, for its digital models, its changeable scales, its honoring known information and data and its conveying uncertainty in models, provides a mathematical framework or platform for researchers to integrate multiple data sources and information with different scales into their prediction models. As a fresher method, stochastic reservoir modeling is on the upswing. Based on related works, this paper, starting with Markov property in reservoir, illustrates how to constitute spatial models for catalogued variables and continuum variables by use of Markov random fields. In order to explore reservoir properties, researchers should study the properties of rocks embedded in reservoirs. Apart from methods used in laboratories, geophysical means and subsequent interpretations may be the main sources for information and data used in petroleum exploration and exploitation. How to build a model for flow simulations based on incomplete information is to predict the spatial distributions of different reservoir variables. Considering data source, digital extent and methods, reservoir modeling can be catalogued into four sorts: reservoir sedimentology based method, reservoir seismic prediction, kriging and stochastic reservoir modeling. The application of Markov chain models in the analogue of sedimentary strata is introduced in the third of the paper. The concept of Markov chain model, N-step transition probability matrix, stationary distribution, the estimation of transition probability matrix, the testing of Markov property, 2 means for organizing sections-method based on equal intervals and based on rock facies, embedded Markov matrix, semi-Markov chain model, hidden Markov chain model, etc, are presented in this part. Based on 1-D Markov chain model, conditional 1-D Markov chain model is discussed in the fourth part. By extending 1-D Markov chain model to 2-D, 3-D situations, conditional 2-D, 3-D Markov chain models are presented. This part also discusses the estimation of vertical transition probability, lateral transition probability and the initialization of the top boundary. Corresponding digital models are used to specify, or testify related discussions. The fifth part, based on the fourth part and the application of MRF in image analysis, discusses MRF based method to simulate the spatial distribution of catalogued reservoir variables. In the part, the probability of a special catalogued variable mass, the definition of energy function for catalogued variable mass as a Markov random field, Strauss model, estimation of components in energy function are presented. Corresponding digital models are used to specify, or testify, related discussions. As for the simulation of the spatial distribution of continuum reservoir variables, the sixth part mainly explores 2 methods. The first is pure GMRF based method. Related contents include GMRF model and its neighborhood, parameters estimation, and MCMC iteration method. A digital example illustrates the corresponding method. The second is two-stage models method. Based on the results of catalogued variables distribution simulation, this method, taking GMRF as the prior distribution for continuum variables, taking the relationship between catalogued variables such as rock facies, continuum variables such as porosity, permeability, fluid saturation, can bring a series of stochastic images for the spatial distribution of continuum variables. Integrating multiple data sources into the reservoir model is one of the merits of stochastic reservoir modeling. After discussing how to model spatial distributions of catalogued reservoir variables, continuum reservoir variables, the paper explores how to combine conceptual depositional models, well logs, cores, seismic attributes production history.