地下水位变化可否使用CEL法模拟？

我使用CEL法，建立了单桩贯入土体的模型，想考虑地下水这个因素，对沉桩阻力和挤土效应的影响，请问可以使用CEL法模拟吗？如何实现，有没有相关文章可阅读？