星间相干光通信中的光学锁相环

在相干体制下,二进制相移键控(BPSK)信号调制的零差接收机可实现理论上的最高灵敏度,是星间相干光通信的研究和应用重点。零差接收机要求本振波和信号波严格相位同步,常用的相位同步技术是光学锁相环(OPLL)。阐述了光学锁相环的基本原理,介绍了近20年来光学锁相环的发展,在此基础上总结了各种类型锁相环的性能和适用范围。回顾了欧洲国家在星间相干光通信计划中通信终端使用光学锁相环的情况。最后对光学锁相环技术进行了总结,对该技术的前景进行了展望。

空间激光跟瞄中指向驱动电机运动特性研究

设计了一种深空非合作目标的激光扫描、捕获、跟踪地面实验装置,通过模拟深空同轨道运动的两颗卫星跟瞄过程,在理论上计算了跟瞄装置中光束指向驱动电机的最小加速度和其在跟踪过程中的运动特性。理论分析与仿真结果表明,当卫星偏离光斑中心一定距离时,指向驱动电机先加速后减速,补偿这个偏心,重新捕获跟踪卫星;重新捕获到跟踪所需时间受电机加速度和望远镜探测精度以及探测器响应处理时间影响,其中探测器精度对重新捕获到跟踪所需时间影响较大,探测器响应处理时间要减小到最小;为了使从捕获到跟踪过程中卫星始终在扫描光斑范围内,经纬仪驱动电机的最小角加速度为25.5°/s2。

High-accuracy same-beam VLBI observations using Shanghai and Urumqi telescopes

The same-beam VLBI observations of Rstar and Vstar, which were two small satellites of Japanese lunar mission, SELENE, were successfully performed by using Shanghai and Urumqi 25-m telescopes. When the separation angle between Rstar and Vstar was less than 0.1 deg, the differential phase delay of the X-band signals between Rstar and Vstar on Shanghai-Urumqi baseline was obtained with a very small error of 0.15 mm rms, which was reduced by 1-2 order compared with the former VLBI results. When the separation angle was less than 0.56 deg, the differential phase delay of the S-band signals was also obtained with a very small error of several mm rms. The orbit determination for Rstar and Vstar was performed, and the accuracy was improved to a level of several meters by using VLBI and Doppler data. The high-accuracy same-beam differential VLBI technique is very useful in orbit determination for a spacecraft, and will be used in orbit determination for Mars missions of China Yinghuo-1 and Russia Phobos-grunt.

Bubble Dynamics and Heat Transfer in Microgravity Pool Boiling

Boiling is an extremely complicated and illusive process. Microgravity experiments offer a unique opportunity to study the complex interactions without external forces, such as buoyancy, which can affect the bubble dynamics and the related heat transfer. Furthermore, they can also provide a means to study the actual influence of gravity on the boiling. Two research projects on pool boiling in microgravity have been conducted aboard the Chinese recoverable satellites. Ground-based experiments both in normal gravity and in short-term microgravity in the Drop Tower Beijing and numerical simulations have also been performed. Steady boiling of R113 on thin platinum wires was studied with a temperature-controlled heating method, while quasi-steady boiling of FC-72 on a plane plate was investigated with an exponentially increasing heating voltage. It was found that the bubble dynamics in microgravity has a distinct difference from that in normal gravity, and that the heat transfer characteristic is depended upon the bubble dynamics. Lateral motions of bubbles on the heaters were observed before their departure in microgravity. The surface oscillation of the merged bubbles due to lateral coalescence between adjacent bubbles drove it to detach from the heaters. Slight enhancement of heat transfer on wires is observed in microgravity, while diminution is evident for high heat flux in the plate case.

Photoionization of 1s electron and corresponding shake-up process in ground and excited lithium atoms

The cross sections of the 18 electron photoionization and corresponding shake-up processes for Li atoms in the ground state 1s(2)2s and excited states 1s(2)2p, 1s(2)3p, 1s(2)3p and 1s(2)3d are calculated using the multi-configuration Dirac-Fock method. The latest experimental photoelectron spectrum at hv = 100 eV [Cubaynes D et al. Phys. Rev. Lett. 99 (2007) 213004] has been reproduced by the present theoretical investigation excellently. The relative intensity of the shake-up satellites shows that the effects of correlation and relaxation become more important for the higher excited states of the lithium atom, which are explained very well by the spatial overlap of the initial and final state wavefunctions. In addition, strong dependence of the cross section on the atomic orbitals of the valence electrons are found, especially near the threshold.

