数字图像处理中多窗口下的自适应中值滤波

根据心理物理学、神经生理学、认知神经学等学科在视觉认知领域的部分研究成果，结合机器视觉、图像处理领域在图像增强方面已经提出的一些方法，提出了结合先验知识的多窗口结构下的分块中值滤波方法，在每一个窗口内单独进行处理与分析，突出了视觉处理目的，减少了运算量， 节省数据存储空间，达到了令人满意的滤波效果，能够在原始图像比较复杂的情况下，较好地对其进行预处理，可以改善、提高后期图像处理过程，如图像分割、图像分析的正确性和有效性。

用光学阵列传感器的机器人物体分类系统

本文介绍用光学阵列传感器的机器人物体分类系统。传感器直接安装在机器人的两个手指上。被抓物体的阴影通过光导纤维传到安放在“安全区”的光敏元件上。计算机识别物体的轮廓后命令机器人抓握物体,并把它运送到指定的地点从而达到物体分类的目的。

PASSIM卷接机组新型电控系统的设计与实现

PASSIM卷接机组原控制系统采用电路板进行逻辑控制,抗干扰性能差,故障率高,造成生产效率降低、原材料消耗增大、维修工作繁重等。为此,设计了一套新型卷接机组电气控制系统。该系统采用PLC进行过程控制,以工控机为上位机完成人机通信,采用交流伺服驱动,并通过高速信号处理专用系统完成重量检测控制及烟支质量检测功能;利用PROFIBUS、CAN及MPI多种总线方式完成各单元间的通讯,实现信号和数据间的传递和共享。改进后的PASSIM机组运行稳定可靠,采样速度快,实时性强,且维修方便。机组的有效作业率由85%左右提高到90%以上,降低了卷烟纸和烟丝等原材料的消耗。

基于FPGA高速并行FFT处理器设计与实现

在研究快速傅里叶变换(FFT)算法的基础上,根据FPGA性能高、灵活性强、速度快的特点,提出了高效的基4-FFT处理器的实现方法。数据存储采用分块存储的方法,大大提高了存取速度。数据寻址采用新型的地址产生方法,可并行产生所需数据地址。同时,在蝶形单元的设计中很好的将并行运算技术和流水线技术相结合了起来,又进一步提高了运算速度。测试结果表明,时钟在50MHz时完成1024点FFT的时间为25.6μs,满足了应用实时性的要求。

基于CAN总线的分布式数据采集与控制系统

随着电子技术和计算机技术的不断发展，工业生产过程的控制系统正在向着智能化、数字化和网络化的方向发展。传统的集散控制方式和计算机分层控制方式已经开始让位于智能终端与网络结合的总线网络控制方式。当今，在工厂中过程控制环境下的分布式自动化系统变得越来越复杂，尤其系统内部的各设备之间需要快速交换大量的信息，以便实现对被控系统更为精确的控制和提供一些辅助的评价函数。这就意味着要不断增加带宽和提高通信速率以满足网络通信的需要。在现有的多种可利用网络设备中，CAN总线以其清晰的定义、极高的可靠性及其独特的设计，被认为是最能有效地解决这一问题的途径之一。而且市场上基于通信技术的产品中，就实时性考虑，由于CAN总线采用的非表意性的通信方式，因此其结构更为简单，实时性更好。基于此背景，我们以CAN总线作为通信媒介，将分布于各控制现场的传感器、执行器和控制器有序地连接起来，构成了一个基于CAN总线的分布式局域网络控制系统。本文首先介绍了基于CAN总线的分布式数据采集与控制系统的总体结构。然后从硬件方面描述了基于CAN总线的通信协议转换单元、数据采集单元和输出控制单元的功能、硬件配置及各单元功能的具体实现过程，给出了各单元的性能指标。软件方面，以C语言作为平台，开发了基于CAN总线的上位计算机管理与监控软件，实现了对整个网络设备的系统管理和系统控制功能。对于该总线系统，作者运用了PID控制和模糊控制算法实现了对水箱液位的控制，达到了理想的效果。基于CAN总线的控制系统很好地解决了集散控制系统难以解决的难题，模糊控制的应用能很好地把总线控制系统应用到具有非线性、大时滞和难于获得精确模型的控制系统中。

