车载容迟网络中一种基于真实轨迹的车辆移动模型

车辆移动模型的选择会严重影响车载容迟网络中关键技术的仿真研究结果.现有车辆移动模型大多属于合成移动模型,不符合真实环境.基于车辆真实移动轨迹建立车辆移动模型,考虑了道路拓扑选择、车辆速度控制、车辆初始位置选择以及路口转向选择等因素.在此基础上,进一步分析了模型对实际车辆行为模拟的接近程度、车辆的分布特性,以及基于该模型的网络连接特性.结果表明,提出的模型能够真实地逼近现实环境中车辆的实际行为.

Effect of poly(propylene carbonate) on the crystallization and melting behavior of poly(beta-hydroxybutyrate-co-beta-hydroxyvalerate)

The crystallization and melting behavior of poly(beta-hydroxybutyrate-co-beta-hydroxyvalerate) (PHBV) and a 30/70 (w/w) PHBV/poly(propylene carbonate) (PPC) blend was investigated with differential scanning calorimetry (DSC) and Fourier transform infrared (FTIR). The transesterification reaction between PHBV and PPC was detected in the melt-blending process. The interaction between the two macromolecules was confirmed by means of FTIR analysis. During the crystallization process from the melt, the crystallization temperature of the PHBV/PPC blend decreased about 8 degreesC, the melting temperature was depressed by 4 degreesC, and the degree of crystallinity of PHBV in the blend decreased about 9.4%; this was calculated through a comparison of the DSC heating traces for the blend and pure PHBV. These results indicated that imperfect crystals of formed, crystallization was inhibited, and the crystallization ability of PHBV was weakened in the blend. The equilibrium melting temperatures of PHBV and the 30/70 PHBV/PPC blend isothermally crystallized were 187.1 and 179 degreesC, respectively.

Structure and properties of the ordered double perovskites Sr2MWO6 (M = Co, Ni) by sol-gel route

The single-phase double perovskites Sr2MWO6 (M=Co, Ni) were prepared by sol-gel method. Crystal Structure, magnetic properties and the morphology of Sr2CoWO6 and Sr2NiWO6 were investigated. X-ray powder diffraction (XRD) analysis shows single phase structure for Sr2MWO6 (M=Co, Ni) without any traces of impurities and the crystal structure of all the samples belongs to the tetragonal I4/m space group. SEM image for Sr2MWO6 (M=Co, Ni) indicate that the grains are homogeneous and connect each other very well. The Neel temperature for Sr2CoWO6 and Sr2NiWO6 are 23 K and 59 K, respectively. Magnetic measurements showed that the magnetic moment in these double perovskites originates mainly from the interactions between Ni ions and Co ions.

A study on the catalytic behaviors of concentrated heteropolyacid solution - I. A novel catalyst for isobutane alkylation with butene

The catalytic behaviors of a novel liquid acid catalyst (composed of heteropolyacid and acetic acid) for alkylation of isobutane with butene was investigated. As a solvent acetic acid had a synergistic effect. It enhanced the acid strength of HPA and its stability. The conditions for the formation of the catalytically active phase were studied systematically. The content of crystal water of HPA and the quantity of solvent affect the formation of active phase and the catalytic activity. Catalytically active phase consists of HPA, acetic acid and hydrocarbon produced from the reaction, as well as traces of water from the crystal water of HPA. This catalyst system is comparable to the sulfuric acid in catalytic activity.

Radiation effects on the immiscible polymer blend of nylon1010 and high-impact polystyrene (HIPS) I: Gel/dose curves, mathematical expectation theorem and thermal behaviour