卫星表面光学反射特性建模与仿真研究

在研究空间相机成像效果探测问题时,针对提高探测相机成像真实效果的需要,研究了卫星表面光学反射特性的建模与仿真方法。利用双向反射分布函数(BRDF)仿真卫星表面的光学反射特性,将卫星表面物理特性与其光学反射特性有机地结合起来。在VC.Net环境下开发卫星表面光学反射特性仿真软件,并将软件集成到空间探测相机成像效果仿真平台上。利用平台进行成像仿真实验,生成真实感仿真图片。仿真结果表明,BRDF能够精确仿真卫星表面的光学反射特性,使成像效果仿真更加准确和逼真。

COMPARATIVE-STUDIES ON PERIODATOCUPRATE(II, III) AND TELLURATOCUPRATE(II, III)

The crystal structures, electronic spectra, and Cu2p XPS of Cu(III) complexes Na4H[Cu(H2TeO6)(2)]. 17H(2)O and Na4K[Cu(HlO(6))(2)]. 12H(2)O have been described. The characterizations of a Cu(III) atom in a complex are as follows: (i) In a square-planar coordination, the average bond length of Cu-O is 0.183 nm, shorter than the 0.190-0.200 nm found for a Cu(II) complex. (2) The ''blue shift'' occurs for d-d transitions in the electronic spectrum of the Cu(III) complex compared to those of its related Cu(II) complex, resulting from the higher valence state. (3) Cu(III) compounds with CuO4 square-planar coordination are expected to be diamagnetic whereas Cu(II) compounds to be paramagnetic. (4) Comprehensive investigations on Cu2p XPS show that the binding energy of Cu2p(3/2) of a pure Cu(III) compound is about 2.0 eV higher than that of its corresponding Cu(II) compound: the shake-up satellites do not appear in the Cu2p XPS for a pure diamagnetic Cu(III) compound, the same as found for a diamagnetic Ni(II) compound: the FWHM of the signal of Cu2p XPS may become broader for Cu(III) compound because its core hole's lifetime shortens due to the higher valence state of copper. (C) 1995 Academic Press, Inc.

X-RAY PHOTOELECTRON-SPECTRA OF SOME BIOINORGANIC COMPLEXES OF THE RARE-EARTHS WITH N-ACETYLALANINE

X-Ray photoelectron spectra of some bioinorganic complexes of La, Ce, PT, Nd, Sm and Eu with N-acetylalanine have been measured and the 3d5/2 and 3d3/2 main peaks and their satellites have also been assigned. ne spin-orbit splitting between the 3d5/2 and 3d3/2 core-level of the rare earth ion in these complexes becomes slightly larger than that of the free rare earth atom due to the effect of the crystal field. The satellite for the 3d main peaks of La in the solid state complex are in higher binding energy region and may be attributable to the L --> 4f charge-transfer shake-up process. The satellites for the 3d main peaks of Ce, Pr, Nd, Sm and Eu are in the lower binding energy region and may be attributable to the 4f --> L charge-transfer shake-down process.

Near-Space Wide-Swath Radar Imaging With Multiaperture Antenna

Near-space, defined as the altitude region between 20 and 100 km, offers many capabilities that are not accessible for low Earth-orbit (LEO) satellites or airplanes because it is above storm and not constrained by orbital mechanics and high fuel consumption. Hence, a high flying speed can be obtained for the maneuvering vehicles operating in near-space. This offers a promising solution to simultaneous high-resolution and wide-swath synthetic aperture radar (SAR) imaging. As such, one near-space wide-swath SAR imaging technique is presented in this letter. The system configuration, signal model, and imaging scheme are described. An example near-space SAR system is designed, and its imaging performance is analyzed. Simulation results show that near-space maneuvering vehicle SAR indeed seems to be a promising solution to wide-swath SAR imaging.

Vertical structure and evolution of the Luzon Warm Eddy

