Formalization and specification for controlling process modularity in mobile computing

新的计算模式，普适计算和全局计算，正在作为高度分布式和移动计算的计算模式展现出来。这篇论文探讨了在抽象层面上支持这些新型计算模式的适合的形式化基础，关注在进程移动单位上的控制, 以便在分布式与移动计算环境下更好地协调进程的移动性。 论文的第一部分概述了针对分布式、移动计算的现有进程演算模型中的进程移动单元，并且设计了一种在此方面更优、更具弹性的进程框架。为了表示这种进程框架，我们提出了一种新的、针对移动和分布式系统的进程演算,这种进程演算的优点是动态、弹性的控制进程的移动单元；具体的思路就是扩展π- calculus以及其支持分布式和移动性的变体。我们把这种新的演算叫做Modular π-calculus。我们通过这种演算的提出来说明进程框架提供了一种针对移动进程更为合适的协调机制以及编程模型，例如移动的代理和动态组件载入的支持。之后，我们通过讨论互模拟的几种提法来具体说明能够反映演算设计的进程描述的关键，之后我们讨论了它们的具体性质。 本文的第二部分提出了一个对进程模型的行为和性质进行推理的规约框架。首先，提出了一个对Modularπ-calculus中进程的系统性质进行规约的模态逻辑。为了更好的理解该逻辑，文中对由这个逻辑推出的进程等价的特征进行了研究，并且证明了该逻辑的区分能力介于互模拟和结构一致之间。接下来关于这个规约框架的自动化，本文针对该逻辑和Modular π-calculus的有限控制子集，提出了模型检测算法，并且给出了算法正确性的证明。同时文中贯穿了一些实际且直观的例子，以展现本文提出的一组框架即演算、逻辑和模型算法的有效性。

BPEL引擎监控平台的设计与实现

WS-BPEL作为WEB服务复合事实上的标准正引起越来越多的关注，它可以用来整合物理上分布的服务，构建复杂的业务流程。复合服务的执行依赖于分布的自治的第三方服务，其实际行为可能与预先指定的行为约束不一致；此外，由于复合服务的动态性和灵活性，其本身在运行期间可能发生动态调整，以适应不断变化的业务环境，因此需要对复合服务的行为在运行期进行连续的监控，保证复合服务的正确执行。同时，通过对复合服务的运行时监控，可以获取执行期间的各种性能信息，从而为流程的优化和动态的服务选择提供依据。 本文首先对业务流程运行时的监控需求进行了分析，然后针对这些监控需求，结合BPEL流程的特点，设计了一种新型的基于事件的监控需求描述语言Mon-BPEL，它可以精确的刻画功能性和非功能性的监控需求。相对于其他基于事件的监控需求描述语言(如 Event Calculus，Past LTL等)，Mon-BPEL可以方便的刻画流程执行事件序列之间的时序关系，同时，它具有简洁易懂，用户友好等特点。 此外，本文研究了基于AOP技术的业务流程监控机制。论文设计了一种针对BPEL语言的AOP扩展方案，以支持对业务流程运行期的动态监控。通过将监控需求映射为方面，从而分离了关注点，降低了业务流程开发的复杂度，提高了流程的可维护性。接着给出了Mon-BPEL监控需求表达式的的自动机语义，并基于自动机理论对监控需求表达式进行计算，检查指定的监控需求是否违背。最后，论文探讨了Aspect的动态编织机制，将监控功能引入BPEL引擎。 论文最后给出了监控平台在OnceBPEL引擎中的实现，并将前面提到的研究成果引入其中,通过实验对系统的功能及性能进行了测试与评估。

Estimation of Transverse Thermal Conductivity of Doubly-periodic Fiber Reinforced Composites

For steady-state heat conduction a new variational functional for a unit cell of composites with periodic microstructures is constructed by considering the quasi-periodicity of the temperature field and in the periodicity of the heat flux fields. Then by combining with the eigenfunction expansion of complex potential which satisfies the fiber-matrix interface conditions, an eigenfunction expansion-variational method (EEVM) based on a unit cell is developed. The effective transverse thermal conductivities of doubly-periodic fiber reinforced composites are calculated, and the first-order approximation formula for the square and hexagonal arrays is presented,which is convenient for engineering application. The numerical results show a good convergency of the presented method, even through the fiber volume fraction is relatively high. Comparisons with the existing analytical and experimental results are made to demonstrate the accuracy and validity of the first-order approximation formula for the hexagonal array.