柴达木盆地地震勘探关键技术研究

The rugged surface topography determined the seismic data acquisition construction conditions and the seismic wave explosive and receiver quality in Qaidam Basin. This dissertation systematically researched the seismic acquisition, imaging process and the attribute analysis techniques of complicated oil and gas reservoir. The main research achievements and cognitions are as follows: 1. Through the stimulation effects research and analysis from the aspect of lithologic water-containing differences, it’s specific that stable hydrous sand layer can effectively enhance the stimulation effects combined with the corresponding field tests. The seismic data S/N ratio has been improved due to the combination explosive stimulation. Through the fold number and maximum offset analyses of target horizon, the complicated geometry has been optimized and the S/N ratio of seismic data has been improved, which made an important basis for improvement of 3D seismic data. 2. It has been proved that the first arrival refraction static correction method under the model constraint of fine surface survey is suitable to the Qaidam Basin of western areas by the real seismic data processing. Although the refraction horizon of near surface has some changes in a certain extent, it’s steady basically. The refraction horizon can be continuously traced in sections, so it’s qualified for the refraction static correction method on the whole. 3. The research is based on the curved-ray pre-stack time migration techniques of rough topography, and improved the imaging precision of complex areas. This techniques adopted the constant and variable velocity scanning mode and enhanced the velocity analysis precision. The 3D pre-stack time migration techniques reasonably solved the imaging and velocity multiple solutions problems of steep-dip faults and the intersections of horizontal layers. What’s more, fine velocity analysis and mute are very important to enhance the imaging precision of the seismic data in complicated Wunan areas. 4. The 3D seismic data edge-preserving processing methods have been realized due to the image process techniques. Because this method uses the large range filter, it can attenuate the noise maximally. The faults, break points, lithologic pinchout points and lithologic body of small scale such as river will not be influenced by blur because of the edge-preserving characterization of the method which is really an effective assistant technique of low S/N ratio seismic data attribute analysis. 5. The use of spectral decomposition technique can effectively identify the reservoirs. The special geology body which will not be identified (or without obvious characters) in the seismic profile may be found through the details changes of different frequencies in the amplitude profiles.

层析静校正技术研究和应用

Static correction is one of the indispensable steps in the conventional onshore seismic data processing, particularly in the western part of China; it is theoretically and practically significant to resolve the issue of static correction. Conventional refraction static correction is put forward under the assumption that layered medium is horizontal and evenly distributed. The complicated nature of the near surface from western part of China is far from the assumption. Therefore, the essential way to resolve the static correction problem from the complex area is to develop a new theory. In this paper, a high-precision non-linear first arrival tomography is applied to solve the problem, it moved beyond the conventional refraction algorithm based on the layered medium and can be used to modeling the complex near surface. Some of the new and creative work done is as follows: One. In the process of first arrival tomographic image modeling, a fast high-order step algorithm is used to calculate the travel time for first arrival and ray path and various factors concerning the fast step ray tracing algorithm is analyzed. Then the second-order and third-order differential format is applied to the step algorithm which greatly increased the calculation precision of the ray tracing and there is no constraint to the velocity distribution from the complex areas. This method has very strong adaptability and it can meet the needs of great velocity variation from the complicated areas. Based on the numerical calculation, a fast high-order step is a fast, non-conditional and stable high-precision tomographic modeling algorithm. Two, in the tomographic inversion, due to the uneven fold coverage and insufficient information, the inversion result is unstable and less reliable. In the paper, wavelet transform is applied to the tomographic inversion which has achieved a good result. Based on the result of the inversion from the real data, wavelet tomographic inversion has increased the reliability and stability of the inversion. Three. Apply the constrained high-precision wavelet tomographic image to the static correction processing from the complex area. During tomographic imaging, by using uphole survey, refraction shooting or other weathering layer method, weathering layer can be identified before the image. Because the group interval for the shot first arrival is relatively big, there is a lack of precision for the near surface inversion. In this paper, an inversion method of the layer constraint and well constraint is put forward, which can be used to compensate the shallow velocity of the inversion for the shot first arrival and increase the precision of the tomographic inversion. Key words: Tomography ，Fast marching method，Wavelet transform, Static corrections, First break