Eddies are frequently observed in the northeastern South China Sea (SCS). However, there have been few studies on vertical structure and temporal-spatial evolution of these eddies. We analyzed the seasonal Luzon Warm Eddy (LWE) based on Argo float data and the merged data products of satellite altimeters of Topex/Poseidon, Jason-1 and European Research Satellites. The analysis shows that the LWE extends vertically to more than 500 m water depth, with a higher temperature anomaly of 5A degrees C and lower salinity anomaly of 0.5 near the thermocline. The current speeds of the LWE are stronger in its uppermost 200 m, with a maximum speed of 0.6 m/s. Sometimes the LWE incorporates mixed waters from the Kuroshio Current and the SCS, and thus has higher thermohaline characteristics than local marine waters. Time series of eddy kinematic parameters show that the radii and shape of the LWE vary during propagation, and its eddy kinetic energy follows a normal distribution. In addition, we used the empirical orthogonal function (EOF) here to analyze seasonal characteristics of the LWE. The results suggest that the LWE generally forms in July, intensifies in August and September, separates from the coast of Luzon in October and propagates westward, and weakens in December and disappears in February. The LWE's westward migration is approximately along 19A degrees N latitude from northwest of Luzon to southeast of Hainan, with a mean speed of 6.6 cm/s.

Anti-cyclonic eddies northwest of Luzon in summer-fall observed by satellite altimeters

Anti-cyclonic eddies northwest of Luzon of the Philippines in summer-fall are identified in the merged data products of satellite altimeters of Topex/Poseidon, Jason-1 and European Research Satellites. The generation and propagation of the anti-cyclonic eddies, which are confirmed by satellite ocean color data, are found to be a seasonal phenomenon that is phase-locked to the onset of the southwesterly monsoon and the relaxation of the cyclonic wind curl in the northeastern South China Sea. The eddies originate from northwest of Luzon in summer, move across the northeastern South China Sea to reach the China continental slope in fall, and propagate southwestward along the continental slope in fall-winter, inducing shelfbreak current variations in the western South China Sea in fall-winter. The anti-cyclonic eddy discovered by Li et al. (1998) in the northern South China Sea is found to originate from northwest of Luzon and carry primarily the South China Sea waters. It does not appear to be an eddy shed from the Kuroshio in the Luzon Strait area as alluded by Li et al. (1998) and others.

多源多时相SAR资料反演水下地形的同化模型

合成孔径雷达(SAR)是重要的微波传感器之一，在海洋遥感领域得到了广泛应用。浅海水下地形SAR探测是SAR海洋应用的重要部分之一，发展业务化水下地形SAR探测技术是水下地形SAR探测研究的最终目标。现有的水下地形SAR探测是基于单景SAR影像开展的，探测结果不能全面反映水下地形。因此，发展一种基于多源多时相SAR影像的水下地形反演方法是该研究方向的前沿问题。 本文基于水下地形SAR成像机制和袁业立（1997）导出的海波高频谱解析表达式，首次建立了基于多源多时相SAR影像的水下地形反演的同化模型与算法，并通过该模型在台湾浅滩水下地形SAR探测中的应用检验了该模型和算法的可行性。本文还开展了已有的浅海水下地形SAR探测技术（浅海水下地形单景SAR影像反演）在台湾浅滩水下地形探测中的应用研究，以及反演初始水深对水下地形SAR反演结果影响的研究。此外，本文还对中国近海水下地形SAR成像能力进行了分析。 在浅海水下地形SAR探测技术应用于台湾浅滩研究方面，基于7景SAR影像开展了台湾浅滩单景SAR影像反演计算，反演结果表明同一水下地形不同时相SAR影像反演结果不能全面反映真实地形特征。 在反演初始水深对水下地形SAR反演结果影响研究方面，通过不同初始水深条件下台湾浅滩水下地形SAR反演结果的比较分析，可以得出：初始水深影响水下地形反演结果的准确性，初始水深的选取需要较粗的实际水深的支持； 在多源多时相水下地形SAR反演同化模型研究方面， 建立了多源多时相水下地形SAR反演的同化模型与算法，通过基于多源多时相SAR影像的台湾浅滩五种不同情形水下地形反演计算，证明了多源多时相水下地形SAR反演同化模型和算法是可行的，并给出了基于多源多时相SAR影像的水下地形反演计算中SAR影像的选取原则。 在中国近海浅海水下地形SAR响应能力分析方面，基于对中国近海潮汐潮流状况、水下地形概况和已有中国近海SAR影像中水下地形特征等的分析，总结出中国近海适合开展浅海水下地形SAR探测的区域主要有：渤海海域、苏北近海海域、台湾海峡海域和南沙群岛海域。

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

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.

磁层-电离层耦合现象的分析研究

