10月 282013
 

多年来,经全球学者的广泛研究。已经获知莫霍面是一个全球性的地球内部分界面。莫霍面及其分布特征反映了地壳的演化和壳幔的相互作用。对应不同的地球动力学过程和构造单元,莫霍面性质与分布特征有所不同。稳定地块的莫霍面呈一级速度间断面,分布较稳定,通常具有较强的地震反射;活动造山带的莫霍面常呈现出速度梯度变化,分布不稳定,地震反射通常较弱,有时缺失。莫霍面作为建立地壳和岩石圈结构三维几何模型的重要组成部分,是了解区域构造演化的关键。莫霍面深度(通常表征地壳厚度)记录了地壳的生长与经历的地球动力学过程,因而莫霍面深度是人们认识地球演化的关键证据。由于地震波在莫霍面上传播特征和产生的强弱反射最为明显,能够为人们直观识别,因此地震学探测(主动源和被动源)一直是探知莫霍面深度的主要方法。青藏高原,作为地球上海拔最高的大陆地区,地壳最厚达80—90 km。其巨厚的地壳厚度一直是研究青藏高原深部结构和大陆动力学的重要参数,为全球地学界所关注。自20世纪50年代至今,大量被动源地震观测和主动源地震探测资料的获得揭示了青藏高原的地壳与上地幔结构。

 Posted by at 上午 8:01
10月 282013
 

深地震测深(Deep Seismic Sounding),也称之为宽角反射与折射(wide.angle reflection and refraction),是用于研究深达几十公里的地壳上地幔结构的一种有效地震探测方法。目前来说,特别对确定莫霍面深度是比较有效的方法。在20世纪60年代的国际上地幔计划、70年代的地球动力学计划以及80年代的岩石罔计划的实施过程中,深地震测深研究在全球获得了大量的成果,几乎遍布所有构造单元.比较精细刻画了地球内部结构,反映出地壳上地幔复杂的结构及不均匀性。

“深地震反射剖面”(Deep Seismic Reflection Profiling)是研究岩石圈精细结构最有效的方法,能获取地壳内部和莫霍面的细微变化。青藏高原的深地震反射剖面探测可追溯到1958年柴达木盆地的反射地震法地壳调查(曾融生等,1961),首次揭示了柴达木盆地的地壳厚度。INDEPTH项目首次将多次覆盖的深反射地震方法应用于青藏高原的探测研究,发现喜马拉雅山北坡莫霍面深度可达72.76km(Zhao,et a1.,1993)。1992年,“格尔木.额济纳旗”地学断面项目组在青藏高原北缘的祁连山进行了深地震反射剖面探测试验,揭示出横过盆山结合部位莫霍面由北面河西走廊44 km向祁连山下加深到50km(吴宣志等,1995)。1998.2000年中国地质科学院岩石圈研究中心在跨越西昆仑.塔里木、阿尔金山.塔里木盆山结合带完成2条深地震反射剖面探测。发现莫霍面的强反射特征,为研究青藏高原北缘碰撞变形的深部过程提供了有力证据(高锐等,2000,2001a,2001b)

宽频地震观测
近年来,随着数字化观测技术的迅速发展,宽频地震观测(Broadband Seismic Observation)方法被广泛采用,研究地震观测台阵下方地壳和上地幔结构,结果层出不穷,特别是接收函数逐渐成为研究地壳上地幔结构的一种简易高效的方法。青藏高原宽频地震观测始于1991年,曾融生与吴大铭的中美合作计划中,沿格尔木到日喀则一线布置了11台宽频带地震仪,揭示了剖面沿线的莫霍面深度变化(曾融生等,1992;吴庆举等,1998)。1992.1994年,中法合作进行了青藏高原天然地震合作研究,沿着青藏高原主干公路布设了宽频地震仪(姜枚,1994)。1995.1998年,又完成了青藏高原北缘的观MY(姜枚等,1999)。INDEPTH.III“德庆.龙尾错”宽频地震观测穿越了班公湖.怒江缝合带,并进入了羌塘地体,对高原内部深部构造特征进行了研究(Tilmann,et a1.,2003;吴庆举等,2004)。1997.1998年,大陆学者与中国台湾学者合作进行塔里木.西昆仑北部的宽频地震观测,发现西昆仑下莫霍面的叠置现象(Kao,et a1.,2001)。进入21世纪后,在青藏高原又进行了大量的宽频地震观测,其中包含了穿越喜马拉雅的“Nepal—Tibet”计划和“Hi—Climb”项目(Pelkum S.。et a1.,2005;钱辉等,2007;李海鸥等,2008;姜枚等,2008),以及穿越青藏高原东缘龙门山断裂带的“川西高原”和“林芝.永川”观测(xu,eta1.,2007;王椿镛,2008)。目前,横跨青藏高原各主要地体都实施了宽频带地震观测(图1),累积长度达10000余公里。

 Posted by at 上午 7:36