钻井废水可生化性研究及处理技术探索

钻井废水是油气井开采钻探过程中产生的废水，钻井废水成分复杂，有机物浓度高、色度高、悬浮物浓度高，水质变化大，排放点分散，不经处理排放会污染环境，破坏生态。随着石油工业的不断发展和国家环保法律法规的日益严格，钻井废水的治理也越来越受到重视。如何采用经济有效的方法处理废弃钻井液，对油气井开采业的可持续发展具有重要意义。本论文以遂宁磨153 井的钻井废水为主要研究对象，在对废水进行絮凝沉降预处理和生物法处理探索的基础上，针对钻井废水可生化性差的特点，采用水解酸化和Fenton 试剂改善钻井废水的可生化性，对反应过程进行了比较详细的考察，对可生化性改善的机理进行了探索。主要研究结论如下：1 用PFS 和PAC 配制的混合混凝剂对钻井废水COD 的去除效果比较显著，在最佳条件下COD 的去除率可达75%，且絮体沉降速度较快，出水pH 保持中性；2 水解酸化法处理钻井废水可显著改善废水的可生化性。经48 小时水解酸化处理，钻井废水的理论BOD5可提高约22 倍，表观BOD5/COD值由0.004 提高到0.034。用接触氧化反应器处理经水解酸化处理后的废水，处理效果比较稳定，COD平均去除率达35.5%；3 研究了Fenton反应中各影响因子对废水COD去除率、BOD5/COD的影响并分析其作用机制，确定了最佳条件：初始pH为4.0，H2O2/Fe2+(摩尔浓度比)为20，H2O2/COD（质量浓度比）为1，反应时间为2 个小时。此条件下，废水的COD去除率约为40%，BOD5/COD值从0.002～0.003 提高至0.15～0.2，可生化性得到很大提高。本论文的主要创新点在于：1 以成分复杂、水质变化大的气井钻井废水为研究对象，从理论BOD 和表观BOD 两方面对水解酸化过程中废水可生化性的变化进行了分析；2 对Fenton 试剂改善钻井废水可生化性的过程、主要影响因素进行了比较详细的考察。本论文的研究成果，可为生物法处理钻井废水的深入研究提供理论依据。Drilling wastewater is produced in the process of oil-gas well drilling,because of its complicated composition, high concentrate of organic compound andsuspended solid, high chroma, levity of water quality and decentralization ofdischarge point, it pollutes environment seriously if discharged without treatment.With the development of petroleum industry and the issuing of more strict laws forenvironmental protection, it has been paid more and more attention on drillingwastewater treatment. It is of great importance for the sustainable development ofoil-gas well drilling to treat drilling wastewater by economical and effective methods.In this paper, drilling wastewater of Mo No.153 well in Suining was studied asthe main object. On the basis of research on pre-treatment with flocculant andbiological treatment, and according to the character of poor biodegradability, thedrilling wastewater was treated by hydrolytic acidification and Fenton’s reagent toimprove its biodegradability. The process and mechanism of biodegradabilitychanging were investigated. The primary conclusions are:1 It is effective to treat drilling wastewater with mixing PFS and PAC asflocculant. The removal rates of COD came up to 75% under optimal conditions, thesedimentation rate of flocculation is rapid, and the pH value of treated water remainedneutral;2 The biodegradability of drilling wastewater was highly improved afterhydrolytic acidification process. The theoretic BOD5 of drilling wastewater increasedby 22 times and its detected BOD5/COD ratio increased from 0.004 to 0.034 afterhydrolytic acidification for 48 hours. The wastewater after hydrolytic acidificationwas treated by biological contact oxidation reactor. Stable treatment performance was achieved, and the average removal rates of COD came up to 35.5%;3 The effects of various affection factors on the removal efficiency of COD andBOD5/COD radio in treating drilling wastewater by Fenton’s reagent wereinvestigated and the mechanism was analyzed. The optimal conditions were: initialpH of solution was 4.0, the molar ratio of H2O2 and Fe2+ was 20, the concentrationratio of H2O2 and COD was 1 and the reaction time was 120 min. Under the aboveconditions, the removal efficiency was about 40% and the ratio of BOD5 and CODincreased from 0.002 ¡« 0.003 to 0.15 ¡« 0.2. The biodegradability of drillingwastewater was greatly improved.The innovations of this thesis are:1 The drilling wastewater was taken as the research object which hascomplicated composition and variational water quality, and the changes ofbiodegradability were analyzed from theoretic BOD and detected BOD aspects duringhydrolytic acidification process;2 The biodegradability changing process and primary affection factors of drillingwastewater treating by Fenton’s reagent were investigated.The results of this study could provide theoretic foundation for further researchon biological treatment of drilling wastewater.

