Page 156 - 《水资源与水工程学报》2024年第6期
P. 156
1 2 & ' ( ) & * + , - 2024 $
5
珋
2
珋
珋
2
s 1 s s A ω p μ s ø m =>(28), FÆ(26)`EF
2
k( + )+k =μ s p 珋 s-
z
r
2
2
2
2
z r r z ω +p í?@AßBCD1o, ¤ª(29)。
(14) 珌 s=C K(Nr)+C I(Nr)+ 珌 s (29)
01 0 0 02 0 0
珋 s(l,z,p)=0 (15) .:C 、C ªçs;N = 槡 p μ s /k; 珌 s 1
r
0
01
02
珋
s(r,0,p) =0 (16) oGª。 ¬fbw¶=Ó¤¥(30)~(32)。
z
A ω M
珋
s(r,M,p) 珌 s = 2 2 (30)
=0 (17) p +M
z
- 珌 sNI(Na)-mI(Nl)
0
0 1
(0<z ≤ d)
0
2 0
0
r { -Q 0 0 0 1 0 1 0 0 0
珋
s(a,z,p) C = N(I(Nl)K(Na)+I(Na)K(Nl))
01
=
2a(M -d) π kp (d<z<M) ( 31)
r
- 珌 sNK(Na)+mK(Nl)
(18) C = 0 1 0 2 0 0
02
0
1
0
0
0
1
0
0
0
L9(16)、(17)_` Fourier ¡ N(I(Nl)K(Na)+I(Na)K(Nl))
(32)
, ;
M 珋 scos( δ z)dz (n ≠ 0) ¬çs C ,C ·¸=Óª(23)、(29)
∫
01
02
{ 0 ., ³´ Fourier ¡¢b·R
珌 s= M (19)
∫ ,Ǥ¥ Laplace §¨KÃ
ª。
珋 sdz (n=0)
0 珌 s 2 ∞ A ω
珋
.: 珌 s s Fourier ¡:ß; δ= 珓 s(r,z,p)= M + M ∑ 珌 scos( δ z)+ ω +p
2
2
n=1
n π /M,n=1,2,3……(n ø¯2s)。 (33)
/ n ≠ 0 Æ, ³´(19) (14)、(15)、 FÆ Laplace §¨ªÎno, HؽÃäª。
( 18) _`¤¥(20)~(22)。 [23]
I£J4 uK·«z¹¤a,Stehfest s£
2
2
珌
珌
2s s ( δ k +μ s p) 珌 s
z
2
r 2 +r -r =0 (20) ¢}×äªIJKLM ëLM»R,
r r k EF³´ Stehfest s£¢}b Laplace ¢
r
珌 s(l,p)=0 (21)
·RQ[äªÇ¤¥(34)。
珌
s(a,p) Q sin( δ d) N
= (22) ln2
r 2a(M -d) π kp δ s(r,z,p)= t ∑ V 珓 s(r,z,p)-Asin( ω t)
i
r
ø m =>(22), FÆ(20)°*í? i= 1
1
( 34)
@AßBCD1o, j(23) 䤪。
.:p=i(ln2)/t,i ø1~N,N ¯Ns, Xòø
珌 s=CK(Nr)+CI(Nr) (23)
1 0 2 0
8~16;V Stehfest s£¢}äbOs。
i
.:C、C ªçs;K 0 ?@Aß- 2 ÙB
0
2
1
3.2 ABwvxvycJxz{
2
槡
CDs;N= (p μ s +k δ )/k;I0 ?@Aß / l → ∞、A=0、 ω =0 Æ, ©ª«ª®
r
0
z
- 1 ÙBCDs。 ³´@AßBCDsR®,
ÕðñKèÎ
'OZ
¬fbw¶=ӪǤa(24)~(25)。
[Kå`PQMª, FÆ#"BCDsR®
-mI(Nl)
C = 1 0 (24) I(Nl)、K(Nl) → ∞, (34) Ç®
1
N(I(Nl)K(Na)+I(Na)K(Nl)) 0 0
1
0
0
1
Q K(Nr) 2M
0
0
mK(NI) s= ( - ·
1 0
C = ( 25) 2 π kMpaNK(Na) π (M -d)
2 r 0 1 0
N(I(Nl)K(Na)+I(Na)K(Nl))
1
1
0
0
∞ K(Nr)sin( δ d)cos( δ z)
.:K 1 ?@Aß- 2 ÙBCDs;I 1 ? ∑ 0 0 ) (35)
1
1
@Aß- 1 ÙBCDs。 n=1 NaK(Na) δ
1
¬(35) P Stehfest s£¢:, l
/ n=0 Æ, ³´(19) _`¤¥
4 [20] ¥^'OZ[K`PQM©ª«ªÎ
珌
2
珌
2s s 2 μ s p 珌 s -A ω Mp μ s 2
r +r -r = r (26)
2
2
2
r r k k( ω +p) XQ。 c/
'O©ì a → 0 Æ, L9BCDs
r
r
珌 s(l,p)=0 (27) R®limNaK(Na)=1, (35) ÇR%®
i
i
i
a →0
珌
s(a,p) Q Q 2M
= (28) s= (K(Nr)- ·
0
0
r 2akp 2 π kMp π (M -d)
r
r