陆地三分量三维地震勘探技术研究与应用

With the development of oil and gas field exploration, it becomes harder to search new reserves. So a higher demand of seismic exploration comes up. Now 3C3D seismic exploration technology has been applied in petroleum exploration domains abroad. Comparing with the traditional P-wave exploration, the seismic attributes information which provided by 3C3D seismic exploration will increase quickly. And it can derive various combined parameters. The precision of information about lithology, porosity, fracture, oil-bearing properties, etc which estimated by above parameters was higher than that of pure P-wave exploration. These advantages mentioned above lead to fast development of 3C3D seismic technology recently. Therefore, how to apply the technology in petroleum exploration field in China, how to obtain high quality seismic data, and how to process and interpret real data, become frontier topics in geophysical field nowadays, which have important practical significance in research and application. In this paper, according to the propagation properties of P-wave and converted wave, a study of 3C3D acquisition parameters design method was completed. Main parameters included: trace interval, shot interval, maximum offset, bin size, the interval of receiving lines, the interval of shooting lines, migration aperture, maximum cross line distance, etc. Their determination principle was given. The type of 3C3D seismic exploration geometry was studied. By calculating bin attributes and analyzing parameters of geometry, some useful conclusions were drawn. With the method in this paper, real geometries for continental lithology stratum gas reservoir and fractured gas reservoir were studied and determined. In the static method of multi-wave, the near surface P-wave, S-wave parameter investigation method has been advanced, and this method has been applied for the patent successfully; the near surface P-wave, S-wave parameter investigation method and the converted refraction wave first arrival static techniques have been integrally used to improve the effectiveness of converted wave static. In the aspect of converted wave procession, the rotation of horizontal component data, the calculation of converted wave common conversion bin, the residual static of converted wave, the velocity analysis of the common conversion point (CCP), the Kirchhoff pre-stack time migration of converted wave techniques have been applied for setting up the various 3C3D seismic data processing flows based on different geologic targets, and the high quality P-wave, converted-wave profiles have been acquired in the actual data processing. In the aspect of P-wave and converted-wave comprehensive interpretation, the thoughts and methods of using zero-offset S-wave VSP data to calibrate horizon have been proposed; the method of using P-wave and S-wave amplitude ratio to predict the areas of oil and gas enrichment has been studied; the method of inversion using P-wave combined with S-wave has been studied; the various P-wave, S-wave parameters(velocity ratio, amplitude ratio, poisson ratio) have been used to predict the depth, physical properties, gas-bearing properties of reservoirs; the method of predicting the continental stratum lithology gas reservoir has been built. The above techniques have all been used in various 3D3C seismic exploration projects in China, and the better effects have been gotten. By using these techniques, the 3C3D seismic exploration level has been improved.

东海陆架盆地复杂天然气藏地震勘探判别技术研究