Solitary magnon excitations in a one-dimensional antiferromagnet with Dzyaloshinsky-Moriya interactions

We investigate solitary excitations in a model of a one-dimensional antiferromagnet including a single-ion anisotropy and a Dzyaloshinsky-Moriya antisymmetric exchange interaction term. We employ the Holstein-Primakoff transformation, the coherent state ansatz and the time variational principle. We obtain two partial differential equations of motion by using the method of multiple scales and applying perturbation theory. By so doing, we show that the motion of the coherent amplitude must satisfy the nonlinear Schrodinger equation. We give the single-soliton solution.

Traveling wave front for the Fisher equation on an infinite band region

Instead of discussing the existence of a one-dimensional traveling wave front solution which connects two constant steady states, the present work deals with the case connecting a constant and a nonhomogeneous steady state on an infinite band region. The corresponding model is the well-known Fisher equation with variational coefficient and Dirichlet boundary condition. (c) 2006 Elsevier Ltd. All rights reserved.

Algorithm Studies on How to Obtain a Conditional Nonlinear Optimal Perturbation (CNOP)

The conditional nonlinear optimal perturbation (CNOP), which is a nonlinear generalization of the linear singular vector (LSV), is applied in important problems of atmospheric and oceanic sciences, including ENSO predictability, targeted observations, and ensemble forecast. In this study, we investigate the computational cost of obtaining the CNOP by several methods. Differences and similarities, in terms of the computational error and cost in obtaining the CNOP, are compared among the sequential quadratic programming (SQP) algorithm, the limited memory Broyden-Fletcher-Goldfarb-Shanno (L-BFGS) algorithm, and the spectral projected gradients (SPG2) algorithm. A theoretical grassland ecosystem model and the classical Lorenz model are used as examples. Numerical results demonstrate that the computational error is acceptable with all three algorithms. The computational cost to obtain the CNOP is reduced by using the SQP algorithm. The experimental results also reveal that the L-BFGS algorithm is the most effective algorithm among the three optimization algorithms for obtaining the CNOP. The numerical results suggest a new approach and algorithm for obtaining the CNOP for a large-scale optimization problem.

Application of altimetry data assimilation on mesoscale eddies simulation

Mesoscale eddy plays an important role in the ocean circulation. In order to improve the simulation accuracy of the mesoscale eddies, a three-dimensional variation (3DVAR) data assimilation system called Ocean Variational Analysis System (OVALS) is coupled with a POM model to simulate the mesoscale eddies in the Northwest Pacific Ocean. In this system, the sea surface height anomaly (SSHA) data by satellite altimeters are assimilated and translated into pseudo temperature and salinity (T-S) profile data. Then, these profile data are taken as observation data to be assimilated again and produce the three-dimensional analysis T-S field. According to the characteristics of mesoscale eddy, the most appropriate assimilation parameters are set up and testified in this system. A ten years mesoscale eddies simulation and comparison experiment is made, which includes two schemes: assimilation and non-assimilation. The results of comparison between two schemes and the observation show that the simulation accuracy of the assimilation scheme is much better than that of non-assimilation, which verified that the altimetry data assimilation method can improve the simulation accuracy of the mesoscale dramatically and indicates that it is possible to use this system on the forecast of mesoscale eddies in the future.

