含油气盆地分析课件 盆地埋藏史
- 格式:doc
- 大小:220.50 KB
- 文档页数:5
负指数模型中参数使用1stopt 解方程组得到
其中a=0.48,c=0.000410714560430877,简化取c=0.0004
公式推导得知
⎰
2
h 1
h (1-ψ(h ))dh=
⎰
4
h 3
h (1-ψ(h ))dh
假设沉积水深一直为230m 。
故在此时h1,h2,h3已知仅h4未知,因手中1stopt 版本较低,无法解决带积分的隐函数求解,故使用matlbe ,程序如下。
例如:当恢复层一的地表真实厚度时,h1=5421m ,h2=6643m ,h3=250m 首先定义一个M-file 文件如下:
function f= hanshu(h1,h2,h3,h4) syms y1 y2 y3 h global y4 z
y1=1-0.48*exp(-0.0004*h)/*负指数压实模型*/ y2=int(y1,h,h1,h2)/*岩层1现今埋深骨架体积*/ y3=int(y1,h,h3,h4) /*岩层1原始地表条件骨架体积*/ y4=y2-y3
z=solve(y4,x4)/*骨架体积不变*/ eval(z)
主程序如下: Syms x4;
Hanshu (5421,6643,250,x4);
故得到岩层1层底界(盆地基底面)的沉降历史:
线性模型(未用)
负指数模型(使用)
)
(h *0.0004-exp *48.0y =
时间(Ma)基底深度
40.3 -230
37.6 -1945.8538
33.1 -3062.4
28.4 -4656.1
24.6 -5322.5
17.5 -6231.9
15.4 -6016.9
10.3 -6504.752
0 -6643
由于均恒作用产生的基底沉降表征公式为:
时间(Ma)该时间段内沉降深度(m)累计沉降深度(m)
40.300
37.6872.3996288872.3996
33.1625.50180081497.901
28.4961.89416682459.796
24.6409.79298612869.589
17.5569.54905743439.138
15.4-134.1528143304.985
10.3306.84290993611.828
087.126057523698.954
时间(Ma)
基底构造沉降累计量该时间段内基底构造沉降
(m)
40.300 37.61073.4541711073.454171 33.11564.49857491.0443992 28.42196.304404631.8058332 24.62452.911417256.6070139 17.52792.76236339.8509426 15.42711.915174-80.84718595 10.32892.924264181.0090901
02944.04620751.12194248
基底沉降历史
基底构造沉降史:
由基底构造沉降历史图件分析可知:早期为较快速沉降,晚期沉降速度较慢,二者之间存在一个明显的不整合,故判断为裂谷型盆地,早期为裂陷阶段,快速沉降,其沉降受构造(以断层为主)控制,晚期为热沉降,沉降速度较慢,二者之间存在明显的不整合!此类盆地具有较高的地温梯度,类似于我国东部渤海湾及松辽等。
埋藏史
层4底界
层6底界(上剥蚀)。