This paper studies the radiation properties of the immiscible blend of nylon1010 and HIPS. The gel fraction increased with increasing radiation dose. The network was found mostly in nylon1010, the networks were also found in both nylon1010 and HIPS when the dose reaches 0.85 MGy or more. We used the Charleby-Pinner equation and the modified Zhang-Sun-Qian equation to simulate the relationship with the dose and the sol fraction. The latter equation fits well with these polymer blends and the relationship used by it showed better linearity than the one by the Charleby-Pinner equation. We also studied the conditions of formation of the network by the mathematical expectation theorem for the binary system. Thermal properties of polymer blend were observed by DSC curves. The crystallization temperature decreases with increasing dose because the cross-linking reaction inhibited the crystallization procession and destroyed the crystals. The melting temperature also reduced with increasing radiation dose. The dual melting peak gradually shifted to single peak and the high melting peak disappeared at high radiation dose. However, the radiation-induced crystallization was observed by the heat of fusion increasing at low radiation dose. On the other hand, the crystal will be damaged by radiation. A similar conclusion may be drawn by the DSC traces when the polymer blends were crystallized. When the radiation dose increases, the heat of fusion reduces dramatically and so does the heat of crystallization. (C) 1999 Elsevier Science Ltd. All rights reserved.

Synthesis and characterization of thermotropic liquid crystalline poly(ester imide)s

Two closely series of poly(ester imide)s had been synthesized by solution polycondensation of p-phenylenebis(trimellitate) dianhydride with aliphatic diamines. The differential scanning calorimetry (DSC) traces of the most poly(ester imide)s exhibited two endotherms representing the solid state to anisotropic phase transition (T-m1) and the anisotropic to isotropic melt transition (T-m2), respectively. Observation under polarizing microscope and wide-angle X-ray diffraction (WAXD) measurements suggested that the anisotropic phase formed above the melting paints (T-m1) had a smectic character. The thermogravimetric analyses (TGA) revealed that the thermal stabilities of the poly(ester imide)s were up to 350 degrees C. (C) 1999 John Wiley & Sons, Inc.

Crystal structure of polypropylene filled with rare earth oxides

The effect of a fine powder of Y2O3, Nd2O3, and Ho2O3 on the crystal structure of isotactic polypropylene (iPP) was studied with WAXD and DSC techniques. The results showed that the addition of the three rare earth oxides (REOs) can increase the crystallite size of the alpha-form crystal and the degree of crystallinity of iPP at an annealing temperature of 120 degrees C and that both Y2O3 and Nd2O3 are the beta-nucleator of iPP. REOs enhance the overall growth rate of the spherulites of iPP. All the iPP samples filled with REOs which were crystallized isothermally at 132 degrees C from the melt exhibited their melting peaks of the beta-form on the DSC heating traces, indicating that the REOs are the nucleating agents for both the alpha- and beta-forms of iPP under isothermal conditions. (C) 1996 John Wiley & Sons, Inc.

Synthesis and properties of thermotropic liquid crystalline copolyesters containing bis(4-hydroxyphenyl)methanone(V)

A new series of fully aromatic thermotropic liquid crystalline copolyesters were synthesized by thermal polycondensation of p-hydroxybenzoic acid (PHB), terephthalic acid (TPA), resorcinol (RES) and bis (4-hydroxyphenyl) methanone (BHP). The DSC traces of copolyesters exhibited two or three endotherms, the low temperature endotherm(Tm) of which indicated the transition from solid state to anisotropic phase. Observation under polarizing microscope and WAXD measurements suggested that the LC-phase formed immediately above Tm had a nematic character.

碎屑岩储层低电阻率油层识别方法与应用——以济阳坳陷为例