基于变分的主动轮廓图像分割方法

图像分割是图像处理中很重要的一个问题，是计算机视觉的基础。因为它能够简化信息的存储和表示，从而能够对获取的图像内容进行智能解释，所以在很多应用问题中，图像分割是必不可少的过程，如医学图像处理，环境三维重建及自动目标识别等。图像分割的方法有很多种，如边缘检测，阈值，区域融合，分水岭及马尔可夫随机场等。虽然这些方法有其各自特点，但是它们在图象分割过程中不能充分将图像底层信息与高层信息结合，从而无法模拟人类视觉系统智能性。当图像底层信息不足时，这些仅基于数据驱动的分割模型无法达到令人满意的结果。尽管某种具体图像分割方法不可能满足所有图像分割要求，但利用尽可能多的高层与底层信息，将图像分割成有意义和人们所期望的区域始终是研究者所追求的目标。图像分割问题的数学建模和计算中有两个关键因素。第一是建立合适的分割模型将分割边界和分割区域的作用有效结合。第二是利用最有效的方法将分割边界和分割区域的几何特征统一到分割模型中。基于变分原理的主动轮廓图像分割将图像视为连续函数。这就使得研究者可以从连续函数空间角度来研究图像分割问题。这同时也为研究者提供严格的数学工具，如微分几何、泛函分析和微分方程等。为此它能很好的解决上述两个问题。第一，Mumford-Shah(M-S)模型为基于变分的主动轮廓分割模型提供了一完整的数学理论框架，并且Mumford-Shah模型从信息论的角度也能得到合理解释。第二，水平集方法能有效的表示分割边界和分割区域的几何特征。与其它方法相比，变分主动轮廓在理论和实际计算过程中都具有显著的优势。首先它能直接处理和表示各种重要的几何特征，如梯度、切向量、曲率等，并且能有效模拟很多动态过程，如线性和非线性扩散等。再则其可以利用很多已有的丰富数值方法进行分析和计算。本文基于变分原理与偏微分方程方法，利用主动轮廓模型具有结合底层图像信息与高层先验知识的特点，将特定先验知识与主动轮廓分割模型进行有效结合以弥补底层图像信息的不足，从而使主动轮分割廓模型具有更强的智能性。本文主要从两点对变分主动轮廓分割模型展开了研究：1、演化轮廓的形状约束；2、演化轮廓的梯度下降流约束及其滤波实现。其主要工作包括以下四个方面的内容：第一，基于M-S模型和样条曲线的开边界检测。本章通过对演化轮廓引入合理边界条件，利用样条曲线表示待检测的开曲线，将一般开曲线的检测问题变为合理的图像分割问题，从而达到一般开曲线检测目的。此方法称为开扩散蛇模型。一般开曲线的检测具有很多应用领域，如：河流、道路、天际线、焊缝等检测。第二，方差主动轮廓模型。在目标跟踪应用中，跟踪目标的主要运动形式表现为平移。本章将此作为一种先验知识与主动轮廓模型结合，提出了一种方差主动轮廓模型(HV)。此模型的特点是轮廓在演化过程中具有平移优先和快速的良好特性。它比已有的主动轮廓模型更适于自动目标跟踪领域。第三，基于M-S模型和隐式曲面变分方法的一般梯度下降流滤波器。本章为一般梯度下降流求取提供了统一框架及解决方法。首先本章将H0梯度下降流和一般梯度下降流统一到Mumford-Shah模型框架中，从而将一般梯度下降流的求取转换为一个极小化泛函问题，并利用隐式曲面变分方法对此极小化泛函进行求解。另外本章从滤波器设计角度出发，通过对H0梯度下降流滤波可以得到一般梯度下降流。滤波器的实现体现了内嵌于一般梯度下降流的先验属性。根据此思想，本章将对应于HV和H1主动轮廓的內积空间顺序组合，对H0梯度下降流进行顺序滤波，提出了一种既具有全局平移优先性又具有局部光滑速度场的主动轮廓，称为HV1主动轮廓。它将H0，H1和HV主动轮廓统一起来。第四，形状保持主动轮廓模型及其应用。针对某些特定目标的检测，本章提出了形状保持主动轮廓模型。此模型能够达到分割即目标的目的，同时能够给出目标的定量描述。基于此模型，本章实现了具有椭圆、直线和平行四边形轮廓特征目标的检测。椭圆形状约束用于眼底图像分割。直线和平行四边行分别用于自动目标识别中的天际线检测和机场跑道跟踪。

