文档介绍:该【自然与环境数学模型Read 】是由【玉面小达摩】上传分享,文档一共【4】页,该文档可以免费在线阅读,需要了解更多关于【自然与环境数学模型Read 】的内容,可以使用淘豆网的站内搜索功能,选择自己适合的文档,以下文字是截取该文章内的部分文字,如需要获得完整电子版,请下载此文档到您的设备,方便您编辑和打印。第五章 自然与环境的数学模型§ 放射性计年模型一. :原子中核电荷数(质子数)相同,但具有不同的质量(中子数)的元素称为同位素。它们的化学性质相同,在周期表中处于同一个位置。有相同的元素符号,但在左上角注名质量数,左下角为质子数。如:123H,12C,13C,14C,234U,238U等。1H,H,。不稳定的同位素具有放射性,通过放射粒子而变为同一元素的不同的同位素或者不同元素的同位素。称之为同位素的蜕变。14148787238234230Th6C7N,37Rb38Sr,UU03. 放射性同位素衰变的数学模型:粒子同质、以相同的概率衰变,粒子群体规模充分大,没有粒子迁移,粒子的衰变率反映群体衰变现象的平均效应,单位时间内衰变正比于群体的数量,比率系数为常数。参量,变量:为衰变系数,N(t)为t时刻粒子数。模型:dN/dt=-N(t),N(0)=N0有解:N(t)=N0e-t称使得N(T1/2)=N0/2的时间T1/2为半衰期:T1/2=ln2/同位素半衰期(年)衰变系数(/年)××10-××10-×10-×10--46C×,又得到定常速率的补充,则放射性元素的变化可以由如下模型来描述dN/dt=-N(t)+(t)=b/+[N(0)-b/]e-t,模型有一个平衡状态N*=b/。当N(0)>N*时,元素不断减少,当N(0)<N*时,,N(t)N*表明放射性元素最终将会达到这个平衡状态,我们称这个现象为放射性元素的平衡。:湖南省长沙市马王堆一号墓于1972年8月出土。考古学家利用出土的木炭确定墓葬的年代。考古中的14C计年法:利用死亡的生物体内14C的含量测定死亡的时间。14C的半衰期(5730)相对于地质年龄很短,但适用于考古学。原理:14N蜕变为14C。,产生中子, 14CO2被植物吸收,进入食物链。生物体内的14C不断衰变为14N成为气体散失。活体还会不断吸收14C,使体内的14C维持平衡。生物死亡后,不再吸收14C,则体内的14C将由于蜕变,在体内的比例数将逐渐降低。分析:记N(t)为t时刻生物体内14C的粒子数。如果生物体在t=0时死亡,则有N(0)=N*(活体内14C的平衡量)。如果在t=t1时物体被发掘,由N(t)=N*e-λt可以得到t1=ln[N*/N(t1)]/λ因为14C与14N具有相同的原子质量,无法使用质谱仪测量样本中14C的含量,只能用计数管测量14C的放射蜕变物的速率。假设:1. 马王堆墓葬的年代生物体中 14C的平衡数量与现代生物体中 14C的平衡数量相同。。,。于是,t1=ln[N*/N(t)]/λ=ln[N*’/N’(t)]T1/2/ln211=ln()×5730/=:1650乌斯赫基于神学的计算世界在-4004年诞生。1850汤姆逊基于地球冷却学说,地球的年龄在2000万年到4000万年之间。1896被贝克勒尔否定。1903居里利用放射性测定地质年代。1907伯尔伍德得到了放射性地质年代的数据。通过岩石的地质年表及同位素的演化规律确定更可靠的地球年龄。可用于地质计年的同位素:,,,。难点:无法确定地球形成初期岩石中元素的含量N*,因此无法运用公式t=ln[N*/N(t)]/λ确定岩石的形成期。:。 , 不再衰变。。。即N(t)+D(t)=N0+,N(t)表示在t时刻放射性元素含量,D(t)表示在t时刻衰变后稳定产物含量,N0,D0表示在t=0时的初始值,由放射性元素的衰变规律,N(t)=N0e-t,得到D(t)=D0+N(t)et-N(t)=D0+(et-1)N(t)岩石放射性计年的等时线模型记m=et-1,r0=D0则对于固定的t和不同的D0,N(t)和D(t)将分布在同一条直线—等时线上。(IsochronDiagram)D(t)=mN(t)+r0,如果同时在岩石所含的n种不同的矿石中观测到不同元素的含量(Ni,Di),i=1,2,,n,则利用最小二乘法就可以给出m,r0的最小二乘估计。从而得到岩石生成的时间t=1/ln(m+1)=T1/2ln(m+1)/(m+1)T1/2难点:无法测得指定元素的绝对含量,用质谱仪只能测的两种元素的相对含量。因此,改用N(t)表示在t时刻放射性元素与稳定元素的相对含量,用D(t)表示在t时刻衰变后的稳定产物与另一种稳定元素的相对含量。岩石年龄的铷87Rb— (87Rb)只因放射性而衰减为稳定的锶 (87Sr),(87Rb)的半衰期很长,为 ×109年, (87Rb)和锶(87Sr)的含量不同,(86Sr)的含量是稳定的, (87Rb)和锶(87Sr)的相对含量 Rb=87Rb/86Sr和Sr=87Sr/86Sr。、白云母和黑云母的岩石。用质谱仪进行分析,测的铷和锶的相对含量如下:=,r0=,从而有t=××.610×9年=×108年。地球的年龄约为44~46亿年(×109),宇宙约为150~200亿年(2×1010)。问题:P143第3题:铀–钍计年法§:问题:建模描述湖泊污染的状况。背景:湖泊:提供水源,水产养殖,交通运输,,,:水体覆盖面积大,污染源复杂,,不利于水体的更新和自净..假设:,以污染物的含量标志污染的状况 .,单流出, . . .三. 建模物理理想模型: 池水含盐问题数学模型:湖水体积:V,污染物浓度: P(t),流入速度:r,流入污染物浓度:PI(t),流出速度:r,流出污染物浓度:P(t).则有V/rdP/dt=PI(t)-P(t)其中=V/:得P(t)=(Ps-K)e-t/情形I:自由倾倒PI=K,P(0)=Ps.+(t)K(t),称K为饱和污染状况。当Ps<K时,P(t)增加,当Ps>K时,P(t)减少。称(t)=P(t)/K为湖水在时刻t的污染水平。不难得到(t)=(Ps-K)/Ke-t/+1当=1时,称为饱和水平;当>1时,称为超饱和状态,P(t)将会下降。=0(一池清水),则t时刻的污染水平为(t)=P(t)/K=1-e-t/给定<1,记T为达到水平污染的时间,则有T=ln(1-),当=?时,有T1/2=ln(2),有T1/2=21年。对于苏比利尔湖,有T1/2=132年。一般来说,对于PsK,若给定s:=Ps/k<<1则有T=ln[(1-s)/(1-)]=0且Ps>0,则P(t)=Pse-t/(t)=P(t)/Ps,=ln[Ps/P(t)]=ln(1/a)=1/2,则有T1/2=,在完全断绝污染物流入的前提下,:PI(t)=K0e-,污染状况以强度a逐年得到控制.-(0)=K0,模型:dP/dt=[-P(t)+K0e]/则模型有解P(t)=K0(e-at-ae-t/)/(1-a)=K0e-at(1-ae(a-1)t/)/(1-a)由此不难证明,dP/dt<0,而且有P(t)0(t)。它表明只要控制污染的力度足够大湖水的污染程度将会不断得到改善。情形III混合情形:在初期,湖泊属于自由污染阶段,当湖水被污染到一定的水平,将对污染源加强管理和控制,降低排污量。模型:如果在初期我们有Ps=0,PI=K1,则湖水将在t=-ln(1-)(<1)达到水平的污染,即有P(t)=K1。此后对污染源加强管理,将排污量降低为PI(t)=K2<K1tt则模型有解K1(1et/),ttP(t)(tt)/K2,tt(K1K2)e五. 讨论1. 蒸发与渗漏:输出正比于湖水表面积 , 因此正比与湖水体积的 2/3次方 。离散动态:变连续的微分方程为离散的差分方程。污染物的变化:例如:DDT被动物吸收,溶解于脂肪中;有机磷引起水藻激增,贮存于水藻体内,当水藻腐烂时,有机磷又游离于水中。 4. 扩散过程:污染物浓度是空间点的函数p(t,x)污染物流速 J(t,x)=-D污染物扩散方程 p(t,x)/ t=-当D是常数时, p(t,x)/ t=Dxp(t,x) ,其中D表示污染物扩散系数x·J(t,x) , 因此 p(t,x)/ t= x·(DDxp(t,x) 。xp(t,x))。问题:P143第4题:伊利湖和安大略湖的污染。