10月 262013
 

1.什么是地震勘探,它主要包括哪几个阶段?
答:1)由人工激发的地震波(弹性波),穿过地下介质运动、遇到弹性分界面返回地面,用仪器接收地震波,得到地震记录。对接收到的地震记录进行处理、解释,从而就能了解地下介质的情况,这个过程叫地震勘探
2)它由三个部分组成:野外资料采集、室内资料处理、室内资料解释。

2.什么是弹性介质?什么是均匀介质?
答:弹性介质:物体受力后,发生形变,但当外力撤消后,即能恢复原状的性质。
均匀介质:波速是常数(定值)的介质
附加:(塑性介质:物体受力后,发生形变,但当外力撤消后,不能恢复原状的性质。)

3.什么是杨氏模量、剪切模量及泊松比的物理意义?
答:(1) 杨氏模量(E):简单拉伸或压缩时,弹性体的相对伸缩△L/L与应力P之比 E=P/(△L/L) ;不同的物体E是不同的,在线性弹性极限范围内,物体的弹性形变满足虎克定律(应力∝应变)
(2) 切变模量(剪切模量)μ:它是简单剪切力作用时的切应力P 与剪应变tgθ之值,即有 μ=P/tgθ=P/( △L/L)
(3) 泊松比(σ):弹性体内发生纵向伸长(或缩短)时,伴随产生的横向相对收缩(或膨胀) △d/d与纵向相对伸(缩) △L/L之比值,称泊松比.
σ=(△d/d)/( △L/L) 它是表示形变变化调整的一种尺度.

4.试叙纵波和横波的传播特点。
答:纵波是弹性介质发生体积形变所产生的波动(体积变化),是一种胀缩力形成的波质点的振动(位移)与波的传播方向一致,可在任何介质中传播
横波是弹性介质发生切变时所产生的波动(形状变化),是旋转力作用形成的波 ,质点的振动(位移)与波传播的方向垂直,只在弹性固体中传播,即横波不通过液体、气体,因为剪切模量=0,纵波速度比横波速度大,所以远离震源时总是纵波先到达检波器 。

5.什么叫射线?
答:射线就是波从一点到另一点传播的路径,它代表了波传播的方向。
射线永远垂直于波前。

6.振动图和拨剖面有什么区别?
答:(1)振动图:固定一点(X=X1)→U=U(t) →振动图
描述参数:视周期T* 视振幅A* 初至t1 延续度△t
(2)波剖面:固定某一时刻(t=t1)→U=U(X)→波剖面
参数描述: 波峰:波剖面中最大正位移;波谷:波剖面中最
负位移;视波长λ*两个相邻波峰或波谷的距离,它表示波在一个视周期这传播的距离。

7.从反射和折射波形成的机制,分析反射波和折射波形成的条件是什么?
答:1)形成反射波的条件:上、下介质界面必须是一个波阻抗界面,即波阻抗差不为零。
2)形成折射波的条件:1》下面介质的波速要大于所有上面介质的波速
2》入射角是以临界角I 入射,达到临界角i=arcsin(V1/V2).

8. 影响地震波速度的主要地质因素有哪些?
答:1.岩石密度、地质年代对地震波速度的影响;(成正比);
2.地层的埋藏深度对速度的影响 ;(成正比);
3.岩石的孔隙度对速度的影响;(孔隙度大,则速度就小);
4.岩石中的孔隙充填物对速度的影响
9.什么叫时距曲线?