基于多Agent制造过程的建模方法

分析了制造系统与制造过程之间的关系;论证了从过程的角度对制造进行建模更恰当;结合Agent和π演算的特点,给出Agent制造系统描述模型及基于π演算的单个Agent的BDI模型,并指出Agent和π演算结合的制造过程模型有利于进行优化目标在不同制造过程层次的分解,不论从方法的角度还是实现的角度,都适合复杂系统建模。Agent和π演算相结合可以有效分析并解决离散事件的建模与仿真中的问题。

油气勘探中的弱磁异常提取与区域场分析

Petroleum and natural gas is an important strategic resources. The short of the reserves will block the development of economy and threaten the safety of nation, along with the main oil fields of our country coming to the height of power and splendor of the exploitation and exploration. Therefore, it makes a great sense to inaugurate new explorative field and increase the reserves of petroleum and natural gas. Magnetic exploration is a main method of geophysics exploration. the developing observation apparatus and the perfect processing method provide wide space for magnetic exploration in these years. The method of magnetic bright spot is an application of magnetic exploration. The vertical migration of the hydrocarbon changes physical and chemical environment above the hydrocarbon reservoir, the new environment make tervalent iron translate into bivalent iron, that produce small scale magnetic anomaly, that is magnetic bright spot. The method of magnetic bright spot explores oil and gas field by the relation between the hydrocarbon and magnetic anomaly. This paper systemically research to pick-up and identify magnetic bright spot combining an oil field item, then point out advantaged area. In order to test the result, the author use the seismic information to superpose the magnetic bright spot, that prove the magnetic bright spot is reliable. then, the author complete a software to pick and identify the magnetic bright spot. The magnetic basement is very important to research forming and evolvement of the basin, especially, it is a crucial parameter of exploring residual basin in the research on pre-Cenozoic residual. This paper put forward a new method to inverse the interface of the magnetic layer on the basis of previous work, that is the method of separation of magnetic field step by step. The theory of this method is to translate the result of magnetic layer fluctuation to the result of magnetization density change, and the magnetic layer is flat, the paper choose thickness of magnetic layer as unit thickness, and define magnetic layer as a unit-thickness layer in order to convenient calculation, at the same time, define the variational magnetization density as equivalent magnetic density. Then we translate the relation between magnetic field and layer fluctuation to the relation between magnetic field and equivalent magnetic density, then, we can obtain the layer fluctuation through calculating equivalent magnetic density. Contrast to conventional parker method, model experimentation and example checkout prove this method is effective. The merit of this method is to avoid flat result in a strongly fluctuant area because of using a uniform average depth, the result of this method is closer to the fact, and this method is to inverse equivalent magnetic density, then translate equivalent magnetic density to layer fluctuation, this lays a foundation to inverse variational magnetic density in the landscape orientation and portrait.

基于量子蒙特卡罗的地震波反演方法