The magnetosphere-ionosphere coupling is mainly manifested by the trans- porting processes of energy into the ionosphere , the energy is carried by solar wind and firstly accumulate at the magnetosphere, and the coupling processes also significantly include the interaction between the magnetosphere and ionosphere for mass and energy. At the quiet condition, energy is delivered by the large-scale convection of the geomagnetic field; the huge energy from solar wind bulk will be injected into and consumed at the near magnetosphere and ionosphere by the geomagnetic storm and substorm activities. Aurorae and FACs (Field-aligned currents) are the important phenomena in the coupling processes. In the present work, firstly, we analyze the activity characteristics of auroral precipitating particle, secondly, we study the distribution characters of large-scale field aligned currents (LS FACs) at storm-time using the observations from different satellites at different altitudes. Finally, we investigate the evolution of the geomagnetic field configuration at the nightside sector on the onset of the expansion phase in a substorm event, the substorm event happened at 0430UT to 0630UT on 8th Nov. 2004. The main results as follows： At the first, the data of the estimated power input (EPI) of auroral particles from NOAA/POES (Polar orbiting environmental satellite) for some 30 years have been analyzed. The variation tendencies of the EPI generally coincide with aa, AE and Dst indices. The annual variation of EPI shows equinox peaks and an asymmetric-activity with a higher peak in the winter-hemisphere than in the summer-hemisphere. The diurnal UT variations are different from north and south hemisphere: for north hemisphere, the peak appears at 1200UT, and the relative deviation is 22% to the daily average of the north hemisphere. For south hemisphere, the maximal deviation is 22% at 2000UT. So the diurnal variation of EPI is more dominant than the annual variation which maximal deviation is 3% to 12% for different seasons. Studies on correlations of the hourly average of EPI, Pa, with AE and Dst indices show a correlation coefficient r=0.74 of Pa and AE, and r=-0.55 of Pa and Dst. The hourly EPIs for north and south polar regions, NPa and SPa, show a north-south asymmetry with a higher correlation of SPa and AE (or Dst). Time delays of EPI with respect to magnetic indices are examined, the maximum correlation coefficient of Pa with AE (r=0.78) occurs when the time delay =0, suggesting a synchronous activity of auroral electrojet and auroral precipitating particles, while =1-2h, the correlation coefficient of Pa with Dst is maximum (r=0.57), suggesting that the activity of auroral particle precipitating may influence the ring current on some extent. Sencondly, we use the high-resolution magnetic field vector data of the CHAMP satellite to investigate the distribution of large-scale FACs during the great magnetic storm on 7th to 8th Nov. 2004. The results show that, whether in the northern or southern hemisphere, the number and density of large-scale FACs during the main-phase are more and bigger than these during the recover-phase, and the number of large-scale FACs in morning sector obviously is more than that in afternoon sector. In terms of the magnetic indices, we find that large-scale FACs in morning sector significantly affected by the substorm activities, while in afternoon sector the large-scale FACs mainly indicate the fluctuations of the ring-current in storm time. Accordingly to the former studies, similarly, we find that in the morning sector, the scale of the large-scale FACs move to the high-latitude region, and in the afternoon sector, large-scale FACs distinctly expand to the low-latitude region. During the time periods that the NOAA/POES auroral precipitating particle power data temporally correspond to the large-scale FACs, the more the power of auroral particle is, the more and bigger the number and density of FACs are. At the same time, we use the magnetic field vector data of POLAR obtain a good form of region 1, region 2, and three pieces of cusp FACs during a single transit at 1930UT-2006UT on 07th. And the characteristics of simultaneous electric field and energy particles observations on Polar are coincide with the five FACs pieces. Finally, by means of the observation of Cluster 4 and Goes 10、 Goes 12, we analyze the evolution process of the change of the magnetic field configuration at night sector at the expansion phase of a substorm event which happened during 0430UT to 0630UT on 8th Nov. 2004, we find that the times of the beginning of the polarizations of magnetic field are observed from Goes 10 to Goes 12 then to Cluster 4. So, at the synchronous orbit ( 6.6 RE) to 10RE distance scale of the neutral sheet, the current disruption spread tailward. Simultaneously, the strengthen of the FACs deduced from these satellites’ magnetic field observations are almost consistent with the times of polarizations, as well as the high energy particles injection and the electric field dominant variation. The onset times determined by the magnetic field polarizations from these satellites are all ahead of the onset time that confirmed from the auroral electrojet indices. So, these characters of different observations can be used as the criterions to determine the onset time for the substorms of such type as we studied.

沾化凹陷泥质岩油气藏成藏特征研究

