This paper is concerned with a three-dimensional time harmonic model of open shell structures buried in poroelastic soils. It combines the dual boundary element method (DBEM) for treating the soil and shell finite elements for modelling the structure, leading to a simple and efficient representation of buried open shell structures. A new fully regularised hypersingular boundary integral equation (HBIE) has been developed to this aim, which is then used to build the pair of dual BIEs necessary to formulate the DBEM for Biot poroelasticity. The new regularised HBIE is validated against a problem with analytical solution. The model is used in a wave diffraction problem in order to show its effectiveness. It offers excellent agreement for length to thickness ratios greater than 10, and relatively coarse meshes. The model is also applied to the calculation of impedances of bucket foundations. It is found that all impedances except the torsional one depend considerably on hydraulic conductivity within the typical frequency range of interest of offshore wind turbines.