答:1》表示地震波的传播时间和爆炸点与检波器之间的距离X的关系曲线,即T—X曲线,称为时距曲线。
2》共炮点时距曲线:由一点激发,若干接收点接受收
10.水平共面炮点反射波时距波与工中心点反射波时距曲线有何异同?
答:两者曲线形式一样,但物理含义不同
共反射点 to含义不同 动校正含义不同
CSP 一段界面 炮点处H回声时间 各道反射时间与炮点处to
时间之差
CRP 一个反射点 M点处回声时点 各道反射时间与M点to之差

11.什么叫视速度?它有何用途?
答:所谓视速度即沿任意方向测波前时所测速度。
用途:用来定性判断时距曲线的弯曲程度。(判断原则:视速度大ˉ斜率小ˉ曲线平缓 ;视速度小ˉ斜率大ˉ曲线弯曲)

12.平均速度是怎么表示的?平均速度和均方根速度的定义是什么?
答:平均速度:波垂直穿过地层的总厚度与总的传播时间之比(在课件中注意)
均方根速度:把层状介质的波的高次曲线看成是二次曲线,此时波所具有的速度叫均方根速度。

13.地震资料采集的工序有哪些?
答:踏勘工区,布置测线,试验工作

14.什么叫观测系统?
答:激发点与接收点在测线上的相对位置关系—观测系统

15.什么是地震地址条件?好的条件有哪些?
答:地震地质条件分为:1)浅层地震地质条件; 2)深层地震地质条件。
一般好的浅部地震地质条件是;
1.低速带 (风化层)厚度小,横向稳定
2.潜水面浅,含水丰富,岩层塑性大(泥岩)能激发出较强能量、频率适中的地震波。
3.主要地震界面上次一级地震界面少,可减少界面的透射损失。
4.剖面上无高速层
好的深部地震地质条件是:
1.地震层位与地质层位一致
2.具有较好的地震标准层。

16.陆地上和海洋上地震勘探的异同什么?
答:相同点:原理相同
不同点:1.)装备不同;2)海洋多采用空气枪、点火花,陆地上多采用炸药震源、可控震源;3)检波器不同,海洋多采用漂浮检波器

17.什么是多次覆盖?
答:所谓多次覆盖(Multifold):即对地下同一反射点,进行重复多次观测,目的是突出反射波,压制干扰波,提高信噪比。

18.水平界面情况下,一次反射波和多次反射波叠加效应是怎样的?
答:一次反射波效应增强,多次反射波叠加效应减弱。

19.共反射点(共中心点)叠加法对哪些波有作用?
答:共反射点(共中心点)叠加突出一次反射波,压制多次反射波。
(多次覆盖对于一次反射波来说,相当于不同位置相同时间波的同相叠加,叠加后能量增强。 这就回答了为什么多次覆盖能突出一次反射波,提高资料信噪比。)

20.为什么要进行数据处理?
答:野外地震资料中包含着有关地下构造和岩性的信息,但这些信息是叠加在干扰背景上且被一些外界因素所扭曲,信息之间往往是互相交织的,不宜直接用于地质解释。因此,需要对野外采集的地震资料进行室内处理。

21.地震资料结实的主要内容是什么?
答:构造的解释、岩性的解释、储存研究。

22.地震分辨率是指什么?
答:地震分辨率分为垂直分辨率和水平分辨率
垂直分辨率指在纵向上能分辨岩层的最小厚度;
横向分辨率指在横向上确定地质体(如断层点、类灭点)位置和边界的精确程度。

23.地震同相轴(波)对比的原则是什么?
答:同相性、振幅显著增强、波形相似、时差变化规律。

24.地震剖面上产生假象的原因有哪些?
答:1)与速度有关的假象; 2)与几何因素有关的假象;
3)处理参数不当产生的假象; 4)静校正参数不当产生的假象.

25地震资料构造解释的主要步骤是什么?
答:1)掌握地质资料,统观全局,研究批剖面的结构;2)重点研究标准层的反射同相轴;3)相位对比;4)波组对比;5)研究异常波;剖面间的对比。

26.地层剖面上断层的主要表现特征有哪些?
答:1)同相轴错断;
2)反射波同相轴数目突然增减或消失,波组间隔变化;
3)反射波同相轴形状突变,出现反射零乱或空白带;
4)同相轴分叉,合并,扭曲,强相位转换等
5)异常波的出现如绕射波,断面反射波等。