So far, there is no methods of logging interpretation effective enough to identify a low resistivity payzone since its resistivity value almost equals to that of an aquifer although many low-resistivity payzones have been found in lots of petroliferous basins worldwide. After a thorough study on those technical difficulties of the logging interpretation for the low-resistivity payzones, some corresponding resolutions have been put forward in this paper. In order to reveal its microscopic mechanism, researches on the discovered low-resistivity payzones have been carried on with analyses of core and lab test data, thus main influencing factors of the low-resistivity reservoirs have been pointed out including conductivity minerals, clay minerals, fluids, porosity and pore structure. In order to make clear the degree of influence of those reservoir factors on resistivity logging(log), lab studies and numeral simulations have been done with the typical core and formation water samples, therefore, their influence degrees have ascertained quantitatively or semi-quantitatively. The distribution law and possible distribution areas of the low-resistivity payzones in Jiyang Depression have been figured out firstly after the macroscopic geology origins (sedimentation, dynamic accumulation process, diagenesis etc.) in the area have been studied. In order to resolve the problem of difficult logging-interpretation, methods of interpretation and identification have been brought forward creatively according to the low-resistivity payzone type ascribed to macroscopic geology laws and to the combined features of logging traces, after a systemic summary of different responses of logging caused by different microscopic mechanism. Those methods have been applied in Dongying and Huimin Sag of Shengli Exploration Area, precision of identification of the low-resistivity payzones improved effectively and good economic attraction prove their great prospect.

叠前地震数据规则化、重建及噪音压制

I address of reconstruction of spatial irregular sampling seismic data to regular grids. Spatial irregular sampling data impairs results of prestack migration, multiple attenuations, spectra estimation. Prestack 5-D volumes are often divided into sub-sections for further processing. Shot gathers are easy to obtain from irregular sampling volumes. My strategy for reconstruction is as follows: I resort irregular sampling gathers into a form of easy to bin and perform bin regularization, then utilize F-K inversion to reconstruct seismic data. In consideration of poor ability of F-K regularization to fill in large gaps, I sort regular sampling gathers to CMP and proposed high-resolution parabolic Radon transform to interpolate data and extrapolate offsets. To strong interfering noise--multiples, I use hybrid-domain high-resolution parabolic Radon transform to attenuate it. F-K regularization demand ultimately for lower computing costs. I proposed several methods to further improve efficiency of F-K inversion: first I introduce 1D and 2D NFFT algorithm for a rapid calculation of DFT operators; then develop fast 1D and 2D CG method to solve least-square equations, and utilize preconditioner to accelerate convergence of CG iterations; what’s more, I use Delaunay triangulation for weight calculation and use bandlimit frequency and varying bandwidth technique for competitive computation. Numerical 2D and 3D examples are offered to verify reasonable results and more efficiency. F-K regularization has poor ability to fill in large gaps, so I rearrange data as CMP gathers and develop hybrid-domain high-resolution parabolic Radon transforms which be used ether to interpolate null traces and extrapolate near and far offsets or suppress a strong interfere noise: multiples. I use it to attenuate multiples to verify performances of our algorithm and proposed routines for industrial application. Numerical examples and field data examples show a nice performance of our method.

基于显式分形插值的高保真地震数据重建方法研究

