
含油气盆地分析课件盆地埋藏史.docx
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底界(上剥蚀) 。