27.断点在平面图上的组合应遵循的一般规律有哪些?
答: 1)先主后次;2)先简单后复杂;3)同一断层在平行的时间剖面上性质相同;4)同一断块内,地层产状的变化应有规律。 5)断层两侧波组具明显特征。6)断点组合要遵循断裂力学机制的规律;7)要尽可能弄清控制断层的构造性质和其成因类型。 8).断点的组合有认识-修改-再认识的过程。

28.纵向分辨率和横向分辨率的含义是什么?
答:纵向分辨率指在纵向上能分辨岩层的最小厚度;
横向分辨率指在横向上确定地质体(如断层点、类灭点)位置和边界的精确程度。

29.高分辨率工作中何为“四高四小”工作方法?
答:所谓“四高”即高时间采样率、高宽频道接收、高次覆盖、高检波带频率;
“四小”即小道距、小偏移距、小组合激距、小药量激发。

30.VSP的概念是什么?
答:在地表设置震源激发地震波,在井内安装检波器接收地震波,即在垂直方向上观测一维人工场,然后对所观测到的资料经过校正、叠加、滤波等处理,得到垂直地震剖面。

31.VSP资料有哪些用途?
答:1)求准速度;2)标定地震地质层位;3)识别多次波;4)提出反褶积因子;5)预测井底下反射层的深度;6)计算吸收衰减系数;7)计算反射界面的倾角;8)提取纵、横波的速度比及泊松比等参考数,岩性解释。

32.三维地震资料与二维地震资料构造解释的异同。
答:1)二维与三维的解释流程和解释内容基本一致
2)剖面的显示三维比二维丰富(信息丰富)
3)不能按二维剖面解释的方法对三维资料作解释,对数据加以利用
4)可利用三维可视化工具对三为地震资料进行解释验证

 Posted by at 下午 1:43
10月 262013
 

介绍

基于python的地球物理模型与解释程序,正在开发中。

详情

The PyGMI software is designed to be modular. Independant groups of modules are stored in separate directories under the packages directory.

To add your own modules:

  1. add your own subdirectory under the packages directory
  2. create a menu.py file in this directory.You may want to copy an existing one from another directory as a template
  3. create an iodefs.py file in this directory if you wish to program a custom import.

https://code.google.com/p/pygmi/wiki/Developers

 Posted by at 上午 1:44
10月 242013
 

xyzgms 1.9.1 geosoft

XYZGMS 生成所有模型,模型中的第一个点必须放置在必须的文件X=0 (GM-SYS 模型坐标),GM-SYS模型坐标将计算到第一个点的距离,-x0参数可重新指定第一个点的值,这个坐标也可以用GM-SYS的工具软件CVTGMS 来修改。

XYZGMS 计算坐标变换,用于实际坐标系下模型的定位与输出。如果模型用gmsmark, gmsprof, or gms GXs生成,并且数据指定有投影信息,则这些信息也存储在模型中,坐标信息存储在.ECS文件中。

XYZGMS计算剖面的方位角和指定相对走向angle=90, 磁倾角magnetic inclination=60, 偏角declination=20 和场强field intensity=56,000 作为默认值. 你可以在GM-SYS中为每个模型改为剖面比较合适的值.

输入的最大列数: 256

最大组数(data fields): 50

模型最大观测点数:16,000

 

geosoft.xyz文件创建一个GM-SYS模型

用法:XYZGMS in out -l=或 -f  -zg= 或 -zm=  [可选参数]

其中:in 基点文件名.xyz

out 输出文件名

-l= 线的类型[B,T,L]和线号

-f 输入的是一个flat文件(一条线,无头文件)

-zg= 重力值所在列

-zm= 磁数据所在列

In, out, -l 是必须输入的,-zg=-zm=必须输入一个

可选的参数可能用于指定列号,常数值,数据窗,多个水平层等

 

例如:XYZGMS data_in model_out -l=T1000 -x=1000,12000 -d=5 -eg=400 -zg=8

 

详细的可选参数:

-gzx,gzy,gzz,gxx,gyy,gxy=重力梯度值