Mudstone reservoir is a subtle reservoir with extremely inhomogeneous, whose formation is greatly related to the existence of fracture. For this kind of reservoir, mudstone is oil source rock, cover rock and reservoir strata, reservoir type is various, attitude of oil layer changes greatly, and the distribution of oil and gas is different from igneous or clastic rock reservoir as well as from carbonate reservoir of self-producing and self-containing of oil and gas. No mature experience has been obtained in the description, exploration and development of the reservoir by far. Taking Zhanhua depression as an example, we studied in this thesis the tectonic evolution, deposit characteristics, diagenesis, hydrocarbon formation, abnormal formation pressure, forming of fissure in mudstone reservoir, etc. on the basis of core analysis, physical simulation, numerical simulation, integrated study of well logging and geophysical data, and systematically analyzed the developing and distributing of mudstone fissure reservoir and set up a geological model for the formation of mudstone fissure reservoir, and predicted possible fractural zone in studied area. Mudstone reservoir mainly distributed on the thrown side of sedimentary fault along the sloping area of the petroleum generatiion depression in Zhanhua depression. Growing fault controlled subsidence and sedimentation. Both the rate of subsidence and thickness of mudstone are great on the thrown side of growing fault, which result in the formation of surpressure in the area. The unlocking of fault which leads to the pressure discharges and the upward conduct of below stratum, also makes for the surpressure in mudstone. In Zhanhua depression, mudstone reservior mainly developed in sub-compacted stratum in the third segment of Shahejie formation, which is the best oil source rock because of its wide spread in distribution, great in thickness, and rich in organic matter, and rock types of which are oil source mudstone and shale of deep water or semi-deep water sediment in lacustrine facies. It revealed from core analysis that the stratum is rich in limestone, and consists of lamina of dark mudstone and that of light grey limestone alternately, such rock assemblage is in favor of high pressure and fracture in the process of hydrocarbon generation. Fracture of mudstone in the third segment of Shahejie formation was divided into structure fracture, hydrocarbon generation fracture and compound fracture and six secondary types of fracture for the fist time according to the cause of their formation in the thesis. Structural fracture is formed by tectonic movement such as fold or fault, which develops mainly near the faults, especially in the protrude area and the edge of faults, such fracture has obvious directivity, and tend to have more width and extension in length and obvious direction, and was developed periodically, discontinuously in time and successively as the result of multi-tectonic movement in studied area. Hydrocarbon generation fracture was formed in the process of hydrocarbon generation, the fracture is numerous in number and extensively in distribution, but the scale of it is always small and belongs to microfracture. The compound fracture is the result of both tectonic movement and hydrocarbon forming process. The combination of above fractures in time and space forms the three dimension reservoir space network of mudstone, which satellites with abnormal pressure zone in plane distribution and relates to sedimentary faces, rock combination, organic content, structural evolution, and high pressure, etc.. In Zhanhua depression, the mudstone of third segment in shahejie formation corresponds with a set of seismic reflection with better continuous. When mudstone containing oil and gas of abnormal high pressure, the seismic waveform would change as a result of absorb of oil and gas to the high-frequency composition of seismic reflection, and decrease of seismic reflection frequency resulted from the breakage of mudstone structure. The author solved the problem of mudstone reservoir predicting to some degree through the use of coherent data analysis in Zhanhua depression. Numerical modeling of basin has been used to simulate the ancient liquid pressure field in Zhanhua depression, to quantitative analysis the main controlling factor (such as uncompaction, tectonic movement, hydrocarbon generation) to surpressure in mudstone. Combined with factual geologic information and references, we analyzed the characteristic of basin evolution and factors influence the pressure field, and employed numerical modeling of liquid pressure evolution in 1-D and 2-D section, modeled and analyzed the forming and evolution of pressure in plane for main position in different periods, and made a conclusion that the main factors for surpressure in studied area are tectonic movement, uncompaction and hydrocarbon generation process. In Zhanhua depression, the valid fracture zone in mudstone was mainly formed in the last stage of Dongying movement, the mudstone in the third segment of Shahejie formation turn into fastigium for oil generation and migration in Guantao stage, and oil and gas were preserved since the end of the stage. Tectonic movement was weak after oil and gas to be preserved, and such made for the preserve of oil and gas. The forming of fractured mudstone reservoir can be divided into four different stages, i.e. deposition of muddy oil source rock, draining off water by compacting to producing hydrocarbon, forming of valid fracture and collecting of oil, forming of fracture reservoir. Combined with other regional geologic information, we predicted four prior mudstone fracture reservoirs, which measured 18km2 in area and 1200 X 104t in geological reserves.