The primary approaches for people to understand the inner properties of the earth and the distribution of the mineral resources are mainly coming from surface geology survey and geophysical/geochemical data inversion and interpretation. The purpose of seismic inversion is to extract information of the subsurface stratum geometrical structures and the distribution of material properties from seismic wave which is used for resource prospecting, exploitation and the study for inner structure of the earth and its dynamic process. Although the study of seismic parameter inversion has achieved a lot since 1950s, some problems are still persisting when applying in real data due to their nonlinearity and ill-posedness. Most inversion methods we use to invert geophysical parameters are based on iterative inversion which depends largely on the initial model and constraint conditions. It would be difficult to obtain a believable result when taking into consideration different factors such as environmental and equipment noise that exist in seismic wave excitation, propagation and acquisition. The seismic inversion based on real data is a typical nonlinear problem, which means most of their objective functions are multi-minimum. It makes them formidable to be solved using commonly used methods such as general-linearization and quasi-linearization inversion because of local convergence. Global nonlinear search methods which do not rely heavily on the initial model seem more promising, but the amount of computation required for real data process is unacceptable. In order to solve those problems mentioned above, this paper addresses a kind of global nonlinear inversion method which brings Quantum Monte Carlo (QMC) method into geophysical inverse problems. QMC has been used as an effective numerical method to study quantum many-body system which is often governed by Schrödinger equation. This method can be categorized into zero temperature method and finite temperature method. This paper is subdivided into four parts. In the first one, we briefly review the theory of QMC method and find out the connections with geophysical nonlinear inversion, and then give the flow chart of the algorithm. In the second part, we apply four QMC inverse methods in 1D wave equation impedance inversion and generally compare their results with convergence rate and accuracy. The feasibility, stability, and anti-noise capacity of the algorithms are also discussed within this chapter. Numerical results demonstrate that it is possible to solve geophysical nonlinear inversion and other nonlinear optimization problems by means of QMC method. They are also showing that Green’s function Monte Carlo (GFMC) and diffusion Monte Carlo (DMC) are more applicable than Path Integral Monte Carlo (PIMC) and Variational Monte Carlo (VMC) in real data. The third part provides the parallel version of serial QMC algorithms which are applied in a 2D acoustic velocity inversion and real seismic data processing and further discusses these algorithms’ globality and anti-noise capacity. The inverted results show the robustness of these algorithms which make them feasible to be used in 2D inversion and real data processing. The parallel inversion algorithms in this chapter are also applicable in other optimization. Finally, some useful conclusions are obtained in the last section. The analysis and comparison of the results indicate that it is successful to bring QMC into geophysical inversion. QMC is a kind of nonlinear inversion method which guarantees stability, efficiency and anti-noise. The most appealing property is that it does not rely heavily on the initial model and can be suited to nonlinear and multi-minimum geophysical inverse problems. This method can also be used in other filed regarding nonlinear optimization.

地震波模拟：非规则网格离散的数值解法、吸收边界及其实际应用

This dissertation presents a series of irregular-grid based numerical technique for modeling seismic wave propagation in heterogeneous media. The study involves the generation of the irregular numerical mesh corresponding to the irregular grid scheme, the discretized version of motion equations under the unstructured mesh, and irregular-grid absorbing boundary conditions. The resulting numerical technique has been used in generating the synthetic data sets on the realistic complex geologic models that can examine the migration schemes. The motion equation discretization and modeling are based on Grid Method. The key idea is to use the integral equilibrium principle to replace the operator at each grid in Finite Difference scheme and variational formulation in Finite Element Method. The irregular grids of complex geologic model is generated by the Paving Method, which allow varying grid spacing according to meshing constraints. The grids have great quality at domain boundaries and contain equal quantities of nodes at interfaces, which avoids the interpolation of parameters and variables. The irregular grid absorbing boundary conditions is developed by extending the Perfectly Matched Layer method to the rotated local coordinates. The splitted PML equations of the first-order system is derived by using integral equilibrium principle. The proposed scheme can build PML boundary of arbitrary geometry in the computational domain, avoiding the special treatment at corners in a standard PML method and saving considerable memory and computation cost. The numerical implementation demonstrates the desired qualities of irregular grid based modeling technique. In particular, (1) smaller memory requirements and computational time are needed by changing the grid spacing according to local velocity; (2) Arbitrary surfaces and interface topographies are described accurately, thus removing the artificial reflection resulting from the stair approximation of the curved or dipping interfaces; (3) computational domain is significantly reduced by flexibly building the curved artificial boundaries using the irregular-grid absorbing boundary conditions. The proposed irregular grid approach is apply to reverse time migration as the extrapolation algorithm. It can discretize the smoothed velocity model by irregular grid of variable scale, which contributes to reduce the computation cost. The topography. It can also handle data set of arbitrary topography and no field correction is needed.