-mx,my,mz=磁梯度值

-zeg=重力数据的高程值(向上为正)

-zem=磁数据高程值

-zt=地形值

-h=h1,h2,h3…最多15个水平层的列数据(向下为正),常数值

-eg= 重力观测值高程(向上为正)

-em= 

-t= 地形限制(默认通过整个范围)

-x=min,max X

-y=min,max Y

-s= min,max Station No.

其它

-u= 转换单位(默认为m(输入feet,-u=3.28)

-d= 抽取因子,每d个点取一个值(默认-d=1

-x0= 剖面上第一个点的x坐标,默认是0.

-q 版本信息

-f 输入一个flat文件,包括单一的线,不需要geosoft线的头文件

-j 这是一个时间模型,输出到*.tsr文件,不再是*.sur

 

XYZGMS 是一个命令行驱动的GM-SYS工具,用来读取GEOSOFT格式的带有多个组的 ‘XYZ’ 数据并创建GM-SYS模型。输入的数据组包括地形,,additional global horizons, 重力数据,磁数据和测点高程,是典型的Oasis montaj建模工具GX

值得注意的是,XYZGMS将多个线段转化为一个单一的剖面,黑线代表用户输入的原始数字化剖面,红线代表垂直投影到实际X,Y空间中模型剖面的位置。模型剖面方位角是基于原始剖面中的第一个和最后一个点。剖面位置随后进行拟合使拟合差最小

 Posted by at 上午 3:16
10月 232013
 

在windows上使用基于python的GeoIP数据集GeoLiteCity.dat,要使用这个,先把dat文件下载下来,在这里:http://geolite.maxmind.com/download/geoip/database/GeoLiteCity_CSV/GeoLiteCity-latest.zip,然后需要下载一个可用GeoIP的python库,官方提供了2种选择:http://www.maxmind.com/app/pythonhttps://github.com/appliedsec/pygeoip,我们选择后者pygeoip,因为前者对linux系统支持比较好,而windows还需要编译相关的库,我不想这样麻烦。后者则是完全基于python的,当把后者下载下来后,执行setup.py build和setup.py install,如果正常的话,就会完全安装好。如果没有提示正常安装,则可能是因为setup.py中引用的库不存在:setuptools,这个库是专门用来进行安装一些第三方库的,它的下载地址在:http://pypi.python.org/pypi/setuptools,选择对应的python版本即可。这时就可以使用一些pygeoip的例子代码了,更多的例子代码在这里:http://code.google.com/p/pygeoip/wiki/Usage

另外,可能还需要下载数据库,GeoIP.dat和GeoIPCity.dat可以在官方网站上找到,不过GeoIPISP.dat和GeoIPRegion.dat没有找到,我GG了一下,找到以下两个地址:http://linux16.net/upload/files/file/Any_files/2008/09/GeoIP.tar.gzhttp://torcache.net/torrent/CD948DD092BB435772A8EADA107107E8DE07C03A.torrent,如果这两者都不能下载时,试试用迅雷的离线下载功能吧,蛮好用的。好像还是没有找到GeoIPRegion.dat,不过它的功能好像GeoIPCity.dat都能提供。

import os
import pygeoip
gi = pygeoip.GeoIP(‘GeoIP.dat’)
gic = pygeoip.GeoIP(‘GeoLiteCity.dat’)
gio = pygeoip.GeoIP(‘GeoIPOrg.dat’)
gioo = pygeoip.GeoIP(‘GeoIPISP.dat’)
#gir = pygeoip.GeoIP(‘GeoIPRegion.dat’) file NOT found on web
print gi.country_code_by_name(‘google.com’)
print gi.country_name_by_name(‘google.com’)
print gic.record_by_name(‘google.com’)
print gio.org_by_name(‘baidu.com’)
print gioo.org_by_name(‘baidu.com’)

print
print gi.country_code_by_addr(‘64.233.161.99′)
print gi.country_name_by_addr(‘64.233.161.99′)
print gic.record_by_addr(‘64.233.161.99′)
print gio.org_by_addr(‘64.233.161.99′)
print gioo.org_by_addr(‘64.233.161.99′)

os.system(‘pause’)

 Posted by at 上午 1:41