Based on the fractal theories, contractive mapping principles as well as the fixed point theory, by means of affine transform, this dissertation develops a novel Explicit Fractal Interpolation Function（EFIF）which can be used to reconstruct the seismic data with high fidelity and precision. Spatial trace interpolation is one of the important issues in seismic data processing. Under the ideal circumstances, seismic data should be sampled with a uniform spatial coverage. However, practical constraints such as the complex surface conditions indicate that the sampling density may be sparse or for other reasons some traces may be lost. The wide spacing between receivers can result in sparse sampling along traverse lines, thus result in a spatial aliasing of short-wavelength features. Hence, the method of interpolation is of very importance. It not only needs to make the amplitude information obvious but the phase information, especially that of the point that the phase changes acutely. Many people put forward several interpolation methods, yet this dissertation focuses attention on a special class of fractal interpolation function, referred to as explicit fractal interpolation function to improve the accuracy of the interpolation reconstruction and to make the local information obvious. The traditional fractal interpolation method mainly based on the randomly Fractional Brown Motion (FBM) model, furthermore, the vertical scaling factor which plays a critical role in the implementation of fractal interpolation is assigned the same value during the whole interpolating process, so it can not make the local information obvious. In addition, the maximal defect of the traditional fractal interpolation method is that it cannot obtain the function values on each interpolating nodes, thereby it cannot analyze the node error quantitatively and cannot evaluate the feasibility of this method. Detailed discussions about the applications of fractal interpolation in seismology have not been given by the pioneers, let alone the interpolating processing of the single trace seismogram. On the basis of the previous work and fractal theory this dissertation discusses the fractal interpolation thoroughly and the stability of this special kind of interpolating function is discussed, at the same time the explicit presentation of the vertical scaling factor which controls the precision of the interpolation has been proposed. This novel method develops the traditional fractal interpolation method and converts the fractal interpolation with random algorithms into the interpolation with determined algorithms. The data structure of binary tree method has been applied during the process of interpolation, and it avoids the process of iteration that is inevitable in traditional fractal interpolation and improves the computation efficiency. To illustrate the validity of the novel method, this dissertation develops several theoretical models and synthesizes the common shot gathers and seismograms and reconstructs the traces that were erased from the initial section using the explicit fractal interpolation method. In order to compare the differences between the theoretical traces that were erased in the initial section and the resulting traces after reconstruction on waveform and amplitudes quantitatively, each missing traces are reconstructed and the residuals are analyzed. The numerical experiments demonstrate that the novel fractal interpolation method is not only applicable to reconstruct the seismograms with small offset but to the seismograms with large offset. The seismograms reconstructed by explicit fractal interpolation method resemble the original ones well. The waveform of the missing traces could be estimated very well and also the amplitudes of the interpolated traces are a good approximation of the original ones. The high precision and computational efficiency of the explicit fractal interpolation make it a useful tool to reconstruct the seismic data; it can not only make the local information obvious but preserve the overall characteristics of the object investigated. To illustrate the influence of the explicit fractal interpolation method to the accuracy of the imaging of the structure in the earth’s interior, this dissertation applies the method mentioned above to the reverse-time migration. The imaging sections obtained by using the fractal interpolated reflected data resemble the original ones very well. The numerical experiments demonstrate that even with the sparse sampling we can still obtain the high accurate imaging of the earth’s interior’s structure by means of the explicit fractal interpolation method. So we can obtain the imaging results of the earth’s interior with fine quality by using relatively small number of seismic stations. With the fractal interpolation method we will improve the efficiency and the accuracy of the reverse-time migration under economic conditions. To verify the application effect to real data of the method presented in this paper, we tested the method by using the real data provided by the Broadband Seismic Array Laboratory, IGGCAS. The results demonstrate that the accuracy of explicit fractal interpolation is still very high even with the real data with large epicenter and large offset. The amplitudes and the phase of the reconstructed station data resemble the original ones that were erased in the initial section very well. Altogether, the novel fractal interpolation function provides a new and useful tool to reconstruct the seismic data with high precision and efficiency, and presents an alternative to image the deep structure of the earth accurately.

积分偏移的计算几何与并行实现

