地下水模型Feflow介绍ppt课件.ppt

TrainingGroundwatermodelingusingFEFLOW5 3 FEFLOW FiniteElementsubsurfaceFLOWsystem isoneofthemostoutstandingsoftwarepackagesforthesimulationofflowandmassand orheattransportprocessesinsaturatedandunsaturatedporousmedia History Functionalities 2 Dhorizontalorverticalor3 D1 D 2 Ddiscreteelementsforfracture tube channelflowFlow massand orheattransportDensitydependentflowSaturated variablysaturated unsaturatedflowandtransportGISinterfacefordataimportandexportSeparatetoolsforgeoreferencingimages FEMAP producingplots FEPLOT XPLOT andfittingparametersforunsaturatedflow sorption FE LM Openprogramminginterface Licensingandmaintenance LicensedviahardlockandlicensemanagersoftwareLicensetypes F2 2DflowFM22DflowandmasstransportF32D 3DflowFM32D 3DflowandmasstransportFMH32D 3Dflow massandheattransportFMHESEducationstarterkit 5limitedlicenses foruniversitiesonlySingleseatornetworklicensesOptionalmaintenancecontracthotline supportservicefreeupdateservice DemonstrationExercise 48p User sManual 168p ReferenceManual 278p WhitePapersVol I 360p WhitePapersVol II 112p WhitePapersVol III 42p WhitePapersVol IV 122p OnlineHelp AllmanualsareavailableprintedaswellasontheFEFLOWCD AdobeAcrobatfiles pdf Documentation GraphicalUserInterface Messagebar Workingwindow Infoboxes Shellmenu Structure FEFLOW Shell MeshEditor MeshGenerator ProblemEditor SimulatorKernel Postprocessor InterfaceManager IFM DefaultDirectoryStructure Projectdirectory femdata supermsh import export results plot Finiteelementproblemfiles fem Superelementmeshfiles smh Import exportfiles shp trp lin pow Resultsfiles dar dac Exportedplotdata shp lin trp pnt ASCII trp pnt lin ply pow dat txt rai smh dar Binary shp dxf tif dbf FEFLOWModelfilescanbesavedinBinaryorASCIIformat femfiniteelementmodelfiles dacmodelresultsfile FEFLOWModelfilessavedinASCIIformat smhsuperelementfiles darreducedmodelresultsfiles FEFLOWimportandexport Trainingprogram Fr 23rdMarch2007 8 30IntroductionbyProf Kaden卡登教授9 00FEFLOW地下水计算机模型特点 FEFLOWfeatures 10 00休息 Break 10 15FEFLOW操作与建模 SetupgroundwatermodelwithFEFLOW 12 00午饭 BreakforLunch 13 00FEFLOW操作 StepbystepusingFEFLOW 14 00休息 Break 14 15FEFLOW建立二维地下水模型 Stepbystepsetup2Dgroundwatermodel 16 30第一天培训结束 Trainingprogram Sa 24thMarch2007 8 30建立二维地下水模型 Stepbystepsetup2Dgroundwatermodel 10 00休息 Break 10 15建立三维地下水模型 Stepbystepsetup3Dgroundwatermodel 12 00午饭 BreakforLunch 13 00建立溶质迁移模型 Stepbystepsetupmasstransportmodel 14 00休息 Break 14 15答疑讨论16 30第二天培训结束 FEFLOWTrainingExercise骆俊峰柏林WASY水资源规划及系统研究所 23rd 25thMar 2007 JinanChina 中国济南 2007年三月23 25日 Modelbuilding Analysis Groundwatermodeling Problem InitialconditionsBoundaryconditionsMaterialproperties ModelboundariesMeshing GriddingTemporaldiscretization Simulation ComparisonwithmeasurementdataCheckforplausibility Target GIS Problemspecification DefinemodelareaDatapreparationandevaluationTopography elevation landuse riveretc GeologyHydrology climatedataGroundwaterextraction artificialrechargeObservationdateforgroundwater riverconsideredContaminantdata2D or3DModeling Constructhydrogeol Model Evaluateavailablegeol DataIdentifyandspecifyaquifersystemunconfinedorconfinedaquiferonelayerormulti layersModeldefinetop bottomelevationforeachlayerconsidered unit mormNN Identifyhydraulicconductivityformodellayer unit 10 4m s Porosity unit 1 e g 0 2 20 Naturalgroundwaterrecharge unit 10 4m d 模型区概况Overviewofthemodelarea 垃圾堆放场 污水处理厂 姆格尔湖 抽水井 农田 地表水 城区 森林 图列 Modelscenarioflowmodel 地表径流 outflow GW recharge GW flowdirection GW extraction N 入流边界 出流边界 地下水流向 地下水补给 地下水抽水井 GW recharge Source10 4m d Wells1000m d Fixedhead32 0m 1stBC Fixedhead46 0m 1stBC 初始水头 demo head ini trp 渗流系数 conduc2d trp 含水层底板标高 bot san1 trp 含水层顶板标高 demo relief trp 含水层顶板Aquifertop 含水层底板Aquiferbottom 渗流系数Conductivitym s 定水头边界 1stBC 模型数据Modeldata 地下水补给 定水头边界 1stBC 地下水抽水井 2D地下水流动模型2DModelscenarioflowmodel Stepsto2Dmodelbuilding Supermeshdesign 模型区初步圈定划分Meshgeneration 生成有限单元网格Modeltypedefine 定义模型类型Flowdata 设置流动模型参数Flowinitials设置初始水头 demo head ini trpFlowboundaries设置边界条件Flowmaterials设置模型流场介质水动力参数Transportdata TransportinitialsTransportboundariesTransportMaterialsReferencedata 2D地下水流动模型建模主要步骤 模型区离散化模型区初步圈定划分 Supermesh 有限单元网格 Finiteelement 地下水模型建模 底图文件 model area lin demo wells pnt mass src lin 文件夹 WASY FEFLOW5 3 demo excercise Meshgeneration Transportmapping TMesh Triangle interface Advancingfront 2Dmodelprojections Horizontalprojection Verticalprojection Axisymmetricprojection 2Dtransientsaturatedflowmodel butwithconstantboundarycondition 定义地下水模型类 2D非稳态饱和流动模型 但边界条件为常数 即不随时间变化 饱和流 非饱和流 非稳态流 定义模型类2D平面潜水流动模型 饱和流 流动问题 非稳态流 平面潜水 设置时间控制Temporal controldata 自动时间步长 初始步长 终止时间 初始时间 设置边界条件1st类边界条件定水头边界 单位 m 2nd类边界条件 单位 2D为m d 3D为m d 出流边界设为 值入流边界设为 值3rd类边界条件复合边界 水头 单位 m transferrateor leakagefactor 10 4m d 4th类边界条件 单井 边界条件 单位 m d 抽水井设为 值注水井设为 值AlltheBCcanbedefinedasconstantortime variable andcanbesetwithconstraints 地下水模型建模 设置流动模型参数 初始条件 地下水模型建模 设置初始条件 demo head ini trp 地下水模型建模 设置流动模型参数 初始条件 BuildComputerModel 设置流动模型参数 边界条件 定水头边界抽水井边界 GW recharge Source10 4m d Wells1000m d Fixedhead32 0m 1stBC Fixedhead46 0m 1stBC 初始水头 demo head ini trp 渗流系数 conduc2d trp 含水层底板标高 bot san1 trp 含水层顶板标高 demo relief trp 含水层顶板Aquifertop 含水层底板Aquiferbottom 渗流系数Conductivitym s 定水头边界 1stBC 模型数据Modeldata 地下水补给 定水头边界 1stBC 地下水抽水井 2D地下水流动模型2DModelscenarioflowmodel Boundaryconditionsflow Note ByfluxBC valuesmeanoutflowfluxes whereas valuescorrespondinfluxesBywellBC valuesforpumpingwellsand valuesforinjectionwells 水流模型边界条件类型 3rdKindBoundaryCondition CauchyType q hOW hGW q Infiltrations Exfiltrationsrate m d Transferrate Leakage Factor 10 41 d hOWWatertableinRiverhGWWatertableinGroundwaterForTransferrate kf dkfconductivityofcolmationbedlayerdthicknessofcolmationbedlayer 3rdKindBoundaryCondition CauchyType ConstraintConditions Example Wellwithminimumhydraulicheadconstraint Qbc m d hmin m Q2 m d Qbc h m hmin Boundarycondition 4thkind Constraintcondition Resultingflux Boundarycondition 1stkind Note thealgebraicsignsforconstraintconditionsarecontrarytotheboundaryconditions i e influxesarepositive outfluxesnegative t0 0 t1 t0 输入 内插流场模型参数 渗流系数k值 conduc2d trp含水层底板标高 bot san1 trp含水层顶板标高 demo relief trpPorosity地下水补给 蒸发量GWrecharge 地下水模型建模 设置流场参数 地下水模型建模 输入 内插流场模型参数 渗流系数k值Conductivity 10 4m s 输入 内插流场模型参数 含水层底板标高 m 输入 内插流场模型参数 含水层顶板标高 m 输入 内插流场模型参数 地下水补给蒸发 4 输入 内插流场模型参数 地下水补给蒸发量 10 4m d 地下水补给蒸发量4 10 4m d 146mm a Runsimulation Saveresultsin dacfile 2 Startsimulation 3 Viewresults 3Dmodel Slice1 Slice2 Slice3 Slice4 Layer1 sand Layer2 clay Layer3 sand Fromgeologicallayeringtoa3Dmodel demo relief trp Initialhead demo head ini trp bot san1 trp bot san2 trp bot clay trp Conductivity condu2d trp Conductivity 10 6 10 4m s Conductivity 10 10 4m s Porosity 0 1 Porosity 0 01 Porosity 0 2 Groundwaterrecharge recharge normal year ply Stepsto3Dmodelbuilding SupermeshdesignMeshgeneration3DmodelbuildingLayerconstructSliceelevationspecification3Dmodel 2DsectionvieworcheckModeltypedefineFlowdata Flowinitials sliceFlowboundaries sliceFlowmaterials layerTransportdata Transportinitials sliceTransportboundaries sliceTransportMaterials layerReferencedataObservationpoints sliceObservationpointgroups sliceReferencedistribution 3Dmodelbuilding 1 Dimensiondefinitionaftermeshing 3Dmodelbuilding 2 3Dlayerconfiguration 3Dmodelbuilding 2 3Dlayerconfiguration 3Dmodelbuilding 2 3Dlayerconfiguration Slice2 3Dmodelbuilding 2 3Dlayerconfiguration 3Dmodelbuilding 2 3Dlayerconfiguration Slice3 3Dmodelbuilding 2 3Dlayerconfiguration 3Dmodelbuilding 2 3Dlayerconfiguration Slice4 3Dmodelbuilding 2 3Dlayerconfiguration Modeltypedefinition Modeltypedefinition Setflowdatatolayer1 aquifer1 Import InterpolationofKx value condc2d trp CopyKxtoKy KzSetporosityglobalwith0 1Setgroundwaterrechargewith join funktion recharge normal year ply Setflowdatatolayer2 aquitard SetKx valueglobalwith1e 6CopyKxtoKy KzSetporosityglobalwith0 01 Setflowdatatolayer3 aquifer2 SetKx valueglobalwith10CopyKxtoKy KzSetporosityglobalwith0 2 Setreferencedataallowingtocomparesimulatedandmeasuredvalues Setobservationpoints slice related forobservationwellsSetobservationpointgroups slice related forwaterbalanceSetreferencedata slice related forknownareadataSetsectionfor2Dcrosssectionshow Setreferencedata Forheadobservation points Forfluxobservation nodes e g asaerialreferenceheaddistribution Slice related Modelcalibration verificationbasedon observedvaluesfromobservationwellsWaterbalanceanalysisknownIsolinesofgroundwaterlevelsComparisonofknownwithmodeledvalues Postprocessor LoadaFEFLOW resultfile dac Browseresultssaved Viewresultsselectedtime Waterbalanceanalysis Spatial temporalanalysis Saveresultsas dar Create femfrom dac Modelscenariotransportmodel inflow outflow GW recharge GW flowdirection GW extraction contaminate Contaminanttransportmodeling Processes AdvectionMoleculardiffusionDispersionNumerical artificial dispersion Processes AdsorptionDecayReaction availableupFEFLOW5 2 Conservativetransportconsidering Nonconservativetransportconsideringadditional Contaminanttransport AdvectiveTransport RealTransport Extendaflowmodeltoatransportmodel Transportdata Transportinitials sliceTransportboundaries sliceTransportMaterials layer Onespeciestransport concentration Boundaryconditionsmasstransport Initialconcentrationformasstransport TransportMaterials Porosity unit1 dependingonaquiferandaquitardtypeDispersivity unitm Langitudinaldispersivity alongmassfluxdirection TranservedispersivitynormaltomassfluxdirectionSorption unit1 dependingonaquiferandaquitardtypeDecay 10 41 s Reaction availabletomultispeciesonly Contaminanttransportmodeling Note Redcolouredarerequiredminimally Problems OscillationscanoccurduetotemporalandspatialdiscretizationSolutionforoscillations finemeshing Pecletnumber 2 dampingtechniquesinspace upwind andtime FE BE NOTE Alldampingtechniquesaddnumerical artificial dispersion Contaminanttransportmodeling Numericalmethodslowestnumericaldispersion butloweststabilityagainstoscillationsnoupwind GFEM AB TRnoupwind GFEM FE BEshockcapture SC AB TRshockcapture SC FE BEstre