山东羊庄岩溶水系统饮用水水源地保护区划分探讨

垫塑　Ｑ：垫　Ｓｃｉｅｎｃｅ　ａｎｄ　Ｔｅｃｈｎｏｌｏｇｙ　Ｉｎｎｏｖａｔｉｏｎ　Ｈｅｒａｌｄ　资源与环境　山东羊庄岩溶水系统饮用水水源地保护区划分探讨　卜华　陈占成‘　张良鹏　（１．山东省鲁南地质工程勘察院　兖州　２７２１　００；２．山东省临沂市水利局　临沂　２７６００１）　摘要：羊庄岩溶水系统是我国北方岩溶发育区较为典型的系统之一，多年来为地方工农业发展提供了丰富的岩溶水资源；为解决枣庄市　和滕州市城市供水矛盾，又规划向两市区供给饮用水源　而近几年的水质监洲资料显示，岩溶水已经受到了一定程度的污栗，多项化学指　标含量逐年升高，对城市供水水质安全形成了潜在威胁，本文根据岩溶水系统的水文地质条件和多年动态观测资料，采用Ｖｉｓｕａｌ　Ｍｏｄｆｌｏｗ　及ＭＴ３ＤＭＳ模块分别建立地下水水流和水质模型，模拟预测了规划水濠地未来开采条件下的溶质运移状况，划分了水ｉ啄地一衄，二衄和　准保护区。为加强岩溶水系统的环境管理和水源地水质防护提供了决策依据。　关键词：岩溶水系统　地下水　水泺地　保护区划分　中图分类号：ＴＶ　文献标识码：Ａ　文章编号：ｌ８７４一ｏ９８ｘ（２ｏｏ８）ｏ７（ｂ）一ｏｌ１０－０４　羊庄岩溶水系统位于山东省枣庄市和　滕州市交界处的羊庄盆地，是两市地下水　资源最丰富的地区。以往地下水资源的开　发利用，主要为工业和当地农业及农村用　水；随着枣庄市、滕州市城市规模不断扩　大和经济迅速发展，城市需水量日益递增，　为缓解城区供水紧张局面，分别规划在岩　溶水系统的龙山头和羊庄地段建设水源地　向两城区供给饮用水源。由于盆地内大部　分地区基岩裸露，地表水、地下水转换迅　速，降水入渗条件优越，地质环境相对脆　弱；地表污染物不断往地下水中渗透运移，　导致地下水的环境质量不断向变差方向发　展。开展饮用水水源地保护区划分，对于　加强岩溶水系统的环境保护与污染防治、　保障枣庄市和滕州市的饮水水质安全具有　重要的指导意义。　１岩溶水系统概况及地质环境条件　１．Ｉ岩溶水系统概况　羊庄岩溶水系统又称羊庄盆地岩溶水　系统，范围涉及枣庄市山亭区和滕卅Ｉ市的　Ｉ２个乡镇，总面积６５０ｋｌＴｔ　。盆地四周低山　环绕，中北部地区丘陵起伏，溶沟、溶槽．　溶洞和石芽等溶蚀地貌发育，并伴生有多　处岩溶泉群。总体地势东北高、西南低，下　游羊庄一西集一带分布着面积约Ｉ　１　０ｋｍ　的　Ｉ一第四系；２～奥陶系；３一寒武系下统；４一寒武系中上统；　平原区域，西南部的龙山头地段地势最低，　５一泰山群；６～隐伏地层代号；７一断裂；８一地层界线；　是岩溶水系统地下水与地表水的唯一排泄　９一隐伏地层界线；１　０一辛召断块；１　１～山亭断块；１　２一羊庄断块　出口。区内多年平均降水量７６１．５ｍｌＴｔ，主　图１　羊庄岩溶水系统地质略图　要地表水体为新薛河及其支流，干流长度　２３ｋｉｎ，河宽８０～１２０ｍ，河道迂回曲折，河床　成；羊庄断块则主要由寒武系上统和奥陶　组富水性弱，地下水资源贫乏。羊庄断块　基岩裸露，形成岩溶水与地表水交替循环　系组成（图１）。　子系统处在盆地下游，地形开阔平缓，地下　的外部通道。　１．２．２水文地质条件　网络状岩溶发育程度高、连通性好，储水　区内分布有４２０多个村庄．３个乡镇驻　羊庄岩溶水系统水文地质边界和羊庄　空间大；地下水不但可以得到该区降水的　地和一个规模较小的县级城市；社会经济　盆地的水文边界基本一致，南部．东南部、　入渗补给，上游两个子系统中的地下水和　以粮食作物和林果业为主，山亭城区一带　东北部和西北部±匀－由地表分水岭环绕，西　地表水也均向该区汇集，补给条件优越，岩　分布有造纸、食品加工等部分工业企业；　部受化石沟阻水断裂控制，西南部有唯一　溶水资源丰富。　社会、生态环境优良，近年来经济发展速　的地下水和地表水出口，地下水以河水基　在地貌、构造和地层分布的影响下，　度较快。　流和第四系潜流的方式排出区外　大气降　羊庄断块岩溶水子系统又具有相对完整的　１．２地质环境条件　水是地下水的主要补给来源，总体径流方　边界条件和一定的性特征，地下水具　１．２．１地质条件　向由东北往西南，自然状况下以补给河水、　有统一的径流场和基本完整的潜水一微承　羊庄盆地位于鲁西台隆的南翼边缘，　泉及潜流的形式排泄，目前主要通过人工　压水水面；发育有魏庄．后石湾．西石楼．　内部断裂构造发育，寒武．奥陶系碳酸盐　开采排泄；新薛河河床岩溶发育，河水与地　羊庄和龙山头五处岩溶富水地段（图２），钻　岩地层广泛分布。整个盆地被长龙断裂和　下水存在密切的补排转化关系。　孔单位涌水量均大干１　Ｏ００ｍ　／（ｄ・ｍ）；其中　曹王墓断裂分割成三大断块，由北往南依　整个岩溶水系统可划分为辛召断块、　前三处富水地段已经建设为工业供水水源　、次为辛召断块，山亭断块和羊庄断块　辛　山亭断块和羊庄断块三个子系统。前两个　地，后两处富水地段为目前规划开发的城　召断块主要由太古界太山群和寒武系下统　子系统位于中上游地区，地势高，地形起伏　市供水水源地。　组成‘；山亭断块主要由寒．武系中、上统组　变化大，地下水循环，蓄存条件差，含水岩　ｌ　ｌ　０科技创新导报Ｓｃｉｅｎｃｅ　ａｎｄ　Ｔｅｃｈｎｏｌｏｇｙ　Ｉｎｎｏｖａｔｉｏｎ　Ｈｅｒａｌｄ　维普资讯 http://www.cqvip.com 资源与环境　Ｓｃｉｅｎｃｅ　ａｎｄ　Ｔｅｃｈｎｏｌｏｇｙ　ｌｎｎｅｖ８　根据两市“十一五”及２０２０年供水规　划，枣庄市计划２０ｌ　０年以后开发龙山头水　源地向新城区供水，开采规模４万１ＴＩ　／ｄ；而　滕州市则在２００８年开发羊庄水源地向城区　供水，开采量２０ｌ０年以前４万１ＴＩ　／ｄ，２０１０年　以后增至８万ｍ　／ｄ。　３．２　水源地开采后的地下水流场预测　两规划水源地均位于羊庄断块岩溶水　子系统的下游地区，开采后的区域地下水　流场特征，根据水文地质条件采用Ｖｉ　Ｓｔｌ　ａｌ　ＭＯｄｆｌＯｗ４．０建立地下水水流模型进行预　测，计算范围以所在岩溶水子系统为基础，　面积约１　３４ｋｉｎ　。描述岩溶水渗流特征的偏　微分方程为：　妥（喀ｊ＋毒ｌ　—ｓ　南　ｔ　１　！　ｆ、＝＾　（　ｒ）　ｘ　ｖ《Ｑ　ｒ＝ｇ　ｒ，ｊｔ∈ｒ　，　ｏ　－７　ｌ　一－　２　■～３卜　；４　√　５　一９　６｝一０　７　Ｉ　８　●　９　１一径流补给边界；２　地表分水岭；３一隔水边界；４～径流排泄边界；　５一地质界线；６地层代号；７一富水地段；８一地下水流向；９一岩溶泉　图２　羊庄断块岩溶水子系统水文地质略图　式中：　一渗流区域；ｈ…地下水水位标高　（１ＴＩ）；　一含水屡渗透系数（ｍ／ｄ）｛￡一含水　层的源、汇项（１ＴＩ　／ｄ）；ｈ　一初始水位分布　（１ＴＩ）；　一潜水含水层在潜水面上的重力给　水度；厂一渗流区域的流量边界；　．．一边界　面法线方向的渗透系数（１ＴＩ／ｄ）；ｎ一边界面　的法线方向；日　厂边界单位面积上的流量　（ｍ／ｄ）。　在对计算区进行网格剖分和模型识　别、验证的基础上，将水源地规划开采量　及各源、汇项代入模型。以２００７年８月的　实测水位线作为初始流场预测出水源地开　采后的地下水稳定流场（图３）。预测结果显　示，两水源地开采后将导致地下水位的平　盘式降低，而各自形成的水位降落漏斗范　围则较小，区内地下水的总体流场特征不　会改变。　４水源地保护区划分　４．１划分方案　’　４．１．１保护区分级　羊庄岩溶水系统是一个开放型的地下　水系统，新薛河自上而下贯穿全区，河水与　各岩溶水子系统中地下水均存在着密切的　水力联系；水源地水质不但受所在岩溶水　子系统环境污染的影响，也同时受上游两　个岩溶水子系统中污染质运移、传递的威　胁。鉴于饮用水水质保护需要，将整个岩　溶水系统均划定为水源地保护区，分为一　级．二级和准保护区三个级别。　４．１．２保护区边界确定　级、二级保护区边界，分别以污染　质在ｌ００天和ｌ０００天时间内运移至水源地　的捕获带范围界线为基准进行划定　Ｊ，部分　边界可在保证水质的前提下，本着便于管　理和利于保护的原则结合具体环境条件进　一图３　水源地开采后的地下水预测流场　２地下水资源及水环境质量状况　２．１地下水资源　根据水文地质条件和多年动态观测资　料，分别采用水量均衡法　人渗参数模型　最优回归模型和有限元数值模型对岩溶水　系统地下水资源进行综合评价：岩溶水补　给资源量２４．１４万ＩＴＩ　／ｄ，开采资源量２４．Ｏ４　万ＩＴＩ　／ｄ…，储存量７９５００万ＩＴＩ　。目前实际开　采量１　１．６８万ＩＴＩ　／ｄ，则剩余开采资源量为　１２．３４万ＩＴＩ　／ｄ。　２．２水环境质量状况　参照２０００～２００７年多个时期的水质分　析资料，分别对代表性地段与河段的水体　环境质量进行评价：新薛河及其它地表水　的大部分化学指标含量符合（Ｇ　Ｂ　３　８　３　８—　２００２）中Ｉ类质量标准，部分化学指标含量　符合Ⅱ类或Ⅲ类质量标准；岩溶水的理化　指标和毒理学指标含量符合地下水（Ｇ　Ｂ／　Ｔ１４８４８～９３）中Ⅱ类或Ⅲ类质量标准，但溶　解性总固体　总硬度、氯氮及亚盐等　部分指标含量逐年上升，部分地点细菌学　指标多次严重超标，说明岩溶水已经受到　污染并呈现逐渐加重趋势，污染源主要为　行适当调整；准保护区边界划至岩溶水系　农村、农业面源污染物及上游城镇排放的　统的外围边界，即羊庄盆地的地表分水岭。　污水、废水　４．２划分技术方法　水源地一级、二级保护区边界采用水　３水源地开发利用规划及地下水流场预测　质模型来计算确定。根据地下水渗流特征　３．１水源地开发利用规划　和水源地开采影响范围，『二级保护区应基　科技创新导报Ｓｃｉｅｎｃｅ　ａｎｄ　Ｔｅｃｈｎｏｌｏｇｙ　Ｉｎｎｏｖａｔｉｏｎ　Ｈｅｒａｌｄ　ｌ　ｌ　ｌ　维普资讯 http://www.cqvip.com

