普通物探重力勘探的数据处理PPT学习教案

1、会计学1普通物探重力勘探的数据处理普通物探重力勘探的数据处理第1页/共71页曲线分别围成的面积相等，重心不变；n 各点上的偏差应不超过实测异常的均方误差。第2页/共71页第3页/共71页多次线性内插圆滑第4页/共71页xaaxg10)(xiaaxig10)(第5页/共71页n从中解出系数 a0和 a1即可确定拟合式，从而计算拟合值。2201()()()()miniig i xg i xaa i xg i x0)()(2100ixigxiaaa0)()(2101xixigxiaaai第6页/共71页n这就是最小二乘法得到的三点圆滑公式。同理可得5点、7点等线性平滑公式。n拟合点数越多，得出的曲线

2、越平缓。)()0()(31)0(0 xggxgag)1 ()0() 1(31)0(gggg第7页/共71页1(0)17 (0)12 ( 1)(1)3 ( 2)(2)35gggggg2012( )g xaa xa x1(0)7 (0)6 ( 1)(1) 3 ( 2)(2) 2 ( 3)(3)21gggggggg第8页/共71页012( , )g x yaa xa y)0 , 1 () 1 , 0()0 , 1() 1, 0()0 , 0(51)0 , 0(gggggg1(0,0) ( 2,0)( 1,0)(0,0)(1,0)(2,0)9gggggg(0, 2)(0, 1)(0,1)(0,2)gg

3、gg1(0,0) ( 1, 1)( 1,0)( 1,1)(0, 1)9ggggg (0,0)(0,1)(1, 1)(1,0)(1,1)ggggg第9页/共71页22012345( , )g x yaa xa ya xa xya y52(0,0)(0,0) ( 1,0)(1,0)(0, 1)(0,1)99gggggg1 ( 1, 1)(1, 1)(1, 1)(1,1)9gggg xy-5 -4 -3 -2 -101 2 3 4 5-1-2-3123第10页/共71页第11页/共71页区域异常对局部异常的影响示意图 第12页/共71页第13页/共71页球体异常与单斜异常的叠加 第14页/共71页第

4、15页/共71页别用不同的趋势表示区域异常背景。图解法划分区域异常 第16页/共71页n第17页/共71页2)()()()(LxgLxgxgxgLLxxgg第18页/共71页niirgngrggg1)(1)0()()0()0()(rg)(rgi平均场法示意图 第19页/共71页(x1,y1)(x2,y1)(x3,y1)(x3,y2)(x2,y2)(x1,y2)( , )gg x y ( , )( , )( , )g x yg x yg x y 网格平均法示意图第20页/共71页第21页/共71页趋势分析法概略图 第22页/共71页重力低代表基底下陷的凹陷区。205700002058000020

5、59000020600000206100002062000020630000206400002065000020660000图 5 1 东 营 凹 陷 布 格 重 力 异 常 等 值 线 平 面 图410000041100004120000413000041400004150000416000041700004180000陈家庄滨县沾化县利津县垦利县东营市东营机场 辛镇郝家牛庄广饶县牛头镇博兴县高青平方王滨州市金家庄纯化镇坨庄现河庄王家岗永安镇草桥618.2591.710992.3第23页/共71页第24页/共71页东东 营营 凹凹 陷陷鲁西隆起区鲁西隆起区陈家庄凸起陈家庄凸起青坨子凸起青坨子

6、凸起广 饶 凸 起广 饶 凸 起20570000205800002059000020600000206100002062000020630000206400002065000020660000图 5 1 东 营 凹 陷 布 格 重 力 异 常 等 值 线 平 面 图410000041100004120000413000041400004150000416000041700004180000陈家庄滨县沾化县利津县垦利县东营市东营机场 辛镇郝家牛庄广饶县牛头镇博兴县高青平方王滨州市金家庄纯化镇坨庄现河庄王家岗永安镇草桥618.2591.710992.3第25页/共71页16500000165200

7、00165400001656000016580000166000001662000044000004420000444000044600004480000后坑湾窑多坝沟洞子湾南湖玉门关1650000016520000165400001656000016580000166000001662000044000004420000444000044600004480000后坑湾窑多坝沟洞子湾南湖玉门关1650000016520000165400001656000016580000166000001662000044000004420000444000044600004480000后坑湾窑多坝沟洞子湾南

8、湖玉门关1次趋势面2次趋势面4次趋势面7次趋势面第26页/共71页20570000205800002059000020600000206100002062000020630000206400002065000020660000图 5 3 东 营 凹 陷 剩 余 布 格 重 力 异 常 等 值 线 平 面 图410000041100004120000413000041400004150000416000041700004180000陈家庄滨县沾化县利津县垦利县东营市东营机场 辛镇郝家牛庄广饶县牛头镇博兴县高青平方王滨州市金家庄纯化镇坨庄现河庄王家岗永安镇草桥2057000020580000205

9、9000020600000206100002062000020630000206400002065000020660000图 5 2 东 营 凹 陷 区 域 布 格 重 力 异 常 等 值 线 平 面 图410000041100004120000413000041400004150000416000041700004180000陈家庄滨县沾化县利津县垦利县东营市东营机场 辛镇郝家牛庄广饶县牛头镇博兴县高青平方王滨州市金家庄纯化镇坨庄现河庄王家岗永安镇草桥东营凹陷区域布格重力异常东营凹陷剩余布格重力异常第27页/共71页Ma4Ma401Ma5Ma7Ma2Ma4021442514430144351

10、44401444514450144551446014465144704265427042754280428542904295第28页/共71页第29页/共71页上在边界在场源外空间SMfnuzuyuxuuS.).(|.022222222第30页/共71页)0 ,(21)0 ,(gf第31页/共71页 2/1222)()()0 ,(21),(zyxddgzyxV 2/3222)()()0 ,(2),(zyxddgzzVzyxg第32页/共71页 0202/322)()0 ,(2), 0 , 0(hrrdrdrghhg22)0 ,(), 0(hdghhg第33页/共71页|i| 算范围 n 的取值

11、应根据计算精度而定。nninnihihiitgihgdhhihghg344) 0 ,(1) 0 ,(1), 0(21)21()21(22)0 ,(ihg第34页/共71页201)0 ,(1)0 ,(21)(niirgndrgrg02/322)()(), 0 , 0(hrhrdrrghgnirriihrhrdrrghg02/3221)()(), 0 , 0(第35页/共71页irih)(rg1irir0)(), 0 , 0(iiirgKhg211210K1) 1(11) 1(12122iiKi,.3,.2,.1i第36页/共71页11()()2211222 3/2()()22( , ,0)( ,

12、 , )2()()jiijnnijijin jnzgd dg x y zxyz xy(i,j)第37页/共71页向上延拓使异常曲线变得平滑、平缓，相对突出了深部和区域地质因素的重力效应，压制了浅部和局部重力异常，相当于低通滤波的效果。第38页/共71页航磁异常重力异常第39页/共71页重力异常上延6070km重力异常上延100120km第40页/共71页向上延拓使异常曲线变得平滑、平缓，相对突出了深部和区域地质因素的重力效应，压制了浅部和局部重力异常，相当于低通滤波的效果。第41页/共71页滤波的效果。第42页/共71页2057000020580000205900002060000020610

13、0002062000020630000206400002065000020660000图 5 7 东 营 凹 陷 布 格 重 力 上 延 30KM 异 常 等 值 线 平 面 图410000041100004120000413000041400004150000416000041700004180000陈家庄滨县沾化县利津县垦利县东营市东营机场 辛镇郝家牛庄广饶县牛头镇博兴县高青平方王滨州市金家庄纯化镇坨庄现河庄王家岗永安镇草桥20570000205800002059000020600000206100002062000020630000206400002065000020660000图 5

14、5 东 营 凹 陷 布 格 重 力 上 延 5KM 异 常 等 值 线 平 面 图410000041100004120000413000041400004150000416000041700004180000陈家庄滨县沾化县利津县垦利县东营市东营机场 辛镇郝家牛庄广饶县牛头镇博兴县高青平方王滨州市金家庄纯化镇坨庄现河庄王家岗永安镇草桥上延5km上延30km向上延拓使异常曲线变得平滑、平缓，相对突出了深部和区域地质因素的向上延拓使异常曲线变得平滑、平缓，相对突出了深部和区域地质因素的重力效应，压制了浅部和局部重力异常，相当于低通滤波的效果。重力效应，压制了浅部和局部重力异常，相当于低通滤波的效果

15、。第43页/共71页)0 , 0 , 0(g(0,0,)gh(0,0,)gmhmiiihgzlzg0), 0 , 0()(), 0 , 0(第44页/共71页)()()(00jhihjhzzlmjmjiij 第45页/共71页反复迭代，直到延拓值与观测值的差达到要求精度；下延面上最终的重力异常值即为向下延拓结果 。第46页/共71页20570000205800002059000020600000206100002062000020630000206400002065000020660000图 5 9 东 营 凹 陷 布 格 重 力 下 延 1KM 异 常 等 值 线 平 面 图41000004

16、1100004120000413000041400004150000416000041700004180000陈家庄滨县沾化县利津县垦利县东营市东营机场 辛镇郝家牛庄广饶县牛头镇博兴县高青平方王滨州市金家庄纯化镇坨庄现河庄王家岗永安镇草桥20570000205800002059000020600000206100002062000020630000206400002065000020660000图 5 1 东 营 凹 陷 布 格 重 力 异 常 等 值 线 平 面 图410000041100004120000413000041400004150000416000041700004180000陈

17、家庄滨县沾化县利津县垦利县东营市东营机场 辛镇郝家牛庄广饶县牛头镇博兴县高青平方王滨州市金家庄纯化镇坨庄现河庄王家岗永安镇草桥618.2591.710992.3地面异常下延1km向下延拓相对突出浅部和局部异常，压制深部和区域异常，相当于高通滤波。第47页/共71页波。布格重力异常向下延拓平面图布格重力异常向下延拓平面图王万银 李跃秋王万银王万银收 集附图121.向下延拓深度为16km；2.等值线间距为5mgal。附图12正等值线及注记零等值线及注记负等值线及注记GPGT93-217GPGT93-217莺歌海盆地重、磁力资料处理及综合解释拟 编审 核项目负责图 号比 例 尺资料来源西安工程学院

18、中国海洋石油勘探开发研究中心1:4000000YIN6Y14-1-1112-HO-1X104-QN-1XY25-1-1YIN2LT13-1-1LT19-1-1LT9-1-1LT34-1-1YIN1YINQ2LT35-1-1IIK30-3-1AYX32-1-1LD30-1-1AY19-1-1Y26-1-1Y8-2-1Y7-4-1102-CQ-1X102-HD-1X103T-G-1X107-PA-1X104-QV-1X112AV-1X114-KT-1X117-STB-1X118-BT-1X120-CS-1X121-CM-1XY26-1-1DF29-1-2LG20-1-1LD20-1-1LD21-1

19、-1Y21-1-4DF1-1-6LD28-1-1ALD20-1-1Y35-1-1103T-H-1XLT1-1-1LD8-1-1LD8-1-3LD15-1-4DF1-1-9Y13-1-1112BT-1X118-CVX-1X119-CH-1X115-A-1XDF1-1-9DF1-1-1LD22-1-4油气井及编号油气显示井及编号干井及编号地震测线及编号图 例技术说明140015001600170018001900200021001400150016001700180019002000210010600107 00108 00109000011000106107 00108 0010900110 0

20、0静嘉边水奇英白龙尾岛东方莺歌海崖城三亚海口临高涠洲岛乌石莺歌海盆地重、磁力资料处理及综合解释1:400000001020km204060 km33布格重力异常平面图布格重力异常平面图王万银 李跃秋王万银王万银收 集附图021.海水密度取1.03g/cm ， 海底地形密度取2.20g/cm ；2.等值线间距为2mgal。附图02正等值线及注记零等值线及注记负等值线及注记GPGT93-217GPGT93-217莺歌海盆地重、磁力资料处理及综合解释拟 编审 核项目负责图 号比 例 尺资料来源西安工程学院 中国海洋石油勘探开发研究中心1:4000000YIN6Y14-1-1112-HO-1X104-

21、QN-1XY25-1-1YIN2LT13-1-1LT19-1-1LT9-1-1LT34-1-1YIN1YINQ2LT35-1-1IIK30-3-1AYX32-1-1LD30-1-1AY19-1-1Y26-1-1Y8-2-1Y7-4-1102-CQ-1X102-HD-1X103T-G-1X107-PA-1X104-QV-1X112AV-1X114-KT-1X117-STB-1X118-BT-1X120-CS-1X121-CM-1XY26-1-1DF29-1-2LG20-1-1LD20-1-1LD21-1-1Y21-1-4DF1-1-6LD28-1-1ALD20-1-1Y35-1-1103T-H-

22、1XLT1-1-1LD8-1-1LD8-1-3LD15-1-4DF1-1-9Y13-1-1112BT-1X118-CVX-1X119-CH-1X115-A-1XDF1-1-9DF1-1-1LD22-1-4油气井及编号油气显示井及编号干井及编号地震测线及编号图 例技术说明140015001600170018001900200021001400150016001700180019002000210010600107 00108 00109000011000106107 00108 0010900110 00静嘉边水奇英白龙尾岛东方莺歌海崖城三亚海口临高涠洲岛乌石莺歌海盆地重、磁力资料处理及综合解释

23、1:400000001020km204060 km第48页/共71页当下延深度较大时不宜一次完成，应采取逐次圆滑，逐次延拓的方法。第49页/共71页n另外延拓理论也还有待于进一步发展、完善。第50页/共71页32222( , ,0)( , , )2()()zgd dFg x y zFzyz 122221( , ,0)2()()gd dFzzyz 2 2 (,0)1(, )(,0)2zzxykzkzxyxyzg k kg k kzeg k kezk 第51页/共71页22( ,0)( , )2()zgFg x zFdxz2 |(, )(,0)xk zxxg kzg ke 第52页/共71页xzV

24、xgzzVzgzzzVzzg用导数区分叠加异常 第53页/共71页xzxVg ()()(0)2xgxgxgx两个不同大小、深度球体的异常及其导数第54页/共71页0( )mkkkg xaxmkkkxxzxakxgxV11)()()2()2(51)()(101) 0(xgxgxxgxgxVxz第55页/共71页第56页/共71页第57页/共71页第58页/共71页zzzVg 2/1222)()()0 ,(21),(zyxddgzyxV 2/3222)()()0 ,(2),(),(zyxddgzzyxgzyxVz 2/5222222)()2)(0 ,(21), 0 , 0(), 0 , 0(zdd

25、zgzgzVzzz第59页/共71页n实际计算也是采用各环平均值加权求和的方式。2002/52222)()2()0 ,(21), 0 , 0(), 0 , 0(zrrdrzrrgdzgzVzzz20)0 ,(21)0 ,(drgrg02/52222)()2()0 ,(), 0 , 0(), 0 , 0(hrrdrhrrghghVzzz021)0 ,()0 , 0 , 0()0 , 0 , 0(drrrggVzzz第60页/共71页所以有n这就是由重力观测值换算的理论公式。02222222zVyVxVV22222220ggggxyz)(222222ygxgzggzzzzg第61页/共71页)()0(4)0(2RggRgzz)5(48)2(12)(16) 0(44621) 0(2R