Exploration study proves that East sea shelf basin embeds abundant hydrocarbon resources. However, the exploration knowledge of this area is very low. Many problems in exploration are encountered here. One of them is that the gas reservoir of this area, with rapid lateral variation, is deeply buried. Correlation of Impendence between sandstone, gas sand and shale is very poor. Another problem is that the S/N ratio of the seismic data is very low and multiples are relatively productive which seriously affect reservoir identification. Resolution of the seismic data reflected from 2500-3000 meter is rather low, which seriously affects the application of hydrocarbon direct identification (HDI) technology. This research established a fine geological & geophysical model based on drilling、well logging、geology&seismic data of East sea Lishui area. A Q value extraction method from seismic data is proposed. With this method, Q value inversion from VSP data and seismic data is performed to determine the subsurface absorption of this area. Then wave propagation and absorption rule are in control. Field acquisition design can be directed. And at the same time, with the optimization of source system, the performance of high resolution seismic acquisition layout system is enhanced. So the firm foundation is ensured for east sea gas reservoir exploration. For solving the multiple and amplitude preserving problems during the seismic data processing, wave equation pre-stack amplitude preservation migration and wave equation feedback iteratively multiple attenuation technologies are developed. Amplitude preservation migration technology can preserve the amplitude of imaging condition and wave-field extrapolation. Multiple removing technology is independent of seismic source wavelet and velocity model, which avoiding the weakness of Delft method. Aiming at the complicated formation condition of the gas reservoir in this area, with dissecting typical hydrocarbon reservoir, a series of pertinent advanced gas reservoir seismic identification technologies such as petrophysical properties analyzing and seismic modeling technology、pre-stack/post-stack joint elastic inversion, attribute extraction technology based on seismic non-stationary signal theory and formation absorption characteristic and so on are studied and developed. Integrated analysis of pre-stack/post-stack seismic data, reservoir information, rock physics and attribute information is performed. And finally, a suit of gas reservoir identification technology is built, according to the geological and geophysical characteristics of this area. With developed innovative technologies, practical application and intergrated interpretation appraisal researches are carried out in Lishui 36-1.The validity of these technologies is tested and verified. Also the hydrocarbon charging possibility and position of those three east sea gas exploration targets are clearly pointed out.

多尺度分解与相位反褶积提高海上地震资料分辨率和信噪比

Offshore seismic exploration is full of high investment and risk. And there are many problems, such as multiple. The technology of high resolution and high S/N ratio on marine seismic data processing is becoming an important project. In this paper, the technology of multi-scale decomposition on both prestack and poststack seismic data based on wavelet and Hilbert-Huang transform and the theory of phase deconvolution is proposed by analysis of marine seismic exploration, investigation and study of literatures, and integration of current mainstream and emerging technology. Related algorithms are studied. The Pyramid algorithm of decomposition and reconstruction had been given by the Mallat algorithm of discrete wavelet transform In this paper, it is introduced into seismic data processing, the validity is shown by test with field data. The main idea of Hilbert-Huang transform is the empirical mode decomposition with which any complicated data set can be decomposed into a finite and often small number of intrinsic mode functions that admit well-behaved Hilbert transform. After the decomposition, a analytical signal is constructed by Hilbert transform, from which the instantaneous frequency and amplitude can be obtained. And then, Hilbert spectrum. This decomposition method is adaptive and highly efficient. Since the decomposition is based on the local characteristics of the time scale of data, it is applicable to nonlinear and non-stationary processes. The phenomenons of fitting overshoot and undershoot and end swings are analyzed in Hilbert-Huang transform. And these phenomenons are eliminated by effective method which is studied in the paper. The technology of multi-scale decomposition on both prestack and poststack seismic data can realize the amplitude preserved processing, enhance the seismic data resolution greatly, and overcome the problem that different frequency components can not restore amplitude properly uniformly in the conventional method. The method of phase deconvolution, which has overcome the minimum phase limitation in traditional deconvolution, approached the base fact well that the seismic wavelet is phase mixed in practical application. And a more reliable result will be given by this method. In the applied research, the high resolution relative amplitude preserved processing result has been obtained by careful analysis and research with the application of the methods mentioned above in seismic data processing in four different target areas of China Sea. Finally, a set of processing flow and method system was formed in the paper, which has been carried on in the application in the actual production process and has made the good progress and the huge economic benefit.

地震数据波动方程反演方法研究