ＳＣｉｅｎｃｅ　ａｎｄ　Ｔｅｃｈｎｏｌｏｇｙ　ｌｎｎｏｖａｔｉｏｎ　Ｈｅｒａｌｄ　Ｑ　Ｑ：！　资源与环境　本处在羊庄断块岩溶水子系统范围内，所　以可将前述水流模型的渗流区作为水质模　型的模拟区。　Ｖｉｓｕａｌ　Ｍｏｄｆｌｏｗ中的ＭＴ３ＤＭＳ模块是　基于有限差分数值方法求解地下水中污染　质运移的模拟软件，较全面地考虑了污染　物在地下水中的对流、弥散和化学反应等　过程。据其建立模拟区的平面二维溶质运　移模型为：　喀　Ｃ（ｘ　ｙｌ　，ｆＪ　Ｉ＿日　‘．（　ｊ，：‘　ｃｏ　：　ｊ　＝ｃ【　：ｆ　Ｊ　．　～：∈ｎ　一一　羊庄水源地　龙山头水源地　Ｉ＝Ｅ　，ｆ　２０　＿ｅ　篓　，，∽川　ｒ．Ｉ．　＝ＥＩ　ｆ　２　０　图４部分边界点溶质运移ｌ　００天模拟图　一　ｗ，　：Ｅ　ｒ、ｆ　１Ｉ．　式中：　一渗流区域；ｃ一污染组分浓度　（ｍｇ／Ｌ）；Ｕ　一二维空间上的地下水实际流　速（ｍ／ｄ）；Ｄ。．～水动力弥散系数张量（ｍ　／ｄ）；　有效孔隙度；鼋　一源的补给强度（ｍ／ｄ）；　Ｃ　一源的溶质浓度（ｍｇ／Ｌ）；鼋　一汇的排泄　强度（ｍ＇／ｄ）；Ｐ，、一含水层岩土的干容重；　一　生物降解系数；ｃ　一固体颗粒吸附的溶质　浓度（ｍｇ／Ｌ）；Ｒ一阻滞因子；ｃｎ（　，　，　一污　染组分的初始浓度（ｍｇ／Ｌ）；厂．一一类边　界；Ｃ．（Ｘ，　，ｚ，ｆ）一一类浓度边界值（ｍｇ／Ｌ）；　厂，一二类边界；　（Ｘ，　，Ｚ，ｔ）～二类边界值　（ｍｇ／Ｌ・ｍ／ｄ）；厂　一三类边界；ｑ．一边界上　单位面积的渗透率（ｍ／ｄ）；ｇ．（　，　，ｚ，ｆ）一三　类边界值（ｍｇ／Ｌ・ｍ／ｄ）　ＭＴ３ＤＭＳ提供了多种算法，溶质在区　内岩溶含水层中的运移以对流为主，则求　解对流项可选用欧拉一拉格朗日混合方法　并包括特征线法（ｆｏｒｗａｒｄ～ｔｒａｃｋｉｎｇ　ｍｅｔｈｏｄ　图５部分边界点溶质运移１　０Ｏ０天模拟图　０１。模拟的污染溶质选取污染负荷较　ｏｆ　ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ，Ｍ０Ｃ），求解矩阵方程采　率９＝０．用广义共轭梯度法（Ｇｅｎｅｒａｌｉｚｅｄ　ｃｏｎｊｕｇａｔｅ—　高、危害性较大的亚盐氮作为代表，其　ｇｒａｄｉｅｎｔ　ｓｏｌｖｅｒ　ｐａｃｋａｇｅ，ＧＣＧ）。　在Ⅲ类水标准中的含量不得超过０．０２ｍｇ／Ｌ，　ｍｇ／Ｌ，为　根据部分试验成果，并参照类似条件　目前地下水中的含量略低于０．０ｌ地区的参数值及经验值，给定地下水的弥　了简化计算，模拟时将现状含量作为其初　散系数ＤＬ＝３０ｍ　／ｄ，岩溶含水层的有效孔隙　始浓度，并设定为０，而将Ⅲ类水质标准的　含量定义为０．０１ｍｇ／Ｌ。　将水质模型在水流模型预测的地下水流　场环境中运行，模拟污染质在渗流场的不　同位置和方向往水源地扩散运移的状况，　污染质锋面（亚盐氮浓度为０．０ｌｍｇ／Ｌ）　运移ｌＯ０天（图４）和１　０００天（图５）刚好到达　水源地的试验点，就是沿该径流方向一级　和二级保护区的边界点。　４．３各级保护区范围　根据模拟计算结果和水质保护需要，　饮用水水源地保护区共划分为三级四个　区。一级保护区分为两个，即龙山头水源　地一级保护区和羊庄水源地一级保护区，　其范围分别处在每个水源地的周围。二级　保护区和准保护区为两个水源地的共同保　护区；其中二级保护区主要在羊庄断块岩　溶水子系统范围内，龙山头西部视松散层　的分布及与区外地下水的联系适当延伸至　区外；准保护区范围为二级保护区以外，．到　羊庄岩溶水系统边界之间的区域（表ｌ和图６）。　表１　饮用水水源地保护区划分表　保护　范嘲　ｌ龙山头水源地　级Ｉ　ｔ废水源地　村庄（个）　ｌＯ　耐积（ｋｍ　）　ｌ２．６２　ｌ６　８６　１５．２６　８８４４　级　１～饮用水水源地；２一保护区分区界线；３一一级保护区；　４～二级保护区；５准保护区　图６　饮用水水源地保护区划分图　准保护区　合　３ｌ】　４２３　５３３．６８　６５０　（下转１１　４页）　１　１　２　科技创新导报Ｓｃｉｅｎｃｅ　ａｎｄ　Ｔｅｃｈｎｏｌｏｇｙ　Ｉｎｎｏｖａｔｉｏｎ　Ｈｅｒａｌｄ　维普资讯 http://www.cqvip.com ！塑！　Ｑ：！　Ｓｃｉｅｎｃｅ　ａｎｄ　ＴｅｃｈｎＯｌＯｇｙ　ｌｎｎｏｖａｔｉｏｎ　Ｈｅｒａｌｄ　资源与环境　ｃｌａｍｓ（Ｓｉｌｉｑｕａ　ｐａｔｕｌａ）ｄｕｆｉｎｇｌ９９１　１９９３［Ｊ】．　Ｎａｔ．Ｔｏｘｉｎｓ，２：ｌ９７　２０５（１９９４ｂ）．　【ｌ３］Ｎｏｖａｃｚｅｋ，Ｉ．，ｅｔ　ａ１．Ｕｐｔａｋｅ，ｄｉｓｐｏｓｉｔｉｏｎ，　ａｎｄ　ｄｅｐｕｒａｔｉｏｎ　ｏｆ　ｄｏｍｏｉｃ　ａｃｉｄ　ｆｒｏｍ　有可能通过食物链进行传播，即使ＤＡ没有累　积到致死剂量，也会对踺康造成慢性影响。　然而，目前并无法预测到ＤＡ有何慢性影响。　对于ＤＡ的研究还有待于检测技术的进一步发　展。　Ａｎ　ｏｕｔｂｒｅａｋ　ｏｆ　ｔｏｘｉｃ　ｅｎｃｅｐｈａｌＯｐａｔｈｙ　ｃａｕｓｅｄ　ｂｙ　ｅａｔｉｎｇ　ｍｕｓｓｅｌｓ　ｃｏｎｔａｍｉｎａｔｅｄ　ｗｉｔｈ　ｄｏｍ０ｉｃ　ａｃｉｄ［Ｊ】．Ｎｅｗ　Ｅｎｇ１．Ｊ．Ｍｅｄ．，　３２２：ｌ７７５—１７８０（１９９０）　【６】Ｔｅｉｔｅｌｂａｕｍ，Ｊ．Ｓ．，Ｒ．Ｊ．Ｚａｔｏｒｒｅ，Ｓ．　参考文献　【ｌ】Ｐｅｒｌ，Ｔ．Ｍ．ｅｔ　ａ１．，ｌ９８７，Ｉｎｔｏｘｉｃａｔｉｏｎ　ｆｏｌ　ｌｏｗｉｎｇ　ｍｕｓｓｅｌ　ｉｎｇｅｓｔｉｏｎ　ｉｎ　Ｍｏｎｔｒｅａｌ　ｌｉｖｅ　ｂｌｕｅ　ｍｕｓｓｅｌｓ（Ｍｙｔｉｌｕｓ　ｅｄｕｌｉｓ）［Ｊ】．　Ｃａｒｐｅｎｔｅｒ，Ｄ．Ｇｅｎｄｒｏｎ，Ａ．Ｃ．Ｅｖａｎｓ．Ａ．　Ａｑｕａｔ．Ｔｏｘｉｃｏ１．，２ｌ：ｌ０３　ｌ１８（１９９１）．　Ｎｏｖａｃｚｅｋ，Ｉ．，ｅｔ　ａ１．Ｄｅｐｕｒａｔｉｏｎ　ｏｆ　ｄｏｍｏｉｃ　Ｇｊｅｄｄｅ，ａｎｄ　Ｎ．Ｒ．Ｃａｓｈｍａｎ．Ｎｅｕｒｏｌｏｇｉｃ　【１４］ｓｅｑｕｅｌａｅ　ｏｆ　ｄｏｍｏｉｃ　ａｃｉｄ　ｉｎｔｏｘｉｃａｔｉｏｎ　ａｃｉｄ　ｆｒｏｍ　ｌｉｖｅ　ｂｌｕｅ　ｍｕｓｓｅｌｓ（Ｍｙｔｉｌｕｓ　ｄｕｅ　ｔｏ　ｔｈｅ　ｉｎｇｅｓｔｉｏｎ　ｏｆ　ｃｏｎｔａｍｉｎａｔｅｄ　ｅｄｕｌｉｓ）［Ｊ】．Ｃａｎ．Ｊ．Ｆｉｓｈ．Ａｑｕａｔ．Ｓｃｉ．，４９：　３ｌ２　３ｌ８（１９９８）．　【Ｊ】．Ｃａｎ．Ｄｉｓ．Ｗｅｅｋ．Ｒｅｐ．，ｌ　３：２２４　２２６．　【２］Ｂａｔｅｓ，Ｓ．Ｓ，ｅｔ　ａ１．．１９８９，Ｐｅｎｎａｔｅ　ｄｉａｔｏｍ　Ｎｉｔｚｓｃｈｉａ　ｐｕｎｇｅｎｓ　ａｓ　ｔｈｅ　ｐｒｉｍａｒｙ　ｓｏｕｒｃｅ　ｍｕｓｓｅｌｓ［Ｊ】．Ｎ．Ｅｎｇ１．Ｊ．Ｍｅｄ．，３２２：１７８ｌ—　ｌ　７８７（１９９０）．　【７］Ｘｉ，Ｄ．，ｅｔ　ａ１．Ｄｏｍｏｉｃ　ａｃｉｄ　ｉｓ　ａ　ｐｏｔｅｎｔ　ｎｅｕｒｏｔｏｘｉｎ　ｔｏ　ｎｅｏｎａｔａｌ　ｒａｔｓ［Ｊ１．Ｎａｔ．　Ｔｏｘｉｎｓ．，５：７４　７９（１９９７）．　【８］Ｗｏｒｋ，Ｔ．Ｍ．，ｅｔ　ａ１．Ｄｏｍｏｉｃ　ａｃｉｄ　ｉｎｔｏｘｉ　ｃａｔｉｏｎ　ｏｆ　ｂｒｏｗｎ　ｐｅｌｉｃａｎｓ　ａｎｄ　ｃｏｒｍｏ　『ｌ　５１Ｄｒｕｍ，Ａ．Ｓ．，ｅｔ　ａ１．Ｄｏｍｏｉｃ　ａｃｉｄ　ｉｎ　ｔｈｅ　Ｐａｃｉｆｉｃ　ｒａｚｏｒ　ｃｌａｍ　ｌｉｑｕａ　ｐａｔｕｌａ（Ｄｉｘｏｎ．　ｏｆ　ｄｏｍｏｉｃ　ａｃｉｄ，ａ　ｔｏｘｉｎ　ｉｎ　ｓｈｅｌｌｆｉｓｈ　１７８９）［Ｊ】．Ｊ．Ｓｈｅｌｌｆｉｓｈ　Ｒｅｓ．，ｌ２：４４３　４５０　（１９９３）．　ｆｒｏｍ　ｅａｓｔｅｒｎ　Ｐｒｉｎｃｅ　Ｅｄｗａｒｄ　Ｉｓｌａｎｄ［Ｊ１．　Ｃａｎａｄａ，Ｃａｎ．Ｊ．Ｆ．Ｓｈ．Ａｑｕａｔｉｃ．Ｓｃｉ．，４６：　ｌ２０３　ｌ　２ｌ　５．　【ｌ６］Ｃａｍｐｂｅｌ１，Ｄ．Ａ．，ｅｔ　ａ１．Ａｍｎｅｓｉｃ　ｓｈｅｌｌ—　ｆｉｓｈ　ｐｏｉｓｏｎｉｎｇ　ｉｎ　ｔｈｅ　ｋｉｎｇ　ｓｃａｌｌｏｐ，　Ｐｅｃｔｅｎ　ｍａｘｉｍｕｓ，ｆｒｏｍ　ｔｈｅ　ｗｅｓｔ　ｃｏａｓｔ　『３］Ｔｏｄｄ．Ｅ．Ｃ．Ｄ．１９９０．Ａｍｎｅｓｉｃ　ｓｈｅｌｌｆｉｓｈ　ｐｏｉｓｏｎｉｎｇ　ａ　ｎｅｗ　ｓｅａｆｏｏｄ　ｔｏｘｉｎ　ｒａｎｔｓ　ｉｎ　Ｓａｎｔａ　Ｃｒｕｚ，Ｃａｌｉｆｏｒｎｉａ【Ｃ】．Ｉｎ：　Ｔｏｘｉｃ　ＰｈｖｔＯＤｌａｎｋｔＯｎ　Ｂｌｏｏｍｓ　ｉｎ　Ｓｅａ，　ｏｆ　Ｓｃｏｔｌａｎｄ［Ｊ】．Ｊ．Ｓｈｅｌｌｆｉｓｈ　Ｒｅｓ．，２０：７５　ｓｙｎｄｒｏｍｅ［Ｃ】．Ｉｎ：Ｇｒａｎｅｌｉ，Ｅ．，Ｓｕｎｄｓｔｒｏｍ，　Ｂ．，Ｅｄｌｅｒ，Ｌ．＆Ａｎｄｅｒｓｏｎ，Ｄ．Ｍ．（ｅｄｓ），　Ｔｏｘｉｃ　Ｍａｒｉｎｅ　ＰｈｙｔＯｐｌａｎｋｔＯｎ　Ｐｒｏ～　ＰＰ．６４３　６４９．（Ｓｍａｙｄａ，Ｔ．Ｊ．ａｎｄ　Ｙ．　８４（２００１）．　Ｓｈｉｍｉｚｕ，Ｅｄｓ）．Ａｍｓｔｅｒｄａｍ：Ｅｌｓｅｖｉｅｒ　（１９９３）．　ｃｅｅｄｉｎｇｓ　ｏｆ　ｔｈｅ　ｆｏｕｒｔｈ　ｉｎｔｅｒｎａｔｉＯｎａｌ　ｃｏｎｆｅｒｅｎｃｅ　ｏｎ　Ｔｏｘｉｃ　Ｍａ　ｒｉｎｅ　【９１Ｓｈａｗ，Ｂ．Ａ．，ｅｔ　ａ１．Ｆｅｅｄｉｎｇ　ｄｅｔｅｒｒｅｎｔ　ａｎｄ　ｔｏｘｉｃｉｔｙ　ｅｆｆｅｃｔｓ　ｏｆ　ａＤＯｆｕｃＯｘａｎｔｈｉｎＯｉｄｓ　ａｎｄ　ｐｈｙｃｏｔｏｘｉｎｓ　ｏｎ　ａ　ｍａｒｉｎｅ　ｃｏｐｏｐｏｄ　Ｐｈｙｔｏｐｌａｎｋｔｏｎ，ｈｅｌｄ　Ｊｕｎｅ　２６　３０，ｌ　９８９，　ｉｎ　Ｌｕｒｉｄ，Ｓｗｅｄｅｎ．Ｅｌｓｅｖｉｅｒ　Ｓｃｉｅｒｉｃｅ　Ｐｕｂ　ｌｉｓｈｉｎｇ　Ｃｏ．，Ｉｎｃ．，Ｎｅｗ　Ｙｏｒｋ．ＰＰ．５０４　５０８．　（Ｔｉｇｒｉｏｐｕｓ　ｃａｌｉｆｏｒｎｉｃｕｓ）【Ｊ】．Ｍａｒ．Ｂｉｏ１．，　ｌ　２８：２７３　２８０（１９９７）．　【ｌ０］Ｊｏｎｅｓ，Ｔ．Ｏ．，ｅｔ　ａ１．Ｈａｅｍｏｃｙｔｅ　ｃｈａｎｇｅｓ　ｉｎ　ｔｈｅ　Ｐａｃｉｆｉｃ　ｏｙｓｔｅｒ，Ｃｒａｓｓｏｓｔｒｅａ　ｇｉｇａｓ，　ｃａｕｓｅｄ　ｂｙ　ｅｘｐｏｓｕｒｅ　ｔｏ　ｄｏｍｏｉｃ　ａｃｉｄ　ｉｎ　ｔｈｅ　ｄｉａｔｏｍ　Ｐｓｅｕｄｏｎｉｔｚｓｃｈｉａ　ｐｕｎｇｅｎｓ　ｆ．　【４］Ｄｅｂｏｎｎｅｌ，Ｇ．Ｌ．，Ｌ．Ｂｅａｕｃｈｅｓｎｅ，ａｎｄ　Ｃ．Ｄｅ　ＭｏｎｔｉｇｎＹ．Ｄｏｍｏｉｃ　ａｃｉｄ，ｔｈｅ　ａｌｉｅｇｅ　“ｍｕｓｓｅｌ　ｔｏｘｉｎ．”ｍｉｇｈｔ　ｐｒｏｄｕｃｅ　ｉｔｓ　ｎｅｕ　ｒｏｔｏｘｉｃ　ｅｆｆｅｃｔ　ｔｈｒｏｕｇｈ　ｋａｉｎｉｔｅ　ｒｅｃｅｐ——　ｔｏｒ　ａｃｔｉｖａｔｉｏｎ：ａｎ　ｅｌｅｃｔｒｏｐｈｙｓｉｏｌｏｇｉｃａｌ　ｍｕｌｔｉｓｅｒｉｅｓ【Ｊ】．Ｔｏｘｉｃａｎ，３３：３４７　３５３　（１９９５ｂ）　ｓｔｕｄｙ　ｉｎ　ｔｈｅ　ｄｏｒｓａｌ　ｈｉｐｐｏｃａｍｐｕｓ［Ｊ】．　Ｃａｎ．Ｊ．Ｐｈｙｓｉｌ．Ｐｈａｒｍａｃｏ１．，６７：２９　３３　【ｌ　ｌ】Ｂａｔｅｓ，Ｓ．Ｓ．Ｄｏｍｏｉｃ　ａｃｉｄ　ｐｒｏｄｕｃｉｎｇ　ｄｉａｔｏｍｓ：ａｎｏｔｈｅｒ　ｇｅｎｕｓ　ａｄｄｅｄ！【Ｊ】．Ｊ．　（１９８９）　【５］Ｐｅｒｌ，Ｔ．Ｍ．，Ｌ．Ｂｅｄａｒｄ，Ｔ．Ｋｏｓａｔｓｋｙ，Ｊ．Ｃ．　Ｈｏｃｋｉｎ，Ｅ．Ｃ．Ｄ．Ｔｏｄｄ，ａｎｄ　Ｒ．Ｓ．Ｒｅｍｉｓ．　ｐｈｙｃｏ１．，３６：９７８　９８５（２０００）．　【ｌ　２】Ｗｅｋｅｌｌ，Ｊ．Ｃ．，ｅｔ　ａ１．Ｏｃｃｕｒｒｅｎｃｅ　ｏｆ　ｄｏｍｏｉｃ　ａｃｉｄ　ｉｎ　Ｗａｓｈ．ｎｇｔｏｎ　ｓｔａｔｅ　ｒａｚｏｒ　（上接１　ｔ　２页）　及时处理。　５保护区环境管理措施与水质防护对策　５．１环境管理措施　（１）一级保护区。实行圈网管理或在边　界处设立界桩，警示作业内容；禁止建　设与取水设施无关的建筑物，可能对　地下水产生污染的管．渠通过。　（２）二级保护区及准保护区。在重要性　边界点设立界碑，明确保护范围和保护要　求；禁止建设严重污染企业和废弃物堆放　场，特殊需要的应经环保部门批准，并采取　严格的防污染措施。　（３）保持新薛河及其它地表水体水质不　低于（ＧＢ３８３８～２００２）Ⅲ类标准，严禁在重　要补给区开山取石和开采矿产，禁止利用　污水灌溉农田，施用有毒、高残留化　肥、农药，大力推进植树造林和退耕还林　工程。　５．２水质防护对策　（１）在饮用水水源地、岩溶水系统不同　功能区和新薛河不同河段布设水质监测　点，定期采集水样分析测试，全面掌控区内　水体环境质量状况；发现问题，查明原因，　（２）加强对上游城镇及工业污水、废水　的处理效果监测，建设农村人畜排泄物集　纳．处理工程系统，确保实现达标排放；开　展各类建设项目环境影响评价，杜绝产生　新的污染源；实施一级保护区村庄、与水　源地无关建筑的搬迁规划，改善水源地周　围生态环境。　（３）大力宣传保护饮用水水源安全的重　参考文献　要意义，提高群众的水质防护意识；制定饮　【１】李传谟，康凤新．岩溶水资源及增源增采　用水水源地水质安全应急预案，及时控制　模型【Ｍ】．济南：山东科学技术出版社，　和处理各类水质污染事件的发生。　ｌ９９９。８０　８９．　素，对潜水或微承压水水源地来说，保护区　范围应尽量涉及整个水文地质单元。　（３）水质模拟模型是目前划分大中型地　下水水源地保护区的最先进方法，计算结　果具有较高可信度，但计算时要使其在水　源地开采后形成的稳定流场或一定开采周　期后的非稳定动力场条件下运行。　６结语　（１）羊庄岩溶水系统是一个边界条件完　整、功能区清晰、蓄存条件优良、三水转　化显著的地下水系统，划分饮用水水源地　保护区，对于实现饮用水水源的科学化管　理，保障枣庄市和滕州市的饮水水质安全　具有重要意义。　（２）地下水水源地保护区划分应充分考　虑含水层的开放性、水源地开采影响范围　以及地下水循环交替能力等水文地质因　【２】Ｈ　Ｊ／Ｔ３３８～２００７，饮用水水源保护区划　分技术规范【ｓ】．　１　１　４　科技创新导报Ｓｃｉｅｎｃｅ　ａｎｄ　Ｔｅｃｈｎｏｌｏｇｙ　Ｉｎｎｏｖａｔｉｏｎ　Ｈｅｒａｌｄ