The Second Round of Oil & Gas Exploration needs more precision imaging method, velocity vs. depth model and geometry description on Complicated Geological Mass. Prestack time migration on inhomogeneous media was the technical basic of velocity analysis, prestack time migration on Rugged surface, angle gather and multi-domain noise suppression. In order to realize this technique, several critical technical problems need to be solved, such as parallel computation, velocity algorithm on ununiform grid and visualization. The key problem is organic combination theories of migration and computational geometry. Based on technical problems of 3-D prestack time migration existing in inhomogeneous media and requirements from nonuniform grid, parallel process and visualization, the thesis was studied systematically on three aspects: Infrastructure of velocity varies laterally Green function traveltime computation on ununiform grid, parallel computational of kirchhoff integral migration and 3D visualization, by combining integral migration theory and Computational Geometry. The results will provide powerful technical support to the implement of prestack time migration and convenient compute infrastructure of wave number domain simulation in inhomogeneous media. The main results were obtained as follows: 1. Symbol of one way wave Lie algebra integral, phase and green function traveltime expressions were analyzed, and simple 2-D expression of Lie algebra integral symbol phase and green function traveltime in time domain were given in inhomogeneous media by using pseudo-differential operators’ exponential map and Lie group algorithm preserving geometry structure. Infrastructure calculation of five parts, including derivative, commutating operator, Lie algebra root tree, exponential map root tree and traveltime coefficients , was brought forward when calculating asymmetry traveltime equation containing lateral differential in 3-D by this method. 2. By studying the infrastructure calculation of asymmetry traveltime in 3-D based on lateral velocity differential and combining computational geometry, a method to build velocity library and interpolate on velocity library using triangulate was obtained, which fit traveltime calculate requirements of parallel time migration and velocity estimate. 3. Combining velocity library triangulate and computational geometry, a structure which was convenient to calculate differential in horizontal, commutating operator and integral in vertical was built. Furthermore, recursive algorithm, for calculating architecture on lie algebra integral and exponential map root tree (Magnus in Math), was build and asymmetry traveltime based on lateral differential algorithm was also realized. 4. Based on graph theory and computational geometry, a minimum cycle method to decompose area into polygon blocks, which can be used as topological representation of migration result was proposed, which provided a practical method to block representation and research to migration interpretation results. 5. Based on MPI library, a process of bringing parallel migration algorithm at arbitrary sequence traces into practical was realized by using asymmetry traveltime based on lateral differential calculation and Kirchhoff integral method. 6. Visualization of geological data and seismic data were studied by the tools of OpenGL and Open Inventor, based on computational geometry theory, and a 3D visualize system on seismic imaging data was designed.

基于波动方程的层间多次波预测方法研究

In exploration geophysics，velocity analysis and migration methods except reverse time migration are based on ray theory or one-way wave-equation. So multiples are regarded as noise and required to be attenuated. It is very important to attenuate multiples for structure imaging, amplitude preserving migration. So it is an interesting research in theory and application about how to predict and attenuate internal multiples effectively. There are two methods based on wave-equation to predict internal multiples for pre-stack data. One is common focus point method. Another is inverse scattering series method. After comparison of the two methods, we found that there are four problems in common focus point method: 1. dependence of velocity model; 2. only internal multiples related to a layer can be predicted every time; 3. computing procedure is complex; 4. it is difficult to apply it in complex media. In order to overcome these problems, we adopt inverse scattering series method. However, inverse scattering series method also has some problems: 1. computing cost is high; 2. it is difficult to predict internal multiples in the far offset; 3. it is not able to predict internal multiples in complex media. Among those problems, high computing cost is the biggest barrier in field seismic processing. So I present 1D and 1.5D improved algorithms for reducing computing time. In addition, I proposed a new algorithm to solve the problem which exists in subtraction, especially for surface related to multiples. The creative results of my research are following: 1. derived an improved inverse scattering series prediction algorithm for 1D. The algorithm has very high computing efficiency. It is faster than old algorithm about twelve times in theory and faster about eighty times for lower spatial complexity in practice; 2. derived an improved inverse scattering series prediction algorithm for 1.5D. The new algorithm changes the computing domain from pseudo-depth wavenumber domain to TX domain for predicting multiples. The improved algorithm demonstrated that the approach has some merits such as higher computing efficiency, feasibility to many kinds of geometries, lower predictive noise and independence to wavelet； 3. proposed a new subtraction algorithm. The new subtraction algorithm is not used to overcome nonorthogonality, but utilize the nonorthogonality's distribution in TX domain to estimate the true wavelet with filtering method. The method has excellent effectiveness in model testing. Improved 1D and 1.5D inverse scattering series algorithms can predict internal multiples. After filtering and subtracting among seismic traces in a window time, internal multiples can be attenuated in some degree. The proposed 1D and 1.5D algorithms have demonstrated that they are effective to the numerical and field data. In addition, the new subtraction algorithm is effective to the complex theoretic models.

消除采集脚印的处理方法研究

