地质过程中的定量方法与计算技术-地温场与热史恢复共98张课件.ppt
- 【下载声明】
1. 本站全部试题类文档,若标题没写含答案,则无答案;标题注明含答案的文档,主观题也可能无答案。请谨慎下单,一旦售出,不予退换。
2. 本站全部PPT文档均不含视频和音频,PPT中出现的音频或视频标识(或文字)仅表示流程,实际无音频或视频文件。请谨慎下单,一旦售出,不予退换。
3. 本页资料《地质过程中的定量方法与计算技术-地温场与热史恢复共98张课件.ppt》由用户(晟晟文业)主动上传,其收益全归该用户。163文库仅提供信息存储空间,仅对该用户上传内容的表现方式做保护处理,对上传内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知163文库(点击联系客服),我们立即给予删除!
4. 请根据预览情况,自愿下载本文。本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
5. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007及以上版本和PDF阅读器,压缩文件请下载最新的WinRAR软件解压。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 地质 过程 中的 定量 方法 计算 技术 地温 恢复 98 课件
- 资源描述:
-
1、地温场与热史恢复石油有机成油论干酪根有机成油论有机质未成熟60-180摄氏度 油裂解为气小于60摄氏度 大于180摄氏度生成石油地温场与热史恢复地温场的一般知识不同盆地类型地温场及演化特征热史重建构造热演化法古温标法结合法地温场的一般知识地温场的一般知识1.热源热源A 地幔热源(地幔热源(60%),变化性,变化性B 放射性元素生热放射性元素生热(40%),稳定性稳定性(10Km厚结晶岩表层厚结晶岩表层)C 其他热源,随机性(其他热源,随机性(火山作用,构造作用,化学热火山作用,构造作用,化学热)2.热源的传输持续时间热源的传输持续时间据据 Lachenbruch 的估算:对于的估算:对于100
2、 km的岩石圈,一般热的岩石圈,一般热传输持续时间为传输持续时间为50-100 Ma,因此,岩石圈具有相对较高因此,岩石圈具有相对较高的热惯性,背景热状态持续的时间可与含油气系统的寿命的热惯性,背景热状态持续的时间可与含油气系统的寿命期相当,所以,确定现今热状态作为起点,反推解释油气期相当,所以,确定现今热状态作为起点,反推解释油气形成时的可能热状态通常是合理的。形成时的可能热状态通常是合理的。t=r2/4a,a=32大陆大陆地壳软流圈大洋大洋岩石圈中间层相界低速带莫霍面100500厚度(km)地地球球的的主主要要分分层层及及其其流流变变学学分分界界(Allen,1990)地温场的一般知识3.
3、地温场的形成机制A 热传导型地温场(控制区域地温场)B 热对流型地温场(增温型,冷却型),造成局部异常C 热辐射对地温场的影响:只影响地表温度依靠物体中的微观粒子的热运动传递能量的过程叫热传递由于流体从空间某一区域移动到另一温度不同的区域时,发生的能量转移的过程由电磁波来传递能量的方式叫热辐射地温场地温场稳态地温场非稳态地温场如果地温场中各点的如果地温场中各点的温度不随时间而变化,温度不随时间而变化,则称其为稳态地温场则称其为稳态地温场如果地温场中各点的如果地温场中各点的温度随时间而变化,温度随时间而变化,则称其为非稳态地温场则称其为非稳态地温场仅用简单的数理模型即仅用简单的数理模型即可进行定
4、量描述和计算可进行定量描述和计算数理模型复杂,不一定有数理模型复杂,不一定有解析解,只有数值解解析解,只有数值解地幔增温带(内热带)恒温带(中性带)变温带(外热带)地表主要受太阳辐射的影响而发生的主要受太阳辐射的影响而发生的日变化、月变化、年变化等日变化、月变化、年变化等地球内热和太阳辐射热的平衡地带,地球内热和太阳辐射热的平衡地带,一般很薄,可视为一个面一般很薄,可视为一个面主要受地球内热的影响。一般主要受地球内热的影响。一般而言,越向深处,温度越高而言,越向深处,温度越高地表浅层的温度场结构地表浅层的温度场结构地温梯度:G:为地温梯度;为地温梯度;T:为深度为为深度为H出的温度出的温度;T
5、o:为地表平均温度(或恒温层);为地表平均温度(或恒温层);H:为深度为深度。100HTTGo大地热流:dHdTKqK:为热导率;dT/dH:为地温梯度;单位为毫瓦/平方米(mw/m2),这是目前国际通行的单位,它与以前的热流单位HFU(Heat Flow Unite)的关系是1HFU=41.86(mw/m2),1HFU=1微卡/平方厘米 秒;式中式中“-”号表号表示热流方向与地示热流方向与地温梯度方向相反温梯度方向相反HoHTTHKdHqTTHKdHqdTKdHqdTo00)()(或采用或采用离散形式离散形式的值。及按最小二乘法可求得若有一组观察值qTniKHTKHqTToiiiiio),.
6、,)(,(321qnkHnTTTnTTkHnTkHqiiioiiiiiiii2211)(布拉德法布拉德法(适用于稳定地温场)(适用于稳定地温场)由热流公式,有参数的最小二乘法估计niinktjiniyyxyQXYXYQj1212112)()(亦即考虑既有关系:的最小二乘法估计量。为这时参数的点估计。作为达到最小值的寻求使得iikkQ .,.,2121值。可由一正规方程组求出由jjnktjinktjiQxyxyjkj ,min).,(,012112121最最小小二二乘乘法法原原理理热传导方程热传导方程.1热传导方程的导出热传导方程的导出 考察空间某物体考察空间某物体G的热传导问题,以函数的热传导
7、问题,以函数u(x,y,z,t)表示物体表示物体G在位置在位置(x,y,z)处及时刻处及时刻t的温度。的温度。据据Fourier定律,物体在无穷小时段定律,物体在无穷小时段dt内流过无穷小内流过无穷小面积面积ds的热量的热量dQ与物体温度沿曲面法线方向的方向导数与物体温度沿曲面法线方向的方向导数 成正比,即成正比,即 nudsdtnutzyxkdQ),(在G内任取一闭曲面内任取一闭曲面 ,所围成的面积记为,所围成的面积记为 ,则从,则从t1t2流进流进此曲面的全部热量为此曲面的全部热量为dtdsnuzyxkQtt 21),(曲面的外法向导数另外,另外,由于流入的热量使物体内部发生变化,在时间间
8、隔由于流入的热量使物体内部发生变化,在时间间隔t=t1-t2内(假定内(假定t2t1),物体内点(物体内点(x,y,z)处的温度,由处的温度,由吸收热量为内则有在为比热为密度;假定:温度差为升为 121221tttctzyxutzyxutzyxutzyxu),(),(),(),(dxdydztzyxutzyxucdsdtnukdxdydztzyxutzyxuzyxzyxcQtt),(),(),(),(),(),(121221 从而有dxdydzdttucdtdxdydzzukzyukyxukxdsdtnuktttttt)()()()(212121 假设假设 u关于变量x,y,z具有二阶连续偏导
9、数,关于t具有一阶连续偏导数,利用奥-高公式,有改变积分次序,可得到021dtdxdydzzukzyukyxukxnuctt)()()()()()(zukzyukyxukxnuc 的任意性,可有和由于21tt,上式称为非均匀的各项同性体的热传导方程。上式称为非均匀的各项同性体的热传导方程。),22222222222222zuyuxuazuyuxucknuackck(,有为常数,记及如果物体是均匀的,即 如果内部有附加热源,设单位时间内单位体积中所产生的热量为F(x,y,z,t),则有 ctzyxFtzyxftzyxfzuyuxuatzyxFzuyuxuanu),(),(),(),()其中,(2
10、2222222222222齐次非齐次前述原理与方法仅实用于较简单的一维稳定地温场。对于非稳定地温场,则需要用更复杂的方程来描述。以三维非稳定地温场为例,可用斯托曼方程(Stallman,1963)描述:tTcTvzTvyTvxcqzTyTxTkrrzyxwwo )()()()(222222热导率水的密度和比热容附加热源或汇岩石的密度和比热容上式包含上式包含3 3方面:传导方面:传导+附加附加+对流对流tTcqzTyTxTkrro)(222222若不存在热对流,有若不存在热对流,有kqzTyTxTo222222若为稳定地温场,则为Poisson方程(椭圆方程)若附加热源为零,则有0222222z
11、TyTxTLaplace方程(调和方程)0222222zTyTxTT即用哈密顿算符,有02 T岩石的热导率岩石的热导率定义:单位时间内流过单位面积的热量与温度梯度负值之比。定义:单位时间内流过单位面积的热量与温度梯度负值之比。)(dZdTAdtdQK热导率的因次为:(能量)(时间)-1(距离)-1(温度)-1。在地学研究中通常使用两种单位:法定计量单位和TCU单位。两种单位因次和换算关系如下:法定单位(w/mk);常用计量单位(TCU)(cal/cm cal/cm s s o oc c)or(mcal/cm/cm s s o oc c)1TCU=0.4186w/m k,1TCU=1mcal/c
12、m cm s s o oc c 1 w/mk=2.388TCU AALCCLllAALCt=0lt/8ltC-CrustL-LithosphereA-AsthenosphereT(0C)T(0C)T(0C)1333 0C1333 0C1333 0C莫霍面以上莫霍面以上古滕堡古滕堡面以上面以上l一维热流方程一维热流方程22ZTDtTT:古地温;Z:以岩石圈底界为原点,直至地表的垂直坐标;t:以拉张发生时间为零,至今的时间;D:岩石圈的热扩散系数0 lZtZT),(mZTtT 00),(边边界界条条件件初初始始条条件件mtTtZT 0),()(),(lZTtZTmt10 111lZ)()(110l
13、Z解热流解热流(热传导热传导)方程,可得岩石圈的温度分布方程,可得岩石圈的温度分布)exp(sin)sin()()(),(tnlZnnnnTlZTtZTnnmm211121)热流:表(立叶定理很快可求得地为热松弛时间。根据傅其中lZDl22 )exp(sin)(1221nmtnnnlkTtq 思考题:怎样由T(Z,t)T(Z,t)求得求得q(t)q(t)?当扩张开始当扩张开始t=0,当当t趋于无穷时,地表由于岩石圈冷却收缩趋于无穷时,地表由于岩石圈冷却收缩沉降到最终的深度,其间,地表相对于这个深度的高度为:沉降到最终的深度,其间,地表相对于这个深度的高度为:)(exp)(sin)()(0220
展开阅读全文