druckerprager.docx
《druckerprager.docx》由会员分享，可在线阅读，更多相关《druckerprager.docx（37页珍藏版）》请在冰豆网上搜索。
druckerprager
exa_hydfracture10.1.5 Hydraulicallyinducedfractureinawellbore
液压诱导井眼破裂
unitss,kNandKPaandm
SolidSection,elset=lower_shale,material=Mat_lower_shale
**
*Material,name=Mat_lower_shale下层页岩
*Elastic
0.625e+07,0.19
*DruckerPrager
30.,1.0,30.
*DruckerPragerHardening
32.00e+03,0.00
33.000e+03,0.005
33.500e+03,0.08
*Permeability,specific=9.8
7.6720e-09,0.17
7.6720e-09,0.19
8.12200e-09,0.24
9.0500e-09,0.30
1.8800e-08,0.36
**
*SolidSection,elset=upper_shale,material=Mat_upper_shale
**
*Material,name=Mat_upper_shale上层页岩
*Elastic
0.4137e+07,0.15
*DruckerPrager
28.,1.0,28.
*DruckerPragerHardening
26.207e+03,0.00
26.500e+03,0.005
27.00e+03,0.08
*Permeability,specific=9.8
9.6720e-09,0.17
9.6720e-09,0.19
1.0200e-08,0.24
2.0500e-08,0.28
4.8800e-08,0.32
**
*SolidSection,elset=target_formation,material=Mat_target_formation
**
*Material,name=Mat_target_formation目标地层
*Elastic
1.2e+07,0.22
*DruckerPrager
36.,0.95,36.
*DruckerPragerHardening
37.930e+3,0.
38.20e+3,0.005
38.34e+3,0.02
38.40e+3,0.08
*Permeability,specific=9.8？
？
？
？
？
other》porefluid井眼流动》permeability表示渗透性。
。
。
。
。
0.560e-08,0.2
1.000e-08,0.22
4.836e-08,0.25
4.880e-08,0.29
4.980e-08,0.31
5.280e-08,0.33
5.800e-08,0.35
6.525e-08,0.40
8.100e-08,0.48
*MembraneSection,elset=casing,material=casing
0.01,
*Material,name=casing
*Elastic
2.e+06,0.3
**
*InitialConditions,type=Ratio,user
*InitialConditions,type=Porepressure,user
*InitialConditions,type=Stress,user
**
*elset,elset=initopen,generate
5183,5186
*INITIALCONDITIONS,TYPE=INITIALGAP
initopen
*nset,nset=iniflm,generate
1000076,1000078
*nset,nset=inifle
1000075,1000079
*nset,nset=inifl
iniflm,inifle
*Nset,nset=target_zone
exa_hydfracture10.1.5 Hydraulicallyinducedfractureinawellbore
*Elset,elset=__PickedSurf260_S4,generate
5383,5440,3
*Surface,type=ELEMENT,name=_PickedSurf260
__PickedSurf260_S5,S5
**
*COHESIVESECTION,ELSET=upper_inter,MATERIAL=upper_inter_mat,RESPONSE=TRACTIONSEPARATION,
THICKNESS=SPECIFIED,CONTROLS=VISCO_upper
*SECTIONCONTROLS,NAME=VISCO_upper,VISCOSITY=0.01
*MATERIAL,NAME=upper_inter_mat
*ELASTIC,TYPE=TRACTION
,,
*DAMAGEINITIATION,CRITERION=QUADS
,,
*DAMAGEEVOLUTION,TYPE=ENERGY,MIXEDMODEBEHAVIOR=BK,POWER=
,,
*GAPFLOW
*FLUIDLEAKOFF,USER
,
**
*COHESIVESECTION,ELSET=lower_inter,MATERIAL=lower_inter_mat,RESPONSE=TRACTIONSEPARATION,
THICKNESS=SPECIFIED,CONTROLS=VISCO_lower
*SECTIONCONTROLS,NAME=VISCO_lower,VISCOSITY=0.01
*MATERIAL,NAME=lower_inter_mat
*ELASTIC,TYPE=TRACTION
,,
*DAMAGEINITIATION,CRITERION=QUADS
,,
*DAMAGEEVOLUTION,TYPE=ENERGY,MIXEDMODEBEHAVIOR=BK,POWER=
,,
*GAPFLOW
*FLUIDLEAKOFF,USER
,
**
*COHESIVESECTION,ELSET=middle_inter,MATERIAL=middle_inter_mat,RESPONSE=TRACTIONSEPARATION,
THICKNESS=SPECIFIED,CONTROLS=VISCO_middle
*SECTIONCONTROLS,NAME=VISCO_middle,VISCOSITY=0.01
*MATERIAL,NAME=middle_inter_mat
*ELASTIC,TYPE=TRACTION
,,
*DAMAGEINITIATION,CRITERION=QUADS
,,
*DAMAGEEVOLUTION,TYPE=ENERGY,MIXEDMODEBEHAVIOR=BK,POWER=
,,
*GAPFLOW
*FLUIDLEAKOFF,USER
,
**
*SolidSection,elset=lower_shale,material=Mat_lower_shale
**
*Material,name=Mat_lower_shale
*Elastic
0.625e+07,0.19
*DruckerPrager
30.,1.0,30.
*DruckerPragerHardening
32.00e+03,0.00
33.000e+03,0.005
33.500e+03,0.08
*Permeability,specific=9.8
7.6720e-09,0.17
7.6720e-09,0.19
8.12200e-09,0.24
9.0500e-09,0.30
1.8800e-08,0.36
**
*SolidSection,elset=upper_shale,material=Mat_upper_shale
**
*Material,name=Mat_upper_shale
*Elastic
0.4137e+07,0.15
*DruckerPrager
28.,1.0,28.
*DruckerPragerHardening
26.207e+03,0.00
26.500e+03,0.005
27.00e+03,0.08
*Permeability,specific=9.8
9.6720e-09,0.17
9.6720e-09,0.19
1.0200e-08,0.24
2.0500e-08,0.28
4.8800e-08,0.32
**
*SolidSection,elset=target_formation,material=Mat_target_formation
**
*Material,name=Mat_target_formation
*Elastic
1.2e+07,0.22
*DruckerPrager
36.,0.95,36.
*DruckerPragerHardening
37.930e+3,0.
38.20e+3,0.005
38.34e+3,0.02
38.40e+3,0.08
*Permeability,specific=9.8
0.560e-08,0.2
1.000e-08,0.22
4.836e-08,0.25
4.880e-08,0.29
4.980e-08,0.31
5.280e-08,0.33
5.800e-08,0.35
6.525e-08,0.40
8.100e-08,0.48
*MembraneSection,elset=casing,material=casing
0.01,
*Material,name=casing
*Elastic
2.e+06,0.3
**
*InitialConditions,type=Ratio,user
*InitialConditions,type=Porepressure,user
*InitialConditions,type=Stress,user
**
*elset,elset=initopen,generate
5183,5186
*INITIALCONDITIONS,TYPE=INITIALGAP
initopen
*nset,nset=iniflm,generate
1000076,1000078
*nset,nset=inifle
1000075,1000079
*nset,nset=inifl
iniflm,inifle
*Nset,nset=target_zone
5,6,7,8,9,10,25,28,65,66,67,68,69,70,71,72
73,74,75,76,77,78,79,80,81,82,83,84,85,86,87,88
89,90,91,92,93,94,95,96,173,174,175,176,177,178,179,363
364,365,366,367,368,369,457,458,459,460,461,462,463,464,465,466
467,468,469,470,471,472,473,474,616,617,618,619,620,621,622,623
624,625,626,627,628,629,630,631,632,633,634,635,636,637,638,639
640,641,642,643,644,645,646,647,648,649,650,651,652,653,654,655
656,657,658,659,660,661,662,663,664,665,666,667,668,669,670,671
672,673,674,675,676,677,678,2264,2265,2266,2267,2268,2269,2270,2271,2272
2273,2274,2275,2276,2277,2278,2279,2280,2281,2282,2283,2284,2285,2286,2287,2288
2289,2290,2291,2292,2293,2294,2295,2296,2297,2298,2299,2300,2301,2302,2303,2304
2305,2306,2307,2308,2309,2310,2311,2312,2313,2314,2315,2316,2317,2318,2319,2320
2321,2322,2323,2324,2325,2326
iniflm,
inifle,
*elset,elset=one
5184
*elset,elset=QA
3113,3105,3097,3089
5185,5186,5187,5188
*nset,nset=QA
1001015,1001016,1001017,1001018
**
**BOUNDARYCONDITIONS
**
**Name:
BC-1Type:
Symmetry/Antisymmetry/Encastre
*Boundary
_PickedSet4,XSYMM
**Name:
BC-2Type:
Displacement/Rotation
*Boundary
_PickedSet5,1,1
_PickedSet5,2,2
**Name:
BC-3Type:
Displacement/Rotation
*Boundary
_PickedSet6,3,3
**
*Amplitude,name=volumerate
0.0,0.0,200.0,-1.0
**----------------------------------------------------------------
**
**STEP:
Step-1
**
*Step,name=Step-1,nlgeom=YES,unsymm=YES
*Geostatic
**
**LOADS
**
**Name:
Load-1Type:
Pressure
*Dsload
_PickedSurf260,P,42000.
**Name:
Load-2Type:
Bodyforce
*Dload
_PickedSet8,BZ,-20.
**Name:
Load-3Type:
Pressure
*Dsload
well_bore,PNU,1.
**
*Boundary,user
**_PickedSet5,8,8,1
TOP,8,8,1
BOT,8,8,1
**
**OUTPUTREQUESTS
**
*Restart,write,frequency=0
**
**FIELDOUTPUT:
F-Output-1
**
**MAKECHANGESHERE
**usethevariable=PRESELECTversionofthe*outputkeywordline
***Output,field,variable=PRESELECT
*Output,field
*nodeoutput
POR,
**
**HISTORYOUTPUT:
H-Output-1
**
*Output,history,variable=PRESELECT
*EndStep
**----------------------------------------------------------------
**
**STEP:
PumpStage
**
*Step,name="PumpStage",nlgeom=YES,inc=200,unsymm=YES
*Soils,consolidation,end=PERIOD,utol=500000.
**MAKECHANGESHERE
**Usethe1200sdatalineinplaceofthe50sdataline
**1.,50.,0.1,30.,
1.,1200.,0.01,30.,
**
**Name:
BC-1Type:
Symmetry/Antisymmetry/Encastre
*Boundary,op=new
_PickedSet4,XSYMM
**Name:
BC-2Type:
Displacement/Rotation
*Boundary,op=new
_PickedSet5,1,1
_PickedSet5,2,2
**Name:
BC-3Type:
Displacement/Rotation
*Boundary,op=new
_PickedSet6,3,3
**
*Boundary,user,op=new
TOP,8,8,1
BOT,8,8,1
**Controls,PARAM=FIELD,FIELD=PORE
**,,1.0e-3
**
*cflow,amplitude=volumerate
iniflm,,
inifle,,
**
**FIELDOUTPUT:
F-Output-1
**
*Output,field
**MAKECHANGESHERE
**Removethenset=QA,andelset=QAparametersettings
**NodeOutput,nset=QA
*NodeOutput
POR,U
**elementoutput,elset=QA
*elementoutput
S,LE,PE
GFVR
SDEG
DMICRT
PFOPEN
LEAKVRT
LEAKVRB
ALEAKVRT
ALEAKVRB
FLVEL
**
**HISTORYOUTPUT:
H-Output-1
**
*Output,history,variable=PRESELECT
*nodeoutput,nset=inifl
por
*elementoutput,elset=one
LEAKVRT
LEAKVRB
ALEAKVRT
ALEAKVRB
SDEG
*nodeprint,nset=inifl
por
*EndStep
**----------------------------------------------------------------
**
**STEP:
Hold
**
*Step,name=Hold,nlgeom=YES,unsymm=YES
hold
**MAKECHANGESHERE
**Usethe7200sdatalineinplaceofthe100sdataline
**Soils,consolidation,end=SS,utol=1.E6
**100.,100.,0.005,1500.,10.0
*Soils,consolidation,end=SS,utol=1.0e+6
100.,7200.,0.005,1500.,10.0
**
**BOUNDARYCONDITIONS
**
**Name:
BC-4Type:
Velocity/Angularvelocity
*Boundary,type=VELOCITY
_PickedSet273_1,1
_PickedSet273_2,2
_PickedSet273_3,3,3
**
*cflow
inifl,,0.0
**
**OUTPUTREQUESTS
**
*Restart,write,frequency=1
**
**
**FIELDOUTPUT:
F-Output-1
**
*Output,field
**MAKECHANGESHERE
**Removethenset=QA,andelset=QAparametersettings
**NodeOutput,nset=QA
*NodeOutput
POR,U
**elementoutput,elset=QA
*elementoutput
S,LE,PE
GFVR
SDEG
DMICRT
PFOPEN
LEAKVRT
LEAKVRB
ALEAKVRT
ALEAKVRB
FLVEL
**
**HISTORYOUTPUT:
H-Output-1
**
*Output,history,va
升级会员 