Geophysical inversion is a theory that transforms the observation data into corresponding geophysical models. The goal of seismic inversion is not only wave velocity models, but also the fine structures and dynamic process of interior of the earth, expanding to more parameters such as density, aeolotropism, viscosity and so on. As is known to all, Inversion theory is divided to linear and non-linear inversion theories. In rencent 40 years linear inversion theory has formed into a complete and systematic theory and found extensive applications in practice. While there are still many urgent problems to be solved in non-linear inversion theory and practice. Based on wave equation, this dissertation has been mainly involved in the theoretical research of several non-linear inversion methods: waveform inversion, traveltime inversion and the joint inversion about two methods. The objective of gradient waveform inversion is to find a geologic model, thus synthetic seismograms generated by this geologic model are best fitted to observed seismograms. Contrasting with other inverse methods, waveform inversion uses all characteristics of waveform and has high resolution capacity. But waveform inversion is an interface by interface method. An artificial parameter limit should be provided in each inversion iteration. In addition, waveform information will tend to get stuck in local minima if the starting model is too far from the actual model. Based on velocity scanning in traditional seismic data processing, a layer-by-layer waveform inversion method is developed in this dissertation to deal with weaknesses of waveform inversion. Wave equation is used to calculate the traveltime and derivative (perturbation of traveltime with respect to velocity) in wave-equation traveltime inversion (WT). Unlike traditional ray-based travetime inversion, WT has many advantages. No ray tracing or traveltime picking and no high frequency assumption is necessary and good result can be got while starting model is far from real model. But, comparing with waveform inversion, WT has low resolution. Waveform inversion and WT have complementary advantages and similar algorithm, which proves that the joint inversion is a better inversion method. And another key point which this dissertation emphasizes is how to give fullest play to their complementary advantages on the premise of no increase of storage spaces and amount of calculation. Numerical tests are implemented to prove the feasibility of inversion methods mentioned above in this dissertation. Especially for gradient waveform inversion, field data are inversed. This field data are acquired by our group in Wali park and Shunyi district. Real data processing shows there are many problems for waveform inversion to deal with real data. The matching of synthetic seismograms with observed seismograms and noise cancellation are two primary problems. In conclusion, on the foundation of the former experiences, this dissertation has implemented waveform inversions on the basis of acoustic wave equation and elastic wave equation, traveltime inversion on the basis of acoustic wave equation and traditional combined waveform traveltime inversion. Besides the traditional analysis of inversion theory, there are two innovations: layer by layer inversion of seimic reflection data inversion and rapid method for acoustic wave-equation joint inversion.

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

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 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.

改善三维波动方程叠前深度偏移处理效果的方法研究

3D wave equation prestack depth migration is the effective tool for obtaining the exact imaging result of complex geology structures. It's a part of the 3D seismic data processing. 3D seismic data processing belongs to high dimension signal processing, and there are some difficult problems to do with. They are: How to process high dimension operators? How to improve the focusing? and how to construct the deconvolution operator? The realization of 3D wave equation prestack depth migration, not only realized the leap from poststack to prestack, but also provided the important means to solve the difficult problems in high dimension signal processing. In this thesis, I do a series research especially for the solve of the difficult problems around the 3D wave equation prestack depth migration and using it as a mean. So this thesis service for the realization of 3D wave equation prestack depth migration for one side and improve the migration effect for another side. This thesis expatiates in five departs. Summarizes the main contents as the follows: In the first part, I have completed the projection from 3D data point area to low dimension are using de big matrix transfer and trace rearrangement, and realized the liner processing of high dimension signal. Firstly, I present the mathematics expression of 3D seismic data and the mean according to physics, present the basic ideal of big matrix transfer and describe the realization of five transfer models for example. Secondly, I present the basic ideal and rules for the rearrange and parallel calculate of 3D traces, and give a example. In the conventional DMO focusing method, I recall the history of DM0 process firstly, give the fundamental of DMO process and derive the equation of DMO process and it's impulse response. I also prove the equivalence between DMO and prestack time migration, from the kinematic character of DMO. And derive the relationship between DMO base on wave equation and prestack time migration. Finally, I give the example of DMO process flow and synthetic data of theoretical models. In the wave equation prestak depth migration, I firstly recall the history of migration from time to depth, from poststack to prestack and from 2D to 3D. And conclude the main migration methods, point out their merit and shortcoming. Finally, I obtain the common image point sets using the decomposed migration program code.In the residual moveout, I firstly describe the Viterbi algorithm based on Markov process and compound decision theory and how to solve the shortest path problem using Viterbi algorithm. And based on this ideal, I realized the residual moveout of post 3D wave equation prestack depth migration. Finally, I give the example of residual moveout of real 3D seismic data. In the migration Green function, I firstly give the concept of migration Green function and the 2D Green function migration equation for the approximate of far field. Secondly, I prove the equivalence of wave equation depth extrapolation algorithms. And then I derive the equation of Green function migration. Finally, I present the response and migration result of Green function for point resource, analyze the effect of migration aperture to prestack migration result. This research is benefit for people to realize clearly the effect of migration aperture to migration result, and study on the Green function deconvolution to improve the focusing effect of migration.

