為進(jìn)一步提高地下水開采誘發(fā)地面變形模擬的準(zhǔn)確性，以比奧固結(jié)理論為基礎(chǔ)，結(jié)合土體非線性流變理論，將比奧固結(jié)理論中的本構(gòu)關(guān)系推廣到黏彈塑性，并考慮孔隙度、滲透系數(shù)、土體的變形參數(shù)隨有效應(yīng)力的動(dòng)態(tài)變化關(guān)系，建立了地下水開采與地面沉降三維全耦合數(shù)學(xué)模型，并應(yīng)用伽遼金有限單元法對(duì)模型進(jìn)行求解，模擬了河北省滄州市地下水開采引起的地下水滲流場(chǎng)、地面沉降、地裂縫發(fā)生發(fā)展趨勢(shì)。結(jié)果表明：抽水引起地下水位降落漏斗和地面沉降的同時(shí)，將導(dǎo)致土體發(fā)生總體指向地下水位降落漏斗中心的水平位移，水平位移較大的地區(qū)集中分布在地面沉降漏斗的邊緣，是地裂縫的高發(fā)區(qū)。