An isogeometric analysis based on a kind of local bivariate B-spline function and Burton-Miller method is introduced in boundary element method (BEM) for 3D Helmholtz problems.The method avoid efficiently computation of singular and nearly singular integrals aroud singular point of parametric surface.Numerical examples show that the algorithm has good performance in accuracy, efficiency and convergence.The algorithm performs isogeometric BEM (IBEM) with traditional global B-spline shape functions on CPU time and accuracy.The method has potential in engineering applications.