元胞自动机在农村土地利用动态模拟中的应用研究-

1、2007年第10期福建电脑元胞自动机在农村土地利用动态模拟中的应用研究成筠1,2,陶刘强2,刘文宝1【摘要】:本文从元胞自动机的原理、特征出发探讨了元胞自动机(Cellular Automata在模拟农村土地利用演化格局应用中的可能性,并对利用元胞自动机进行农村土地动态模拟的核心问题-转换规则进行了初步探讨。【关键字】:元胞自动机;转换规则;农村土地利用;动态模拟1.引言随着全球变化研究的深入,人们认识到土地利用/覆盖变化(Land Use Changes/Covers,简称LUCC是造成全球气候变化的重要原因。中国作为一个农业大国,对土地利用/土地覆盖的研究已经有很长的时间,但是以前的很多研

2、究多侧重土地资源的调查、分区、分类、规划等等。农村土地动态演化具有明显的复杂性和非线形特征,属于复杂巨系统。显然,靠牛顿力学等确定型数理模型的"移植"来建立的地理模型难以模拟这种复杂现象。系统动态学模型在地球复杂系统研究中应用较为广泛。但系统动态学模型也有"先天不足",表现在缺乏对空间因素的处理,同时,它是基于微分方程的,不利于GIS环境下以离散为特征的计算机建模。近年来,建立于微观动态模拟理论基础上的一些微观离散模拟方法如CA(Cellular Automata等被引入到GIS环境建立时空动态模拟分析引擎。目前的CA都用在模拟城市交通流、城市扩展和林火

3、的模拟和预测。基于此,本文从农村土地利用/覆盖的特点出发,提出用CA模拟农村土地的动态演化机制。2.CA及其在模拟农村土地利用变化中的可行性2.1CA及其特点元胞自动机是一种具有时空计算特征的动力学模型,或者说是一种时间、空间、状态都离散,空间上相互作用及其时间上因果关系皆局部的网络动力学模型。标准的CA系统是一个由元胞(Cell、状态(States、邻域(Neighbors、和规则(Rules构成的四元组。所有离散的元胞构成一个元胞空间;在某一时刻一个元胞只能取有限的离散状态集中的一种;邻域是元胞按照一定的形状划分的元胞集合,它影响该元胞下一时刻的状态;转换规则定义了元胞类型间相互转化的函数

4、。CA的时间也是离散时间,不具有物理意义。用集合语言将CA模型描述为:St+1=f(S t,N,式中,S为有限集合,代表元胞状态;N代表元胞邻域;t表示时间;f为局部转换规则。CA模型以其框架简单、开放和可以模拟十分复杂系统的行为而具有很强的生命力,从目前的研究看,它具有以下特点:(1空间性,反映在以栅格单元空间来定义元胞自动机,因而能很好地和许多空间数据相互兼容。(2离散性,即时间、空间、状态都是离散的。(3同步性,可以将元胞自动机的状态变化看成是对数据或信息的计算或处理,这时的重要特征是计算并行性。(4局部性,每个元胞的状态,只对于其周围半径为r的邻域内的元胞下一时刻的状态有影响。(5高维

5、数,在动力系统中一般将变量的个数称为维数,从这个角度看,元胞自动机的维数是无穷维的。2.2CA在模拟农村土地动态演变中应用的可行性本文主要从以下几个方面阐明CA在模拟农村土地动态演变应用中的可行性:(1CA采用"自下而上"的建模方式,而且没有既定的数学方程,只是一个建模原则,因而具有很好的开放性和灵活性。符合人们认识复杂事物的思维方式。(2CA是一个基于微观个体相互作用的时空动态模型,将地理实体的时空特性统一到模型中,通过研究对象的元胞空间划分、初始状态确定及转换规则,CA就可以自行迭代模拟土地演化。(3CA将时间和空间离散化,适合于建立计算机模型和并行计算特征。(4CA具

6、有不依赖比例尺的概念,因此可以模拟局部的,区域的或大陆级的演化过程。(5CA模型中的元胞和基于栅格GIS中的栅格一样,所以容易将二者集成。3.CA模型在农村土地利用动态模拟中的应用农村土地利用的动态演化属于复杂的地理现象,为了使CA 能够模拟和分析农村土地利用演化,须将标准CA的诸要素进行"地理化"形成地理元胞自动机(GeoCA参见图1:图1地理元胞自动机模型:扩展的元胞自动机模型53.1农村土地动态演变CA模型的扩展(1地理元胞及状态的扩展,标准CA中元胞和元胞的状态是抽象概念,在本模型中必须给元胞及其状态赋予特定的地理含义:元胞为农村土地划分的地块,元胞的状态是土地利用

7、类型,如:城镇建设用地、耕地、林地。(2元胞空间的扩展。元胞空间不再是抽象的空间,而是与笛卡儿坐标系下的实际地理空间相对应,并抽象为二维地域,同时把元胞空间划分为规则的矩形,考虑到元胞空间和笛卡儿坐标系下的地理实体对应,不同的元胞大小表现出不同的地理宏观规律,因此1Km×1Km大小的地理单元为元胞。(3邻域定义。本模型中采用Moore邻域作为元胞的邻域,这个标准CA中采用的邻域定义是一样的。其邻域定义如下: M M oore=v i=(v ix,v iy|v ix-v ox|1,|v iy-v oy|1,(v ix,v iyZ2(4转换规则的扩展。转换规则是CA模型的核心,它决定了C

8、A的动态演化过程。(5时间概念的扩展。CA模型中的模拟时间必须和地理实体演化中的真实时间建立对应关系,否则,地理时空建模就失去价值。一般采用历史数据和其它模型预测数据建立而这之间的联系。3.2农村土地动态演变CA模型的转换规则的确定地理元胞自动机(GeoCA中的转换规则是地理特征和规律在局部和微观上的体现。当然,大部分地理规律是宏观规律,一方面很难确定其在局部和微观上的具体形式,另一方面我们期932007年第10期福建电脑望通过局部间的相互作用而得到某种"突变"。所以,根据模型的规则不能推理到宏观地理规律的。因此,它与宏观的地理规律既有联系又有区别。鉴于农村土地动态变化十分

9、复杂,除受局部个体间相互作用的影响,还受各种区域的和更大尺度因素的影响特点,CA模型中状态转换规则必须兼顾微观、区观和宏观,建立综合的多层次规划,才能有效地模拟农村土地动态演化过程。而且为了更真实地模拟其变化过程,本模型采用随机CA模型,即在元胞状态和邻域构成确定的条件下,元胞的状态转换不是确定的,而是受一定的概率控制。元胞转换规则的扩展是GeoCA和标准CA最根本的区别。首先,从微观层次,主要考虑邻域函数这个指标。邻域函数是用于刻画当前元胞受周围邻域元胞影响程度的指标,可以用下面的公式表示(1(1;这里P(x,y 邻域构形内某种土地类型元胞所占百分比,N表示邻域构形内的元胞个数,X(x,y表

10、示某中土地类型的元胞,表示邻域空间。其次,从中观层次考察耕地、城镇建设用地和其它类型的土地使用之间转换的驱动因子。1从耕地的角度来看:从耕地减少面积统计分析,我国耕地面积减少的主要原因如下:(1城镇建设用地。城镇是带动农村经济的重要推动力,自1992年以来,城镇建设占用耕地的比例有较大幅度的上升而且在较长时期内还有上升的势头。(2乡镇企业的发。乡镇企业发展一直是我国耕地减少的主要原因,尤其是工业型的乡镇企业的发展。(3个人建房。在农村地区,往往对土地资源的重要性认识不足,视土地资源为公有财产,随意乱占、滥用土地,加速了耕地资源的流失。(4农村人口增加。人口密度与土地利用覆盖变化速率呈正相关关系

11、,人口增长速度越快,土地利用变化也越快。由于人口数量大增,在生产条件没有大幅度提高的条件下,导致耕地被大量占用。(5农业结构内部调整。农业内部结构调整在我国耕地减少中占较大比例,特别,自1986年土地管理法的实施,同时随着人们对环境意识的提高,耕地减少的现象开始减少。但自1992年后在失常经济的作用下,加之粮价低下很多优质粮田改为鱼塘,果园以及经济开发区,加速了耕地的流失,影响农业生产。(6自然灾害。长期以来,土地退化(主要包括水土流失、海水入侵、土地沙化、土壤盐碱化、土壤污染等是我国面临的主要问题,每年因此损失的大量耕地资源。2从城镇建设用地角度看。(1人口增加。随着城镇人口的增加,城镇的发

12、展将向四周扩展,且在相当长的时间内城镇用地不被用做其它用途。(2经济繁荣。发达的经济是城镇扩展的必备条件,一般来说经济越发达其城镇规模也越大。(3交通的通达性。交通也是城镇外扩的重要途径,城镇扩展形式与交通呈正相关,城镇将沿着交通线发展。3从林地类型角度看,其发生和发展除受到区域自然环境的直接控制,主要包括气候、地貌、土壤等以及原始林地覆盖类型外,还受人类活动的干扰,主要表现为人类影响的可达性与频率及对林地的利用方式的选择,包括人口分布,交通状况、社会需求、经济结构的调整等。综合对中观层次上的驱动力因子利用专家打分法给出权值Wk,其中k代表各个驱动因子。最后利用公式(2计算土地类型的转换系数。

13、(2其中C t ij是在位置(i,j处t时刻土地类型转换系数;随后采用下面的负幂函数将结果映射为0,1范围内的土地类型转换概率P。P=aexp-?C(3其中a表示标准化常数,?表示一种调节系数。最后,在宏观层次上,考虑时态尺度上的影响因子,如由土地可持续利用导出的土地资源代际分配,控制城市发展进程的规划和国家有意识地为保护耕地、林地而划分的保护区。考虑以上因素,整个模型可以用集合的语言描述为:St+1=f (,St,(N,C,R n,R p其中为元胞空间,可以是邻域、层次区域、和全局区域;St+1和St分别为元胞t+1和t时刻的状态;C 为转换系数;Rn,Rp为自然条件和政策规划的约束;f为规则映射函数。本模型综合考虑了微观、区观和宏观的作用因素,建立了更加全面的转换规则,避