电阻率测井反演方法应用研究

Formation resistivity is one of the most important parameters to be evaluated in the evaluation of reservoir. In order to acquire the true value of virginal formation, various types of resistivity logging tools have been developed. However, with the increment of the proved reserves, the thickness of interest pay zone is becoming thinner and thinner, especially in the terrestrial deposit oilfield, so that electrical logging tools, limited by the contradictory requirements of resolution and investigation depth of this kinds of tools, can not provide the true value of the formation resistivity. Therefore, resitivity inversion techniques have been popular in the determination of true formation resistivity based on the improving logging data from new tools. In geophysical inverse problems, non-unique solution is inevitable due to the noisy data and deficient measurement information. I address this problem in my dissertation from three aspects, data acquisition, data processing/inversion and applications of the results/ uncertainty evaluation of the non-unique solution. Some other problems in the traditional inversion methods such as slowness speed of the convergence and the initial-correlation results. Firstly, I deal with the uncertainties in the data to be processed. The combination of micro-spherically focused log (MSFL) and dual laterolog(DLL) is the standard program to determine formation resistivity. During the inversion, the readings of MSFL are regarded as the resistivity of invasion zone of the formation after being corrected. However, the errors can be as large as 30 percent due to mud cake influence even if the rugose borehole effects on the readings of MSFL can be ignored. Furthermore, there still are argues about whether the two logs can be quantitatively used to determine formation resisitivities due to the different measurement principles. Thus, anew type of laterolog tool is designed theoretically. The new tool can provide three curves with different investigation depths and the nearly same resolution. The resolution is about 0.4meter. Secondly, because the popular iterative inversion method based on the least-square estimation can not solve problems more than two parameters simultaneously and the new laterolog logging tool is not applied to practice, my work is focused on two parameters inversion (radius of the invasion and the resistivty of virgin information ) of traditional dual laterolog logging data. An unequal weighted damp factors- revised method is developed to instead of the parameter-revised techniques used in the traditional inversion method. In this new method, the parameter is revised not only dependency on the damp its self but also dependency on the difference between the measurement data and the fitting data in different layers. At least 2 iterative numbers are reduced than the older method, the computation cost of inversion is reduced. The damp least-squares inversion method is the realization of Tikhonov's tradeoff theory on the smooth solution and stability of inversion process. This method is realized through linearity of non-linear inversion problem which must lead to the dependency of solution on the initial value of parameters. Thus, severe debates on efficiency of this kinds of methods are getting popular with the developments of non-linear processing methods. The artificial neural net method is proposed in this dissertation. The database of tool's response to formation parameters is built through the modeling of the laterolog tool and then is used to training the neural nets. A unit model is put forward to simplify the dada space and an additional physical limitation is applied to optimize the net after the cross-validation method is done. Results show that the neural net inversion method could replace the traditional inversion method in a single formation and can be used a method to determine the initial value of the traditional method. No matter what method is developed, the non-uniqueness and uncertainties of the solution could be inevitable. Thus, it is wise to evaluate the non-uniqueness and uncertainties of the solution in the application of inversion results. Bayes theorem provides a way to solve such problems. This method is illustrately discussed in a single formation and achieve plausible results. In the end, the traditional least squares inversion method is used to process raw logging data, the calculated oil saturation increased 20 percent than that not be proceed compared to core analysis.