Abstract: Most conventional analytical methods used in reservoir pressure behavior analysis are based on regular geometry reservoirs.But the actual carbonate rock reservoirs are almost irregular with different size fractures and vug，which makes pure analytical method helpless.Boundary element method (BEM) is a boundary type numerical method，it changes governing differential equations to integral equations of boundaries.In this paper，BEM is established on the basis of single phase，homogeneous and compressible fluid flowing mathematical model by Laplace alteration.Several pressure behaviors are analyzed with different fractures and vug array mode by associating with curve charts worked out with this method.It offers technical support for resolving reservoir parameters，deliverability evaluation and more flexible theory method for well test analysis.

Key words: beads shaped reservoir, numerical method, pressure dynamic, boundary element method, well test analysis