Seismic technique is in the leading position for discovering oil and gas trap and searching for reserves throughout the course of oil and gas exploration. It needs high quality of seismic processed data, not only required exact spatial position, but also the true information of amplitude and AVO attribute and velocity. Acquisition footprint has an impact on highly precision and best quality of imaging and analysis of AVO attribute and velocity. Acquisition footprint is a new conception of describing seismic noise in 3-D exploration. It is not easy to understand the acquisition footprint. This paper begins with forward modeling seismic data from the simple sound wave model, then processes it and discusses the cause for producing the acquisition footprint. It agreed that the recording geometry is the main cause which leads to the distribution asymmetry of coverage and offset and azimuth in different grid cells. It summarizes the characters and description methods and analysis acquisition footprint’s influence on data geology interpretation and the analysis of seismic attribute and velocity. The data reconstruct based on Fourier transform is the main method at present for non uniform data interpolation and extrapolate, but this method always is an inverse problem with bad condition. Tikhonov regularization strategy which includes a priori information on class of solution in search can reduce the computation difficulty duo to discrete kernel condition disadvantage and scarcity of the number of observations. The method is quiet statistical, which does not require the selection of regularization parameter; and hence it has appropriate inversion coefficient. The result of programming and tentat-ive calculation verifies the acquisition footprint can be removed through prestack data reconstruct. This paper applies migration to the processing method of removing the acquisition footprint. The fundamental principle and algorithms are surveyed, seismic traces are weighted according to the area which occupied by seismic trace in different source-receiver distances. Adopting grid method in stead of accounting the area of Voroni map can reduce difficulty of calculation the weight. The result of processing the model data and actual seismic demonstrate, incorporating a weighting scheme based on the relative area that is associated with each input trace with respect to its neighbors acts to minimize the artifacts caused by irregular acquisition geometry.

数据驱动型多次波衰减方法的研究

In the practical seismic profile multiple reflections tend to impede the task of even the experienced interpreter in deducing information from the reflection data. Surface multiples are usually much stronger, more broadband, and more of a problem than internal multiples because the reflection coefficient at the water surface is much larger than the reflection coefficients found in the subsurface. For this reason most attempts to remove multiples from marine data focus on surface multiples, as will I. A surface-related multiple attenuation method can be formulated as an iterative procedure. In this essay a fully data-driven approach which is called MPI —multiple prediction through inversion (Wang, 2003) is applied to a real marine seismic data example. This is a pretty promising scheme for predicting a relative accurate multiple model by updating the multiple model iteratively, as we usually do in a linearized inverse problem. The prominent characteristic of MPI method lie in that it eliminate the need for an explicit surface operator which means it can model the multiple wavefield without any knowledge of surface and subsurface structures even a source signature. Another key feature of this scheme is that it can predict multiples not only in time but also in phase and in amplitude domain. According to the real data experiments it is shown that this scheme for multiple prediction can be made very efficient if a good initial estimate of the multiple-free data set can be provided in the first iteration. In the other core step which is multiple subtraction we use an expanded multi-channel matching filter to fulfil this aim. Compared to a normal multichannel matching filter where an original seismic trace is matched by a group of multiple-model traces, in EMCM filter a seismic trace is matched by not only a group of the ordinary multiple-model traces but also their adjoints generated mathematically. The adjoints of a multiple-model trace include its first derivative, its Hilbert transform and the derivative of the Hilbert transform. The third chapter of the thesis is the application for the real data using the previous methods we put forward from which we can obviously find the effectivity and prospect of the value in use. For this specific case I have done three group experiments to test the effectiveness of MPI method, compare different subtraction results with fixed filter length but different window length, invest the influence of the initial subtraction result for MPI method. In terms of the real data application, we do fine that the initial demultiple estimate take on a great deal of influence for the MPI method. Then two approaches are introduced to refine the intial demultiple estimate which are first arrival and masking filter respectively. In the last part some conclusions are drawn in terms of the previous results I have got.